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VERIFICATION  AND  VALIDATION  FINAL  REPORT 
FOR  FIRE  AND  SMOKE  SPREAD  MODELING  AND  SIMULATION  SUPPORT 

OF  T-AKE  TEST  AND  EVALUATION 


1.0  PURPOSE 


The  Operational  Requirements  Document  (ORD)  for  the  Auxiliary  Dry  Cargo  Carrier, 
T-AKE,  designates  survivability  as  one  of  the  key  performance  parameters  for  the  ship.  The 
objective  requirements  linked  to  this  parameter  are  multi-fold,  and  include  the  following: 

•  Built-in  redimdancy  so  that  a  single  equipment  failure  or  loss  will  not  degrade  system 
capability.  Redundant  system  functions  will  be  provided  and  will  be  separated  so  as  to 
avoid  the  possibility  of  total  system  loss  due  to  a  single  casualty. 

•  Low  expansion  fire  fighting  foam  for  fire  protection  of  the  helicopter  deck  and  hangar 
and  machinery  space  bilges. 

•  Water-mist  or  other  USCG  approved  total  flooding  fire  suppression  system  for  machinery 
spaces. 

•  Seawater  fire  sprinkler  systems  (wet  firemain/dry  sprinklers)  for  cargo  ordnance  holds. 

These  requirements  are  carried  into  the  System  Specification,  which  states  that  the  ship  shall 
meet  the  above  firefighting  requirements  under  all  loading  conditions.  It  is  the  intention  of  the 
T-AKE  Program  to  use  modeling  and  simulation  (M&S)  to  certify  that  the  ship  meets  these 
requirements.  The  specific  methodology  (tool)  that  will  be  used  in  this  effort  consists  of  a 
combination  of  computer  fire  modeling  using  the  “Consolidated  Model  of  Fire  Growth  and 
Smoke  Transport  (CFAST),  Version  3.1.7”  (Peacock,  et  al.,  January  2000a;  Peacock,  et  al., 
January  2000b)  and  prescriptive  rules  for  fire  spread  determined  with  the  finite  difference  heat 
conduction  model  -HEATING  Version  7.3  (Back,  et  al.,  March  2003).  The  objective  of  this 
document  is  to  demonstrate  that  the  methodology,  including  supporting  models,  provide 
reasonable  and  conservative  results  in  the  context  of  using  for  T-AKE  Vulnerability  Assessment 
Report. 


2.0  TOOLS  SUMMARY 

2.1  Tool/Model:  Fire  Spread  Methodology 


Describing  the  methodology  for  predicting  fire  and  smoke  spread  requires  an  understanding 
of  the  shipboard  damage  zones.  The  zones  are  categorized  in  three  ways: 

Primary  Damage  Area  (PDA(F))  compartments:  PDA(F)  compartments  are  those  that  incur 
significant  damage  from  the  weapon  hit  (e.g.,  failed  deck/bulkhead  and  shock  holing).  From  the 
standpoint  of  fire  and  smoke  spread,  the  vent-connected  spaces  are  included.  This  is  due  to  the 
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fact  that  the  vents  provide  a  means  for  spreading  fire  and  smoke  directly  as  well  as  supplying 
additional  oxygen  for  the  fires. 

Adjacent  to  Primary  Damage  Area  (APDA(F))  compartments:  These  compartments  share  a 
common  boundary  (bulkhead,  overhead  or  deck)  with  a  PDA(F)  compartment.  The 
recommended  definition  for  fire  and  smoke  spread  for  adjacent  compartments,  APDA(F), 
excludes  any  compartments  that  have  vent  openings  to  the  PDA(F). 

Beyond  Adjacent  to  Primary  Damage  Area  (BAPDACFl)  compartments:  These 
compartments  border  APDA(F)  spaces.  They  can  be  otherwise  described  as  “twice  removed” 
from  the  PDA(F).  These  spaces  may  be  a  concern  due  to  potential  fire  spread  later  into  the 
event.  It  is  typically  assumed  that  there  are  no  communicable  openings  (i.e.,  vents,  ducts,  open 
doors  or  open  hatches)  between  BAPDA(F)  and  APDA(F)  spaces.  This  is  based  on  the 
assumption  that  the  ship  is  at  general  quarters  prior  to  the  weapon  detonation. 

The  Methodology  for  predicting  fire  and  smoke  spread  in  a  naval  platform  follows  the  steps 
listed  below.  Steps  (3)  -  (5)  are  described  in  further  detail  in  the  Appendices,  as  they  constitute 
the  phenomenology  that  is  included  in  the  tools  being  verified  and  validated. 

1 .  Establish  the  primary  damage; 

2.  Review  the  primary  damage  information; 

3.  Characterize  the  secondary  damage  (fire  and  smoke  spread)  in  the  PDA(F); 

4.  Characterize  the  secondary  damage  to  PDA(F)  spaces; 

5.  Characterize  the  secondary  damage  to  APDA(F)  and  BAPDA(F) 
compartments;  and 

6.  Document  the  time  dependent  secondary  damage  estimates  in  snapshot  format. 

Steps  3  and  4  involve  characterizing  the  secondary  damage  in  the  PDA(F)  by  two  different 
approaches — CFAST  fire  model  and  prescriptive  rules.  For  conservative  estimates,  the  predicted 
tenability  and  fire/smoke  spread  are  dictated  by  the  approach  that  estimates  the  shorter  time  scale 
(i.e.,  more  conservative)  for  progressive  damage. 


2.2  Overview 


The  current  approach  for  the  characterization  of  secondary  damage  associated  with  a  weapon 
hit  has  evolved  over  several  Total  Ship  Survivability  (TSS)  and  Battle  Damage  Repair 
Assessment  (BDRA)  evaluations.  The  genesis  of  the  current  Methodology  was  rudimentary  fire 
spread  times,  developed  by  the  Naval  Research  Laboratory  (NRL)  and  Hughes  Associates,  Inc. 
based  on  USS  Stark  incident  data  and  testing  (Scheffey,  et  al.,  April  1992).  This  information 
was  implemented  in  a  TSS  model  developed  by  the  Naval  Surface  Warfare  Center,  Carderock 
Division  (NSWC/CD).  The  fire  spread  methodology  was  subsequently  improved  and  the 
framework  was  expanded  to  allow  for  smoke  spread  during  the  evaluation  performed  for  the 
LPD-17  VAR  2  (White,  et  al.,  September  1997).  Most  recently,  the  methodology  was  further 
refined  during  the  CVNX  assessment  (Back,  et  al.,  April  2002).  This  report  reflects  the  most 
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current  methodology  used  for  the  CVNX  assessment  and  for  the  VAR  3  proeess  for  the  LPD-17 
platform  detailed  design.  This  approach  is  proposed  to  be  adopted  for  T-AKE. 


2.3  Intended  Purpose 


The  vulnerability  assessment  process  that  uses  this  methodology  includes  the  simulation  of 
weapon  damage,  characterization  of  secondary  damage  (e.g.,  fire  and  smoke),  and  description  of 
the  battle  damage  repair  activities  and  effectiveness.  These  activities  are  generally  referred  to  as 
a  TSS  and  Battle  Damage  Repair  Assessment  BDRA.  The  overall  objective  of  the  TSS/BDRA  is 
to  identify  specific  weaknesses  and/or  shortcomings  as  well  as  high  payoff  research  and 
development  efforts  that  might  inerease  the  survivability  of  the  particular  ship  platform.  The 
current  TSS/BDRA  process  involves  consideration  of  various  weapon  hits.  For  each  specific 
weapon  threat  scenario  (shot  line),  the  initial  damage  is  identified  by  utilizing  a  blast  damage 
model.  Secondary  damage,  such  as  fire  and  smoke  spread  that  results  from  the  blast  damage, 
must  also  be  characterized.  The  Methodology  described  in  this  document  concentrates  on  the 
fire  and  smoke  spread  methodology  only.  A  framework  is  outlined  by  which  a  conservative 
estimate  of  the  secondary  fire  and  smoke  damage  can  be  determined,  supported  by  supplemental 
fire  test  data  from  current  and  previous  tests. 


2.4  Tool/Model  Accediting  Organization. 

PEO  Ships,  PMS  325B  (TAKE  Program  Office)  80  M  Street,  Suite  610  Washington, 
DC  20003. 

2.5  CFAST 

2.5.1  Required  Operating  Environment 

CFAST  will  run  on  a  386  or  later  IBM  compatible  PC  with  at  least  4  MB  fi-ee  extended 
memory,  at  least  VGA  compatible  video  card  and  a  100%  MS-MS  compatible  mouse. 

2.5.2  Existing/Published  Documentation 

An  extensive  list  of  references  is  provided  at  the  end  of  this  document  (Section  9.0  - 
References)  and  is  reviewed  in  detail  in  Appendix  A. 
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2.6  Heuristics  and  HEATING 
2.6. 1  Required  Operating  Environment 


The  minimal  environment  for  running  HEATING  is  a  computer  system  with  a  FORTRAN  90 
compiler  and  a  text  editor.  This  suffices  to  create  the  executables  and  the  input  files.  HEATING 
has  two  main  formats  for  data  output:  an  ASCII  file  for  results  at  user-specified  nodes  and  a 
binary  file  for  field  data.  The  most  expedient  environment  for  post-processing  the  subsequent 
results  would  be  any  system  that  runs  Tecplot,  a  commercial  scientific  plotting  tool  sold  by 
Amtec,  Inc.  However,  the  h7tec  FORTRAN  code  that  comes  with  HEATING  can  be  used  on 
any  system  with  a  FORTRAN  90  compiler,  a  text  editor,  and  a  scientific  plotter  to  visualize  the 
results.  The  process  entails  employing  h7tec  to  extract  data  fi'om  HEATING'S  binary  data  file 
and  then  editing  the  resulting  ASCII  file(s)  to  match  a  specific  graphing  program's  requirements. 
HEATING  7.3  comes  with  Windows-compatible  executables  of  the  main  simulation  plus  all  the 
supporting  programs  (such  as  h7tec).  Again,  the  resultant  data  is  best  plotted  with  Tecplot  but  it 
is  also  possible  to  use  h7tec  to  extract  the  field  data  and  Wordpad  to  edit  it  so  that  another 
Windows-based  plotter  could  be  utilized.  Node  data  can  be  easily  plotted  with  Excel.  The 
operating  environments  used  for  the  present  verification  are  detailed  in  Appendix  C.  The 
heuristics  derived  from  test  data  do  not  require  a  model  shell  or  computer  fi*amework  per  se. 


2.6.2  Existing/Published  Documentation 


An  extensive  list  of  references  is  provided  at  the  end  of  this  document  (Section  9.0  - 
References).  The  derivation  of  the  fire  spread  rules  (heuristics)  is  described  in  Section  4.0  and 
Appendix  B.  A  review  of  HEATING  is  provided  in  Appendices  C  and  D. 


2.7  Take  Ship  Data  Set 


The  TAKE  General  Arrangement  and  Ventilation  drawings  will  be  used  for  the  Ship  data  set 
when  the  accreditation  for  the  Vulnerability  Assessment  is  accomplished.  For  this  Validation 
experimental  data  from  the  Ex-USS  Shadwell  was  used  to  validate  volumetric  data. 


3.0  RISK  ANALYSIS  AND  ASSESSMENT  SUMMARY 
3.1  Methodology 


The  Methodology  described  in  this  report  outlines  the  fi-amework  through  which  a 
conservative  estimate  of  the  secondary  fire  and  smoke  damage  can  be  established.  Fire  may 
spread  to  compartments  within  the  PDA(F)  via  either  hot  gases  passing  through  an  open  vent  or 
heating  of  an  intact  boundary  (bulkhead  or  deck).  In  this  situation,  fires  resulting  from  the 
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weapon  detonation  in  the  immediate  vicinity  of  the  burst  point  would  spread  to  other 
compartments  within  the  PDA(F).  Both  modes  of  fire  spread  are  examined  when  evaluating  fire 
spread  to  these  PDA(F)  compartments  to  determine  which  mode  results  in  a  faster  (more 
conservative)  fire  spread  time. 


3.2  CFAST 


In  the  first  mode  of  spread,  the  temperature(s)  of  hot  gases  are  determined  using  CFAST. 
Based  on  the  review  of  several  potential  ignition  scenarios,  it  is  conservatively  assumed  that 
ignition  via  hot  gases  will  occur  in  a  PDA(F)  compartment  when  the  accumulated  fire  gas 
temperatures  reach  230®C  (446°F)  in  that  adjacent  space  (Back,  et  al.,  February  2000;  Naval 
Ships’  Technical  Manual,  Dec  2001;  Setchkin,  1949;  Darwin,  et  al.,  1994;  Drysdale,  1999). 

CFAST  was  designed  to  execute  very  quickly  (in  real-time,  or  quicker)  on  affordable  desktop 
computers.  In  addition,  it  was  intended  to  be  relatively  simple  to  use,  which  makes  the  model 
accessible  to  users  (such  as  ship  architects  and  shipboard  damage  control  personnel)  who  are  not 
experienced  fire  modelers.  In  order  to  meet  these  requirements,  the  CFAST  developers  made 
certain  simplifying  assumptions  (in  addition  to  those  which  are  inherent  in  zone  models) 
regarding  the  range  of  phenomena  and  the  types  of  scenarios  that  could  be  modeled.  In  turn, 
these  assumptions  led  to  limitations  on  the  scope  of  the  problems  which  can  be  modeled,  the  ease 
with  which  they  can  be  modeled  and  the  accuracy  of  the  results.  These  limitations  can  be 
characterized  as  inherent  (things  which  CFAST  is  not  capable  of  doing),  or  practical  (things 
which  require  a  priori  knowledge  which  is  not  likely  to  be  available  or  which  require  too  much 
time  and  effort  to  be  feasible).  Inherent  limitations  include  inability  of  objects  to  self-heat  based 
on  heat  transfer  and  spontaneously  ignite.  Practical  limitations  include  inability  to  predict 
phenomena  such  as  time-dependent  solid  phase  combustion  of  Class  A  fuels. 

The  limitations  in  Table  1  were  identified  in  Hoover  et  al.,  June  2000,  using  CFAST 
Version  3.1.4  (with  improved  radiation,  vertical  vent  flow  and  vertical  heat  transport  algorithms). 
One  of  the  most  restrictive  limitations  of  CFAST  is  the  requirement  that  “intelligent”  input 
parameters  be  used  to  run  the  model.  This  imposes  restrictions  on  what  can  be  modeled  and 
what  approximations  must  sometimes  be  made  to  circumvent  those  limits. 


3.3  HEATING 


In  the  second  mode  of  evaluation  for  fire  spread,  the  prescriptive  rules  used  in  conjunction 
with  the  finite  element  model  HEATING  are  not  plagued  with  the  numerous  limitations  inherent 
in  computer  models.  Hence,  this  part  of  the  analysis  provides  for  conservative  estimates  of  fire 
spread  to  APDA(F)  spaces  and  BAPDA(F)  spaces,  which  moderates  the  variability  that  can 
occur  in  the  results  from  the  use  of  computer  modeling  only  (particularly  for  CFAST).  The  use 
of  these  prescriptive  rules  is  conservative  since  it  is  assumed  that  there  is  an  unlimited  amount  of 
fuel  and  air  available  in  each  compartment.  In  the  event  that  there  is  a  small  amount  of  fuel  in 
the  compartment  or  the  compartment  is  closed,  fire  spread  may  be  delayed  or  may  not  occiu  at 
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all.  The  combination  of  the  computer  fire  model  and  prescriptive  rules  provides  an  inherent 
safety  factor  in  the  Methodology. 

Table  1  -  Identified  Limitations  Of  CFAST 


Limitation 

Only  one  internal  ambient  state  is  permitted. 


Compartments  must  be  rectilinear. 


Ceiling,  floor  and  walls  are  each  limited  to  a 
single  set  of  thermo-physical  properties. 

Compartments  have  only  one  wrap-around 
wall,  not  individual  surfaces. 


User-specified  fire  histories  are  required. 


CFAST  does  not  account  for  the  variability  of 
fire  chemistry  with  pyrolysis  and  combustion 
conditions. 


Effects 

May  lead  to  minor  errors  in  predictions  if  there 
are  large  initial  temperature  differences  among 
the  compartments. 

Heat  transfer  via  conduction  and  radiation  will 
not  be  correct  for  surfaces  which  are  not 
rectangular.  Approximation  methods  that 
minimize  this  limitation  are  provided  in  Hoover 
et  al.,  2000. 

Approximations  must  be  made  if  the  ceiling, 
floor  or  walls  are  composed  of  multiple  regions 
having  different  properties. 

May  lead  to  significant  errors  if  there  are  large 
differences  in  the  properties  of  walls.  Work¬ 
around  methods  have  been  devised  in  Hoover  et 
al.,  2000. 

The  user  must  know  (or  be  able  to  accurately 
estimate),  in  advance,  the  development  of  the 
fire. 

Good  subject  matter  experts  should  estimate 
user-specified  heat  release  rates  so  reasonable 
pyrolysis  and  combustion  product 
concentrations  will  be  predicted. 


4.0  SUMMARY  OF  PREVIOUS  DEVELOPMENTS  AND  APPLICATIONS 


The  Methodology  described  in  this  report  has  been  successfully  used  in  other  Navy  damage 
control/fire  protection  applications.  Additionally,  elements  of  the  Methodology,  in  particular 
CFAST,  have  been  used  to  analyze  unique  Navy  situations.  Specifically,  the  techniques  have 
been  used  for  a  broad  range  of  problems  encountered  in  Fire  Hazard  Analyses  (FHAs).  This 
section  briefly  describes  these  applications. 
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4.1  Methodology  Use  in  Platform  Vulnerability  Assessment  Reports  (VARs) 


The  genesis  of  the  rules-based  aspect  of  the  Methodology  was  the  USS  STARK  fire 
dynamics  testing  at  the  NRL  Chesapeake  Beach  Detachment  (Back  et  al.,  1991)  and  ex-USS 
SHAD  WELL  (Scheffey  et  al.,  1999)  fire  test  facilities.  The  heat  transfer  data  from  these  tests 
was  summarized  and  published  in  the  scientific  literature  (Darwin  et  al.,  1994).  Developers  of 
the  Ship  Vulnerability  Model  (SVM)  at  NSWCCD  desired  a  relatively  simple  method  to 
integrate  fire  and  smoke  spread  characteristics  into  their  battle  damage  model.  The  STARK  test 
data  were  combined  with  hand  calculation  techniques  to  develop  general  rules  for  fire  and  smoke 
spread  for  large,  post-flashover  fires  on  Navy  ships.  These  rules,  published  in  1992  (Scheffey  et 
al.,  1992)  were  the  precursor  to  the  Methodology  described  in  this  report.  These  heuristics 
continue  to  be  adequate  for  many  applications,  particularly  where  there  is  significant  variability 
in  potential  fire  growth/spread  scenarios. 


4.1.1  LPD17VAR 


The  Methodology  was  first  developed  and  applied  in  the  LPD  17  VAR  3  (bid  design).  It  was 
recognized  that  a  strictly  rules-based  approach  could  lead  to  over-conservative  predictions 
(i.e.,  very  rapid  fire  and  smoke  spread).  Modeling  of  the  near-field  area  of  the  Primary  Damage 
Area  (PDA)  was  included  to  assess  the  fire  spread  in  vent-coimected  areas  in  the  PDA.  Fire 
spread  heuristics  and  a  library  of  space-specific  fire  growth  curves  were  then  used  in  ADPA 
spaces.  The  LPD  17  VAR  3  fire  and  smoke  spread  analysis  was  subsequently  documented 
(White  et  al.,  September  1997)  and  approved  by  Office  of  the  Secretary  of  Defense,  Live  Fire 
Test  and  Evaluation  (OSD  LFT&E).  The  Methodology  is  currently  being  used  in  the  LPD  17 
VAR  4  (detailed  design). 


4.1.2  CVN(X)VAR 


A  revised  Methodology  was  applied  to  the  CVN(X)  VAR  2.  The  primary  improvement  was 
a  modified  fire  growth  algorithm  for  very  large  spaces.  A  criticism  of  the  previous  Methodology 
approach  was  the  very  rapid  fire  development  in  all  spaces,  independent  of  size.  It  was 
perceived  that  fire  development  (i.e.,  time  from  ignition  to  flashover)  would  be  longer  in  larger 
spaces  compared  to  small  spaces. 

The  Methodology  as  applied  to  CVN(X)  (now  designated  CVN  21)  is  documented  (Back  et 
al.,  April  2002);  the  VAR  is  in  the  final  approval  stages.  It  is  anticipated  that  the  Methodology 
will  be  used  for  the  CVN  21  VAR  3  now  being  initiated. 

A  non-classified  version  of  the  Methodology  was  published  for  general  use  (Back  et  al., 
March  2003).  This  document  forms  the  primary  technical  basis  of  this  report  as  described  in 
Appendix  B.  It  was  during  the  CVNX  VAR  2  and  subsequent  documentation  of  the  general 
Methodology  that  the  PDA(F)  and  ADPA(F)  definitions  were  adopted  which  clarify  the 
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differences  in  the  treatment  of  these  areas  by  the  blast  damage  and  fire  communities  (See 
Appendix  B). 


4.2  Use  in  WET  Predictions 

4.2.1  Ex-USS  CARON 


The  basic  Methodology  was  used  to  predict  unabated  fire  and  smoke  spread  in  weapons 
effects  testing  (WET)  of  the  ex-USS  CARON  sponsored  by  PMS  500.  The  results  are 
documented  in  stand-alone  segments  which  were  incorporated  in  the  ex-USS  CARON  WET  pre¬ 
hit  and  post-hit  reports  (Back  et  al.,  January  2002,  Back  et  al,  July  2003). 


4.2.2  Ex-USS  PETERSON 


The  basic  Methodology  is  being  used  in  pre-hit  predictions  for  the  DD(X)  industry-sponsored 
ex-USS  PETERSON  WET.  A  fire  suppression  algorithm  has  been  added  to  the  Methodology 
since  a  water  mist  system  designed  to  prevent  fire  spread  will  be  installed  on  the  ship.  The  water 
mist  suppression  algorithm  is  based  on  a  steady-state  suppression  model  documented  in  the  peer- 
reviewed  scientific  literature  (Back  et  al.,  November  2000).  Documentation  of  this  prediction  is 
in  preparation.  The  ex-USS  PETERSON  will  also  be  assessed  using  a  network  fire  and  smoke 
spread  model.  This  model  is  being  adopted  by  the  DD(X)  Recoverability  Team  and  will  be  used 
as  a  replacement/improvement  to  the  Methodology  for  DD(X)  VAR  1.  The  ex-USS  PETERSON 
test  data  will  be  used  as  part  of  the  V&V  of  the  network  model.  Floyd  et  al.,  2003,  briefly 
describes  this  new  modeling  approach. 


4.3  CFAST  Modeling  for  Navy  Shipboard  FHAs 


CFAST  has  been  used  in  a  number  of  Navy  ship  platforms  to  model  unique  or  hazardous 
situations.  This  modeling  usually  is  incorporated  into  a  Fire  Hazard  Analysis.  An  FHA  is  used 
to:  inventory  the  fuel  loading;  assess  ignition  sources;  quantify  the  hazard  resulting  from  likely 
fire  scenarios;  determine  the  impact  on  personnel,  equipment,  and  the  ship;  and  develop 
appropriate  mitigation  strategies.  CFAST  has  been  an  important  tool  for  quantifying  fire  hazards 
on  ships.  Table  2  outlines  Navy  shipboard  FHAs  where  CFAST  has  been  used  to  quantify  the 
hazard.  The  wide-ranging  application  of  CFAST  in  these  situations  demonstrates  the  breath  and 
scope  of  the  use  of  this  model  for  Navy  fire  hazard  assessment. 


5.0  COMMENTS  AND  TROUBLE  REPORTS 


The  methodology  documented  in  this  report  reflects  the  latest  approach  to  characterize  fire 
and  smoke  spread  in  the  VAR  3  process  for  the  LPD  17  platform  detailed  design  and 
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CVNX  VAR.  Comments  and  trouble  reports  received  from  users  and  technologists  relying  on 
the  methodology  results  have  been  limited  to  the  conservative  nature  of  the  approach.  This 
criticism  has  been  addressed  by  the  development  of  a  network  fire  model  (Floyd  et  al.,  2003). 
The  network  fire  model  would  provide  a  physics-based  model  that  could  quickly  assess  fire  and 
smoke  spread  for  ship  platforms.  Plans  also  include  accounting  for  fire  detection  and 
suppression  effects.  This  type  of  model  would  provide  better  far  field  secondary  damage 
characterization  capabilities.  It  is  unclear  if  the  network  fire  model  will  replace  the  current 
methodology  or  if  it  will  be  used  in  conjunction  with  the  methodology  or  portions  of  the 
methodology  described  in  this  report. 

Table  2  -  Use  of  CFAST  in  Navy  Shipboard  Fire  Hazard  and  Modeling  Applications 


Client/Platform 

Task 

Reference 

LPD  17 

Analysis  of  smoke  and  heat 
spread  for  well  deck  fire 
fighting 

White  et  al.,  2001 

ONR/DC-ARM 

Modeling  of  the  probability  of 
compartment  flashover 
resulting  from  a  missile  hit 
scenario  (geometry  and  vent 
dependent) 

Back  et  al.,  October  1998 

NSWC/OHIO  and  VIRGINIA 
Class  Submarines 

Fire  modeling/risk  assessment 
of  plastic  waste 
storage/disposal 

Back  et  al.,  December  2000 

LPD  17 

Composite  mast  fire  hazard 
analysis 

White  et  al.,  2000 

NSWC/DD 

USS  RADFORD  AEMS 
composite  mast  fire  hazard 
analysis 

White  et  al.,  1996 

PMS  400/DDG 

Fire  hazard  analysis  of  pipe 
insulation 

Beitel  et  al.,  1996 

PMS  378/C  VN  21 

Fire  hazard  analysis  of  hangar 
bay 

Parker  et  al.,  2003 

Submerged  vessel 

Fire  hazard  analysis  of  dry 
deck  shelter 

White  et  al.,  November  1997 

6.0  ALGORITHM  REVIEW 
6.1  Methodology 


The  Methodology  is  summarized  in  the  flowchart  in  Figure  1.  The  process  involves 
modeling  PDA(F)  compartments  with  the  multi-compartment  zone  fire  model  CFAST.  The  use 
of  CFAST  allows  for  a  more  complete  characterization  of  these  compartments.  The  output  from 
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CFAST  is  used  to  predict  fire  spread  times.  Fire  spread  beyond  these  compartments  to  APDA(F) 
and  BAPDA(F)  compartments  is  also  predicted  using  a  conservative  set  of  prescriptive  rules. 


6.2 

CFAST 

6.2.1 

History 

HAZARD  1  Version  1.0  was  the  prototype  fire  hazard  assessment  model  developed  by  the 
National  Institutes  of  Standards  and  Technology  (NIST)  that  integrated  fire  dynamics  with 
biology,  toxicology  and  human  behavior  (Jones,  1985;  Jones,  1983).  FAST  version  18.3  was 
implemented  in  HAZARD  1.0  in  1988  and  provided  the  user  with  a  method  to  simulate  the 
important  time-dependent  phenomena  involved  in  residential  and  commercial  building  fires.  In 
1990  FAST  Version  18.5  was  integrated  with  Consolidated  Compartment  Fire  Model  (CCFM)  to 
develop  CFAST  1.0  that  was  functionally  equivalent  to  FAST  but  more  modular  like  CCFM. 

The  current  versions  of  CFAST  available  on  the  web  are  Version  5.0  (5/7/03)  and  Version  3.1.7 
(1 1/8/01)  that  use  the  Graphical  User  Interface  (GUI) — FAST. 

In  the  case  of  buildings,  the  assumptions  made  by  CFAST  have  proven  to  be  reasonably  good 
and  the  limitations  have  not  been  unduly  restrictive  (See  Table  1).  However,  naval  fire 
protection  problems  differ  in  many  critical  respects  from  those  found  in  the  civilian  world  and  it 
was  not  clear  how  the  limits  of  CFAST  may  impact  the  utility  of  the  model  for  naval 
applications. 

The  US  Navy  became  involved  in  development  of  CFAST  under  its  Office  of  Naval 
Research  (ONR)  program  conducted  through  the  Science  and  Technology  (S&T)  Damage 
Control  (DC)  Project  at  NRL.  The  emphasis  in  this  program  was  to  evaluate  and,  where 
necessary,  improve  the  capabilities  of  CFAST  in  predicting  fire  and  smoke  transport  in  naval 
platforms.  Specifically,  the  Navy’s  funding  has  been  directed  toward  making  CFAST  more 
useful  for  simulation  of  shipboard  fires  by  adding  capabilities  for  modeling  phenomena  that  are 
absent  from,  or  of  little  significance  to,  building  fires.  In  particular,  these  phenomena  have 
included  mass  transport  via  vertical  gas  flow  through  horizontal  vents  (represented  by  hatches 
and  scuttles)(Bailey,  1994);  energy  transport  via  conduction  through  decks  (vertical  heat 
transfer)  (Bailey,  et  al.,  1998;  Bailey  and  Tatem,  1995);  an  improved  radiation  transport  sub¬ 
model  (Hoover  et  al.,  1996;  Hoover  et  al.,  1997;  Hoover,  1995;  Hoover  and  Bailey,  2000);  and, 
algorithms  for  heat  conduction  through  bulkheads  as  well  as  smoke  movement  along  a  long 
passageway  (Bailey  et  al.,  2002;  Bailey  and  Tatem,  1998). 

Although  initial  Joint  efforts  between  NIST  and  NRL  were  targeted  toward  improved 
algorithm  development,  later  NRL  work  focused  on  how  CFAST’s  limits  may  impact  the  utility 
of  the  model  for  naval  applications.  That  particular  research  has  identified  CFAST’s  limitations 
determined  their  effects  on  shipboard  fire  modeling  and  developed  methods  for  circumventing 
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those  problems  that  have  significant  adverse  effects  (Hoover  and  Tatem,  June  2000;  Hoover  and 
Tatem,  December  2000)  (See  Table  1). 


6.2.2  Validation  of  New  Phenomena 


Vertical  Vent  Flow  -  was  developed  to  handle  vertical  gas  flow  through  horizontal  vents 
(i.e.,  scuttles  and  hatches)  as  opposed  to  doors  and  windows.  This  phenomenon  has  importance 
in  naval  platforms  but  not  necessarily  in  buildings. 

Vertical  (and  Horizontal)  Heat  Transport  -  was  developed  to  correctly  account  for 
conductive  heat  transfer  in  the  vertical  (and  horizontal)  direction  through  metal  decks  in  naval 
platforms.  The  addition  of  this  phenomenon  in  shipboard  environments  is  necessary  since  rising 
temperatures  in  adjacent  steel  compartments  as  a  result  of  fire  can  cause  materials  to  reach 
ignition  levels  even  though  no  breach  has  occurred  in  the  initial  fire  compartment. 

Single  Zone  (SHAFT)  -  When  this  option  is  invoked,  CFAST  uses  a  single  set  of 
parameters  (i.e.,  temperature  and  gas  concentrations)  to  represent  the  entire  compartment 
volume,  rather  than  two  sets  (one  for  each  layer). 

Natural  and  Forced  Ventilation  -  was  developed  to  simulate  the  presence  of  heating, 
ventilation  and  air  conditioning  (HVAC)  systems  and  their  use  to  control  smoke  movement. 

Both  of  these  concepts  are  important  phenomena  that  must  be  accounted  for  in  naval  platforms. 

Radiation  Transport  -  was  developed  to  permit  the  use  of  growing  combustion  gas 
concentrations  to  calculate  the  absorbance  or  radiative  effect  of  the  fire  environment.  In  the  two- 
zone  CFAST  model  the  absorbance  of  each  layer  is  then  calculated  as  a  function  of  carbon 
dioxide,  water  vapor  and  soot. 

Corridor  Flow  —  was  developed  to  predict  the  flow  of  a  ceiling  jet  (hot  or  smoke  layer) 
down  the  corridor  with  time.  Prior  to  the  addition  of  this  phenomenon,  CFAST  assumed  an 
isothermal  upper  layer  would  immediately  cover  the  entire  ceiling  area  along  the  entire  length  of 
the  corridor  and  the  only  transient  change  in  the  corridor  was  the  smoke  layer  depth. 


6.2.3  CFAST  as  Fire  Model  for  PDA(F) 


CFAST  is  the  multi-compartment  zone  fire  model  used  to  characterize  the  environment  in 
PDA(F)  spaces  for  each  hit  scenario.  In  developing  the  input  for  CFAST,  the  ventilation 
conditions  must  be  defined.  These  conditions  include  those  provided  by  natural  and  forced 
ventilation.  Natural  ventilation  may  be  provided  by  holes  in  the  skin  of  the  ship  (resulting  from 
the  weapon  hit)  and  openings  between  PDA(F)  compartments.  At  a  minimum,  oxygen  for 
combustion  is  provided  naturally  through  any  openings  to  weather  that  result  from  the  weapon 
hit.  Forced  ventilation  may  be  provided  in  vital  spaces  where  ventilation  is  not  secured  under 
materiel  condition  ZEBRA. 
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For  the  initial  CFAST  run,  only  the  fuel  load  in  PDA(F)  compartments  is  considered.  It  is 
assumed  that  fires  ignite  in  all  PDA(F)  compartments  that  are  intimate  with  the  compartment 
containing  the  burst  point  and  have  available  fuel  at  the  time  of  the  weapon  detonation.  This  is  a 
very  aggressive  assumption.  Since  specific  details  about  the  fuel  loading  are  not  usually 
available,  fire  growth  rates  associated  with  the  type  of  materials  expected  in  the  spaces  are  used 
for  these  typical  materials.  More  specifically,  it  is  assumed  that  the  fires  follow  “t-squared”  fire 
growth  curves. 


6.3  HEATING 


The  amount  of  time  required  to  spread  fire  through  an  intact  boundary  is  determined  by 
identifying  the  fire  growth  time  and  the  thermal  penetration  time.  The  fire  growth  time  is  a 
function  of  type  of  compartment  (PDA(F)  as  compared  to  APDA(F))  or  BAPDA(F)  and/or 
compartment  size.  The  thermal  penetration  time  is  a  function  of  the  bulkhead  or  deck  thickness. 
The  thermal  penetration  time  is  determined  using  the  finite  difference  heat  conduction  model  - 
HEATING  (Version  7.3)  (Wickstrom,  1999),  which  has  been  previously  Verified  and  Validated 
(Trelles  et  al,  NRL  Report  in  prepziration;  Himt,  et  al.,  NRL  Report  in  preparation).  Fire  spread 
from  APDA(F)  spaces  and  BAPDA(F)  spaces  occurs  via  heat  transfer  through  intact  bulkheads 
and  decks.  This  mode  of  fire  spread  is  the  most  predominant  since  these  spaces  normally  do  not 
have  vent  cormections  to  the  spaces  with  fires. 


7.0  VERIFICATION  AND  VALIDATION  DESCRIPTION  AND  RESULTS 

7.1  Verification 

7.1.1  Methodology 


As  previously  stated,  the  Methodology  is  a  combination  of  CFAST  and  prescriptive  rules  that 
invoke  HEATING.  Its  verification  is  a  composite  of  the  verifications  of  these  two  components 
that  are  described  briefly  below. 

7.1.2  CFAST 


CFAST  is  based  on  solving  a  set  of  equations  that  predict  state  variables  (pressure, 
temperature  and  others)  based  on  the  enthalpy  and  mass  flux  over  small  increments  of  time. 
These  equations  are  derived  from  the  conservation  equations  for  energy,  mass,  and  momentum, 
and  the  ideal  gas  law.  The  conservation  equations  are  always  correct,  everywhere.  Thus  any 
errors  made  by  the  model  cannot  come  from  these  equations,  but  rather  come  from  simplifying 
assumptions  or  from  processes  having  been  omitted.  Examples  of  sources  of  errors  and 
omissions  have  been  highlighted  in  Section  3.0,  Table  1. 
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7.1.3  HEATING 


The  thermal  penetration  time  (tpeneiration)  was  determined  using  the  finite  difference  heat 
conduction  model  —  HEATING  Version  7.3.  In  the  verification  process  for  “HEATING” 
(Appendix  C),  a  series  of  eight  benchmark  reference  cases  were  proposed  by  Wickstrom  and 
Palsson  (1999)  for  use  in  verifying  two-dimensional  heat  conduction  computer  programs 
(Wickstrom  et  al,  1999;  Pilsson  et  al.,  2000).  These  entailed  thermal  calculations  of  increasing 
complexity.  The  reference  cases  were  developed  to  address  the  following  aspects  of  a  numerical 
heat  transfer  calculation: 

•  Convergence; 

•  Variable  convective  and  radiative  boundary  conditions; 

•  Variable  thermal  material  properties; 

•  Multiple  materials; 

•  Phase  change  of  material  constituents;  and 

•  Radiation  and  convection  in  voids. 

Results  from  the  Verification  Process  showed  that  HEATING  compared  successfully  with 
Reference  Cases  1  through  8  of  the  benchmark  proposed  by  the  developer  (Wickstrom  et 
al.,  1999).  It  was  found  that  for  cases  1  through  6, 64  cells  are  sufficient  for  satisfactory 
convergence  with  either  the  analytical  solution  or  with  the  results  calculated  by  the  computer 
code  “Temperature  Analysis  of  Structures  Exposed  to  Fire”  (TASEF).  For  Reference  Case  1, 
two  bottom  cells  suffice  to  yield  excellent  results  for  the  one-dimensional  transfer  case  (7a).  For 
the  two  dimensional  case  (7b),  8  bottom  cells  result  in  differences  smaller  than  3.2%.  For  the 
more  complex  geometry  of  Reference  Case  8,  the  results  for  the  most  refined  mesh  are  in 
excellent  agreement  with  TASEF  when  only  radiative  transfer  is  occurring  in  the  cavities.  When 
convective  heat  transfer  was  added  to  the  cavities,  HEATING  yielded  moderately  lower 
temperatures  than  TASEF. 


7.2 

Validation 

7.2.1 

Methodology 

As  with  “Verification”,  the  “Validation”  for  the  Methodology  is  a  composite  of  the 
“Validations”  for  CFAST  and  HEATING  that  are  described  briefly  below. 


7.2.2  CFAST 


Before  the  Navy’s  involvement  (circa.  1989),  V&V  work  on  CFAST  dealt  with  earlier 
versions  of  CFAST  (I.x  or  earlier).  The  modifications  to  CFAST  through  Navy  involvement  are 
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contained  in  versions  16.4, 2.x,  3.x,  4.0.1  and  5.0.1.  Versions  3.1.4  and  4.0.1  have  undergone 
the  most  extensive  validations,  using  experimental  data  documented  in  various  programs 
eonducted  aboard  ex-USS  SHAD  WELL  (Williams,  et  al,  2001)  in  Mobile  AL.  However,  the 
most  current  publicly  available  version  of  CFAST  is  5.0.1,  for  which  no  validations  have  been 
eompleted  to  date.  The  inherent  difficulty  in  using  any  version  of  the  model  is  that  phenomena 
implemented  and  validated  in  earlier  versions  are  not  re-verified  after  “fixes”  have  been  made  in 
later  versions.  On  the  other  hand,  it  is  very  possible  that  some  of  the  problems  encountered  by 
researchers  using  earlier  versions  of  CFAST  may  have  been  fixed  in  the  latest  version(s).  The 
V&V  requirements,  therefore,  eonsist  of  demonstrating  that  these  tools,  although  validated  over 
a  number  of  versions  of  the  methodology,  still  produce  aceeptable  results  and  that  the  results  are 
applicable  to  a  full-size  ship,  since  SHAD  WELL  represents  a  retired  naval  platform.  For  these 
reasons  (and  others  documented  later  in  this  Section)  CFAST  Version  3.1.7  is  the  fire  model 
recommended  for  assessment  in  the  PDA(F)  for  T-AKE  Test  and  Evaluation. 

The  Methodology  invokes  the  new  phenomena  of  vertical  vent  flow,  single  zone  (SHAFT) 
and  improved  gas  absorbanee  eoefficient  in  CFAST  Version  3.1.7.  The  validation  for  vertical 
vent  flow  using  SHADWELL  experimental  data  from  the  Smoke  Ejeetion  System  tests 
(Appendix  A)  showed  consistent  trends  in  eompartment  temperatures  between  CFAST 
predictions  and  SHADWELL  experiments.  The  absolute  temperature  eomparison  showed  an 
under-prediction  by  the  model  attributable  to  over-prediction  of  gas  flow  through  the  vent.  In  the 
fire  eompartment  over  a  90-minute  period,  the  experimental  temperature  range  was  150  deg  C  to 
400  deg  C.  The  predicted  model  temperature  in  that  eompartment  ranged  from  100  deg  C  to  300 
deg  C.  In  the  vent-connected  eompartment,  the  experimental  temperature  range  was  200  deg  C 
to  700  deg  C  and  the  model-predicted  temperature  range  was  200  deg  C  to  400  deg  C.  Vertieal 
vent  flow  phenomenon  was  the  first  to  be  added  to  enhance  the  performanee  of  CFAST.  Some 
faetors  that  may  have  contributed  to  the  only  fair  comparison  were  added  in  later  modifications 
of  CFAST,  e.g.,  vertical  and  horizontal  heat  transfer  and  long  eorridor  effects. 

CFAST  Version  3.1.7  also  invokes  the  improved  gas  absorbance  coefficient  and  the  single 
zone,  or  SHAFT,  option.  One  of  the  assumptions  originally  made  by  CFAST  was  the  use  of 
constant  (rather  than  time- variant)  extinetion  coefficients  in  the  upper  and  lower  layers  to 
ealculate  the  absorbance  or  radiation  effect  of  the  combustion  gases.  In  order  to  address  this 
deficieney,  CFAST  was  modified  to  permit  the  absorbance  of  each  layer  to  be  calculated  as  a 
function  of  the  time-variant  eoneentration  of  carbon  dioxide,  water  vapor  and  soot.  With  the 
new  algorithm,  the  total  absorbance  is  ealculated  analytically  from  the  gas  absorbance  and 
various  soot  pareuneters  and  the  gas  absorbance  is  estimated  by  interpolation  of  values  stored  in 
lookup  tables. 

The  SHAFT  option  was  used  in  the  validation  of  the  gas  absorbance  algorithm.  When  this 
option  is  invoked  on  a  particular  compartment,  the  model  uses  a  single  set  of  parameters 
(i.e.,  temperature  and  gas  concentration)  to  represent  the  entire  compartment  volume,  rather  than 
two  sets  (one  for  each  layer).  This  is  appropriate  whenever  the  compartment  gases  are  well 
mixed  and  may  therefore  be  represented  by  a  single  zone.  With  the  SHAFT  option  in  effeet,  the 
temperature  at  the  top  of  a  compartment  is  affected  by  heat  conducted  from  the  compartment 
below;  therefore,  vertical  energy  transfer  is  handled  correctly.  With  this  in  mind,  users  must  be 
aware  that  the  original,  hard-coded  extinction  value  used  for  the  upper  layer  in  any  version  of 
CFAST  before  Version  3.0  is  inappropriate  for  compartments  which  accumulate  little,  if  any. 
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combustion  products;  for  example,  those  eompartments  far  from  the  fire.  If  the  SHAFT  option  is 
used  for  one  of  these  eompartments,  the  upper  layer  extinction  will  be  applied  to  the  entire 
compartment  and  the  results  for  the  whole  compartment  will  be  in  error.  This  problem  will,  of 
course,  be  moot  if  CFAST  3.0  or  later  versions  are  used  because,  in  that  case,  the  extinction 
calculation  will  be  transient. 

Complete  validation  of  the  horizontal  heat  conduction  algorithm  requires  comparison  with 
additional  experimental  data.  In  the  Methodology  used  for  the  CVNX  VAR  (and  subsequently 
to  be  used  for  T-AKE),  the  inclusion  of  vertical/horizontal  heat  transfer  in  CFAST  is  not  as 
significant  as  other  discussed  phenomena  since  the  Methodology  invoked  to  spread  fires  through 
intact  bulkheads,  decks  and  overheads  is  via  the  finite  differenee  heat  conduction  model 
HEATING  Version  7.3.  These  give  the  additional  reasons  referred  to  earlier  as  to  why  CFAST 
Version  3.1.7  is  the  recommended  version  for  use  in  the  T-AKE  Test  and  Evaluation. 

The  results  from  the  fire  modeling  are  used  to  prediet  the  temperatures,  interface  heights 
(demarcation  between  smoke  and  relatively  smoke-free  zones)  and  smoke  concentrations  in  the 
PDA(F)  compartments.  The  temperature  data  are  used  to  determine  when  or  if  fire  will  spread 
within  these  compartments  as  well  as  to  other  spaces.  The  interface  height  and  smoke  density 
data,  in  addition  to  the  temperature  data,  are  used  to  determine  how  tenable  these  compartments 
would  be  for  protected  and  unprotected  personnel.  (The  criteria  used  for  tenability  are  not  part 
of  this  documentation).  This  information  is  the  most  basic  output  provided  by  the  fire  model. 

The  temperature  data  have  been  validated  and  verified  extensively  and  the  results  of  this  V&  V 
are  eontained  in  Appendix  A. 


7.2.3  HEATING 


HEATING  has  eompared  suceessfully  with  eight  validation  cases  of  increasing  complexity. 
Validation  cases  were  selected  based  on  the  availability  of  test  data  and  material  properties  for 
the  configurations  tested.  The  validation  basis  varied  from  case  to  case  depending  on  the 
particular  test  available,  and  included  fire  resistance  ratings,  temperatures  at  fixed  locations  at 
fixed  times,  and  transient  temperatures  at  fixed  locations.  Convergence  was  demonstrated  in 
each  case  by  increasing  the  number  of  nodes,  reducing  the  time  step,  reducing  the  tolerance,  and 
using  alternative  solution  methods. 

The  sensitivity  of  each  validation  case  to  the  input  parameters  was  also  investigated  because 
of  the  inherent  uncertainty  in  the  specific  material  properties,  especially  at  elevated  temperatures, 
and  as  well  as  the  boundary  condition  parameters.  An  uncertainty  of  plus  or  minus  25  percent 
was  assumed  where  applicable.  It  was  found  that  the  validation  cases  generally  separated  into 
three  categories:  those  sensitive  to  boundary  condition  parameters,  namely  the  effective 
emissivity  of  the  exposed  surface;  those  sensitive  to  the  material  properties,  in  particular  the 
thermal  conductivity  of  the  dominant  material,  and  those  that  were  insensitive  to  uncertainty  in 
either  parameter  group. 
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8.0  VERIFICATION  AND  VALIDATION  ADEQUACY 


The  methodology  described  in  this  report  is  intended  to  identify  specific  weaknesses  and/or 
shortcoming  that  might  decrease  ship  survivability.  The  methodology  is  summarized  in  the 
flowchart  in  Figure  1  (Section  6).  The  process  involves  modeling  PDA(F)  compartments  with 
the  multi-compartment  zone  fire  model  CFAST.  The  use  of  CFAST  allows  for  a  more  complete 
characterization  of  these  compartments.  The  output  from  CFAST  is  used  to  predict  fire  spread 
time.  Fire  spread  beyond  the  PDA(F)  compartments  to  APDA(F)  and  BAPDA(F)  compartments 
is  predicted  using  a  conservative  set  of  prescriptive  rules  in  conjunction  with  the  finite  element 
model  HEATING  to  evaluate  thermal  penetration  times. 

The  software  that  constitutes  CFAST  is  a  collection  of  data  and  computer  programs  that  are 
used  to  simulate  the  important  time-dependent  phenomena  involved  in  fires.  The  major 
functions  provided  include  calculation  of: 

•  The  buoyancy-driven  as  well  as  forced  transport  of  this  energy  and  mass  through  a  series 
of  specified  compartments  and  connections  (e.g.,  doors,  windows,  cracks,  ducts,  hatches, 
scuttles), 

•  The  resulting  temperatures,  smoke  optical  densities,  and  gas  concentrations  after 
accounting  for  heat  transfer  to  surfaces  and  dilution  by  mixing  with  clean  air. 

As  can  be  seen  from  this  list,  fire  modeling  involves  an  interdisciplinary  consideration  of 
physics,  chemistry,  fluid  mechanics,  and  heat  transfer.  In  some  areas,  fundamental  laws 
(conservation  of  mass,  energy,  and  momentum)  can  be  used,  whereas  in  others  empirical 
correlations  or  even  “educated  guesses”  must  be  employed  to  bridge  gaps  in  existing  knowledge. 
The  necessary  approximations  required  by  operational  practicality  result  in  the  introduction  of 
uncertainties  in  the  results.  The  user  should  understand  the  inherent  assumptions  and  limitations 
of  the  programs,  and  use  these  programs  judiciously  -  including  sensitivity  analyses  for  the 
ranges  of  values  for  key  parameters  -  in  order  to  make  estimates  of  these  uncertainties. 

CFAST  has  been  validated  extensively  by  NIST,  NRL  and  others.  Although  differences 
between  the  model  and  the  experiments  were  evident  in  the  studies,  they  are  typically  explained 
by  limitations  and  assumptions  of  the  model  and  of  the  experiments.  Like  all  predictive  models, 
the  best  predictions  come  with  a  clear  understanding  of  the  limitations  of  the  model  and  of  the 
inputs  provided  to  the  calculations. 

The  use  of  heuristies  for  far-field  predictions  provides  conservative/generic  guidelines  for 
fire  and  smoke  spread  outside  of  the  PDA(F).  These  prescriptive  rules  assume  that  an  unlimited 
amount  of  fuel  and  air  is  available  in  each  compartment  so  predictions  are  worst  case.  In  the 
event  that  there  is  a  small  amount  of  fuel  in  the  compartment  or  the  compartment  is  closed,  fire 
spread  may  be  delayed  or  may  not  occur  at  all.  The  finite  difference  heat  conduction  model  - 
HEATING  -provides  excellent  predictions  for  thermal  penetration  time  through  a  bulkhead  or 
deck.  The  limitations  imposed  by  CFAST  are  tempered  by  the  use  of  worst-case  scenarios  in 
conjunction  with  the  HEATING  model  to  make  the  overall  Methodology  approach  a 
conservative  estimate  of  the  secondary  fire  and  smoke  damage. 
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1.0  CHARACTERIZATION  OF  SECONDARY  DAMAGE  IN  THE  PDA(F) 

1.1  Development  and  Assumptions  of  the  Fire  Model 


In  the  current  methodology  (Back,  March,  2003),  the  multi-compartment  zone  fire  model 
used  for  characterizing  the  environment  in  primary  damage  area  that  includes  vent-cOnnected 
spaces  for  each  hit  scenario  is  Consolidated  Fire  and  Smoke  Transport  (CFAST)  Version  3.1.7 
(Peacock  et  al,  January  2000;  January  2000).  APDA(F)  and  BAPDA(F)  spaces  are  not  included 
in  the  fire  model.  CFAST  is  a  member  of  a  class  of  models  referred  to  as  zone  or  finite  element 
models.  This  means  that  each  room  is  divided  into  a  small  number  of  volumes  (called  layers), 
each  of  which  is  assumed  to  be  internally  uniform.  That  is,  the  temperature,  smoke  and  gas 
concentrations  within  each  layer  are  assumed  to  be  exactly  the  same  at  every  point.  In  CFAST, 
each  room  is  divided  into  two  layers.  Since  these  layers  represent  the  upper  and  lower  parts  of 
the  room,  conditions  within  a  room  can  only  vary  firom  floor  to  ceiling,  and  not  horizontally. 
This  assumption  is  based  on  experimental  observations  that  in  a  fire,  room  conditions  do  stratify 
into  two  distinct  layers.  While  we  can  measure  variations  in  conditions  within  a  layer,  these  are 
generally  small  compared  to  differences  between  the  layers. 

CFAST  is  based  on  solving  a  set  of  equations  that  predict  state  variables  (pressure, 
temperature  and  so  on)  based  on  the  enthalpy  and  mass  flux  over  small  increments  of  time. 
These  equations  are  derived  from  the  conservation  equations  for  energy,  mass,  and  momentum, 
and  the  ideal  gas  law.  The  conservation  equations  are  always  correct,  everywhere.  Thus  any 
errors  made  by  the  model  cannot  come  from  these  equations,  but  rather  come  from  simplifying 
assumptions  or  from  processes  having  been  omitted.  Examples  of  each  source  of  error  will  be 
highlighted  in  the  following  discussion. 


1.1.1  Fires 


CFAST  uses  a  specified  fire.  This  requirement  is  probably  the  most  important  restriction  on 
the  use  of  CFAST.  The  input  heat  release  rate  of  the  foel,  specified  by  the  user,  is  critical  to 
achieving  accurate  model  predictions  of  temperature  (one  of  the  primary  parameters  used  in  the 
validation  process).  This  fuel  is  converted  into  enthalpy  (the  conversion  factor  is  the  heat  of 
combustion)  and  mass  (the  conversion  factor  is  the  yield  of  a  particular  species)  as  it  bums. 
Burning  can  take  place  in  the  portion  of  the  plume  in  the  lower  layer  (if  any),  in  the  upper  layer 
or  in  a  door  jet.  For  an  unconstrained  fire  (i.e.,  not  oxygen-limited),  CFAST  always  uses  the 
specified  heat  release  rate.  Burning  will  all  take  place  within  the  fire  plume.  For  a  constrained 
fire,  the  heat  release  rate  can  be  reduced  below  the  user-specified  level  if  there  is  insufficient 
oxygen  to  support  that  combustion  rate.  Burning  will  take  place  where  there  is  sufficient 
oxygen.  In  the  simulation  analyses,  constrained  fires  are  typically  used  with  conservative  worst 
case  fire  growth  rates.  This  leads  to  oxygen-limited  fires  that  are  as  large  as  can  be  supported. 

In  general  this  results  in  a  conservative  prediction. 

Specification  of  the  heat  release  rate  by  the  user  means  that  this  version  of  CFAST  (3.7.1), 
like  all  current  zone  fire  models,  does  not  include  a  pyrolysis  model  to  predict  fire  growth. 
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Rather  pyrolysis  rates  for  each  fire  modeled  define  the  fire  history.  The  similarity  of  that  input  to 
the  real  fire  problem  of  interest  will  determine  the  accuracy  of  the  resulting  calculation.  The  user 
must  account  for  any  interactions  between  the  fire  and  the  pyrolysis  rate.  Generally,  more  severe 
fires  are  postulated. 

Another  issue  that  must  be  considered  carefully  when  executing  CFAST  is  the  fuel  load  since 
the  user  is  specifying  a  time-dependent  description  of  the  growth  of  the  fire.  Since  CFAST  does 
not  model  the  fire  itself  but  the  effects  of  the  fire  on  the  environment,  the  amount  of  fuel 
consumed  in  the  fire  is  based  on  the  user-specified  fire  growth  descriptions.  Back,  March  2003, 
and  Hoover  and  Tatem,  June  2000,  provide  excellent  guidelines  for  factors  to  be  considered  by 
the  user  in  making  this  selection.  Typically,  aggregate  fuel  loading  is  assumed  to  be  present 
since  it  is  difficult  to  preclude  the  presence  of  combustible  material. 


1.1.2  Vent  Flows 


Flow  through  vents  is  the  dominant  phenomenon  in  a  fire  model  because  it  fluctuates  most 
rapidly  and  transfers  the  greatest  amount  of  enthalpy  on  an  instantaneous  basis  of  all  the  source 
terms.  Also,  it  is  most  sensitive  to  changes  in  the  environment.  There  are  two  primary  means 
for  flow  through  vents.  The  first  is  referred  to  as  horizontal  gas  flow  through  vertical  vents.  It  is 
the  flow  that  is  normally  thought  of  when  discussing  fires.  It  encompasses  flow  through  doors, 
windows  and  other  vertical  openings.  The  other  is  vertical  gas  flow  through  horizontal  vents  and 
can  occur  if  there  is  a  hole  in  the  ceiling  or  floor  of  a  compartment  (e.g.  ship  hatches  and 
scuttles).  This  latter  phenomenon  is  particularly  important  for  shipboard  vent  flow  and  for  roof 
venting  in  structural  fires. 

Flow  through  vents  is  governed  by  the  pressure  difference  across  a  vent.  There  are  two 
situations  which  give  rise  to  flow  through  vents.  The  first,  and  usually  thought  of  in  fire 
scenarios,  is  that  of  air  or  smoke  which  is  transported  from  a  compartment  primarily  by 
buoyancy.  The  second  type  of  flow  is  due  to  expansion  that  is  particularly  important  when 
conditions  in  the  fire  environment  are  changing  rapidly.  Rather  than  depending  entirely  on 
density  differences  between  the  two  gases,  the  flow  is  forced  by  volumetric  expansion.  The 
earlier  versions  of  CFAST  did  not  solve  this  part  of  the  problem  entirely  correctly.  In  most  cases 
the  differences  are  small  except  for  rapidly  changing  situations.  However,  these  small 
differences  become  very  important  if  we  wish  to  follow  flows  due  to  small  pressure  differences, 
such  as  will  occur  in  a  mechanical  ventilation  system.  Atmospheric  pressure  is  about  100,000 
Pa;  fires  produce  pressure  changes  from  1  to  1,000  Pa  and  mechanical  ventilation  systems 
tj^jically  involve  pressure  differentials  of  about  1  to  100  Pa.  In  order  to  solve  these  interactions 
correctly,  we  must  be  able  to  follow  pressure  differences  of  0. 1  Pa  out  of  100,000  Pa  for  the 
overall  problem,  or  10*^  for  adjacent  compartments. 

Natural  ventilation  may  be  provided  by  holes  in  the  skin  of  the  ship  (resulting  from  the 
weapon  hit)  and  openings  between  PDA(F)  compartments.  At  a  minimum,  oxygen  for 
combustion  is  provided  naturally  through  any  openings  to  weather  that  result  from  the  weapon 
hit.  Forced  ventilation  may  be  provided  in  vital  spaces  where  ventilation  is  not  secured  under 
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materiel  condition  ZEBRA.  In  order  to  handle  the  flows  for  natural  and  forced  ventilation,  small 
pressure  differences  are  followed  in  later  versions  of  CFAST. 


1.1.3  Non-Navy-Funded  Verification  and  Validation  of  CFAST 


Jones  et  al.,  2000,  reviewed  a  number  of  studies  that  have  compared  the  predictions  from 
CFAST  to  real-scale  fires.  Nelson,  1989,  used  simple  computer  fire  models  along  with  existing 
experimental  data  to  develop  an  analysis  of  a  large  high-rise  building  fire.  This  analysis  showed 
the  value  of  available  analytical  calculations  in  reconstructing  the  events  involved  in  a  multiple- 
story  fire. 

Bukowski,  1990;  1992;  1995,  has  applied  the  FAST  and  CFAST  models  in  several  fatal  fire 
reconstructions.  Details  of  the  fires  including  temperatures,  vent  flows,  and  gas  concentrations 
were  consistent  with  observed  conditions  and  witness  accounts. 

Several  studies  comparing  model  predictions  with  experimental  measurements  are  available. 
Deal,  1990,  reviewed  four  computer  fire  models  (CCFM,  FIRST,  FPETOOL  and  FAST)  to 
ascertain  the  relative  performance  of  the  models  in  simulating  fire  experiments  in  a  small  room. 
All  the  models  simulated  the  experimental  conditions  including  temperature,  species  generation, 
and  vent  flows  “quite  satisfactorily”. 

Duong,  1990,  studied  the  predictions  of  several  computer  fire  models  (CCFM,  FAST, 

FIRST,  and  BRI),  comparing  the  models  with  one  another  and  with  large  fires  in  an  aircraft 
hanger.  For  a  4  MW  fire  size,  he  concluded  that  all  the  models  are  “reasonably  accurate.”  At  36 
MW,  however,  “none  of  the  models  did  well.”  Beard,  1990;  1992,  evaluated  four  fire  models 
(ASET,  FAST,  FIRST,  and  JASMINE)  by  modeling  three  well-documented  experimental  fires, 
ranging  in  scope  from  single  compartments  to  a  large-department-store  space  with  closed  doors 
and  windows.  He  provides  both  a  qualitative  and  quantitative  assessment  of  the  models  ability  to 
predict  temperature,  smoke  obscuration,  CO  concentration,  and  layer  interface  position  (for  the 
zone-based  models). 

Peacock  et  al.,  1993, 1998,  compared  the  CFAST  model  to  a  range  of  experimental  fires. 

The  model  provided  predictions  of  the  magnitude  and  trends  (time  to  critical  conditions  and 
general  curve  shape)  for  the  experiments  studied  which  range  in  quality  from  within  a  few  %  to  a 
factor  of  two  or  three  of  the  measured  values. 

Other  peer-reviewed  publications  have  also  compared  the  results  of  CFAST  with 
experimental  data.  For  completeness,  these  comparisons  are  reported  here  with  brief  synopses  of 
conclusions  drawn.  More  detailed  discussions  can  be  obtained  from  the  citations. 

Dembsey,  et  al.,  1995,  compared  predictions  from  CFAST  2.1  with  experiments  and  found 
the  predictions  of  upper  layer  gas  temperatures  to  be  higher.  The  authors  conjectured  that  the 
difference  was  due  to  insufficient  heat  transfer  out  of  the  simulated  compartment.  (Version  3.1.7 
of  CFAST  has  an  improved  vertical  heat  transfer  sub-model.) 
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Collier,  1995,  compared  experiments  in  an  actual  multi-room  house.  Comparisons  between 
the  predicted  interface  height  (CFAST  2.x)  and  the  compartments’  temperatures  with  actual  data 
provided  varying  degrees  of  agreement  depending  upon  the  time  and  compartment  being 
compared. 

Chow,  1996  made  multi-room  comparisons  of  the  upper  layer  temperature  and  interface 
height  with  predictions  from  CFAST  2.0.  Again,  there  are  varying  degrees  of  agreement 
depending  upon  the  time  and  the  compartment  being  compared. 

Davis,  et  al.,  1996,  compared  ceiling  jet  temperature  predictions  using  CFAST  2.0  to 
experimental  data  from  a  15  m  tall  hangar.  Results  from  the  predictions  were  consistently  higher 
than  those  of  the  experiments. 

Rockett,  1997,  evaluated  the  accuracy  of  plume  algorithms  in  CFAST  2.0.  He  compared 
CFAST  predictions  and  experimental  values  for  the  vent  flow,  dependent  upon  heat  release  rate. 

Luo,  1997,  compared  CFAST  2.0.1  and  experimental  data  for  interface  height  and  upper 
layer  temperatures  and  showed  the  predictions  to  be  higher  than  the  actual. 

He,  et  al.,  1998,  looked  at  two  methods  for  determining  interface  height  in  CFAST  (version 
unknown)  and  compared  both  methods  to  data  from  large  scale  experiments.  CFAST  predictions 
compared  poorly  with  actual  measurements  by  significantly  over-predicting  upper  layer 
temperatures. 


1.1.4  Navy-Funded  V erification  and  Validation  of  CFAST 

1 . 1 .4. 1  Natural  and  Forced  Vent  Flow  Validation 

1.1. 4. 1.1  General 


Smoke  ejection  system  (SES)  fire  tests  conducted  by  the  Naval  Surface  Warfare 
Center/Carderock  (CD/NSWC),  Code  633,  onboard  ex-USS  SHAD  WELL  were  used  to  validate 
the  concept  of  vertical  gas  flow  through  horizontal  vents  and  smoke  control  through  fan 
pressures  and  closures  (Bailey,  1994).  Because  of  the  presence  of  scuttles  and  hatches,  as  well 
as  the  presence  of  heating,  ventilation  and  air  conditioning  (HVAC)  systems  and  their  use  to 
control  smoke  movement,  both  of  these  concepts  are  important  phenomena  that  must  be 
accounted  for  in  naval  platforms.  The  experimental  data  from  the  SES  tests  were  compared  to 
predictions  from  CFAST  version  1.6.4. 

The  forward  area  of  ex-USS  SHADWELL,  frames  0  to  36,  was  modified  to  resemble  the  DDG 
5 1  design  criteria  (Figure  A-1).  The  test  area  consisted  of  a  machinery  space  connected  to  a 
mock  DC  passageway,  hereafter  referred  to  as  DC  passageway,  by  a  vestibule.  The  vestibule 
was  located  directly  above  the  machinery  overhead  and  was  divided  into  two  compartments  - 
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Fig.  A-1  -  SES  Test  Area 
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inboard  and  outboard  -  by  a  standard  Navy  archway.  The  fuel  for  the  machinery  space  fire  was  a 
diesel  spray,  which  was  ignited  by  a  heptane  pan  fire.  The  machinery  space  was  equipped  with 
both  a  supply  and  exhaust  fan.  The  DC  passageway  subdivisions  also  had  mechanical 
ventilation.  There  were  vertical  thermocouple  trees  in  the  machinery  space,  both  the  inboard  and 
outboard  vestibules,  and  the  DC  passageway.  These  temperature  measurements  from  these 
thermocouple  trees  were  compared  with  the  temperatures  predicted  in  the  compartment  by 
CFAST. 

1.1. 4. 1.2  Limitations: 

1 .  The  model  predicts  the  average  temperatures  in  an  upper  and  lower  layer,  whereas  the 
experiments,  in  some  cases,  measure  a  single  point.  In  other  cases  where  several  points 
are  measured,  an  average  layer  temperature  must  be  determined  based  on  limited 
individual  measurements  at  discrete  locations. 

2.  The  accuracy  of  the  model  predictions  depends  on  accurate  inputs.  The  unavailability 
of  any  experimental  setup  information  may  have  led  in  this  V&V  exercise  to  best 
“guess-timates”  for  input  parameters. 

In  the  mflchinerv  space,  all  of  the  thermocouple  data  were  averaged  to  give  an  average 
compartment  temperature.  These  temperatures  were  compared  to  a  volume-average  of  the 
model-predicted  temperatures.  With  the  horizontal  vent  open  in  the  space,  the  model  and 
experiment  both  showed  that  the  compartment  temperature  increases  when  the  fire  size  is 
increased  and  decreases  when  fire  size  is  decreased.  Though  the  temperature  trends  are 
consistent,  the  absolute  temperature  is  under-predicted  by  the  model  when  the  horizontal  vent  is 
open.  These  observations  imply  that  the  model  is  over-predicting  the  flow  through  the  horizontal 
vent.  This  over-prediction  is  seen  consistently  when  comparing  temperatures  in  the  DC 
passageway. 

The  inboard  vestibule  was  connected  to  the  machinery  space  via  a  horizontal  vent  and  to  the 
outboard  vestibule  by  an  archway.  In  the  inboard  vestibule  the  model  and  experimental  data 
show  a  substantially  higher  temperature  relative  to  the  machinery  space.  Model-predicted  and 
experimentally  calculated  temperatures  also  show  the  same  increasing  and  decreasing  trends. 

The  absolute  temperatures  are  close  during  the  first  40  minutes  of  the  experiment  but  then  begin 
to  deviate  substantially. 

In  the  outboard  vestibule,  the  model-predicted  and  experimentally  determined  temperature 
profiles  were  very  similar  to  their  counterparts  in  the  inboard  vestibule.  Both  the  model  and  the 
experiment,  however,  show  that  the  temperatures  in  the  outboard  vestibule  are  lower  than  those 
in  the  inboard  vestibule  throughout  the  experiment. 

In  the  DC  passageway,  the  model  over-predicted  the  temperatures  relative  to  the 
experimentally  determined  values.  Some  factors  that  may  contribute  to  this  only  fair  comparison 
have  been  addressed  in  later  modifications  to  CFAST,  e.g.,  horizontal  and  vertical  heat  transfer 
and  “corridor”  effects  for  long  passageways.  The  improved  corridor  sub-model,  which  is 
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discussed  later,  is  transient  -  it  only  affects  the  upper  layer  during  the  initial  ceiling  jet  flow. 
After  that,  everything  reverts  to  the  old  algorithm  for  handling  smoke  movement.  Other  factors 
that  still  need  to  be  addressed  are  horizontal  placement  of  meehanical  ventilation  openings  and 
flame  extension  through  horizontal  vents. 


1.1.5  Vertical  Heat  Transfer  Algorithm  (Heat  Conduction  Through  Decks) 
1.1.5. 1  General 


The  vertical  heat  transfer  algorithm  was  added  to  CFAST  to  correctly  account  for  conductive 
heat  transfer  in  the  vertical  direction  through  metal  decks  in  naval  platforms.  The  addition  of 
this  phenomenon  in  shipboard  environments  is  necessary  since  rising  temperatures  in  adjacent 
steel  compartments  as  a  result  of  fire  can  cause  materials  to  reach  ignition  levels  even  though  no 
breach  has  occurred  in  the  initial  fire  compartment. 

Comparisons  were  made  using  CFAST  version  2. 1  (with  the  vertical  heat  transfer  turned  on 
[2.1  +]  and  off  [2.1])  and  experimental  data  fi:om  real  scale  fire  tests  conducted  onboard  ex-USS 
SHAD  WELL  in  the  Internal  Ship  Conflagration  Control  (ISCC)  test  area  in  the  port  wing  wall 
(Figure  A-2)  (Bailey  and  Tatem,1995). 

The  compartments  that  comprised  the  ISCC  test  area  (Figure  A-2)  were  Berthing  2,  Ricer  2, 
CIC  and  Pilot  House.  The  experimental  conditions  were  produced  by  a  post-flashover  fire  in 
Berthing  2,  created  by  diesel  spray  fires  ignited  by  burning  heptane  pool  fires.  The 
measurements  of  interest  in  the  validation  process  were  the  compartment  temperatures  as 
measured  by  thermocouples.  Since  CFAST  is  a  zone  model,  it  divides  the  compartments  into 
two  layers  -  an  upper  and  a  lower.  Each  layer  is  assumed  to  be  of  uniform  temperature  and 
composition.  The  same  limitations  described  in  this  report  under  the  Section  "Natural  and 
Forced  Vent  Flow  Validation”  for  model  and  experimental  temperatures  apply  to  those 
parameters  in  this  validation. 

Model  assumptions  and  experimental  limitations  must  be  specified  for  accurate  interpretation 
of  the  comparisons.  A  limitation  that  existed  in  the  experiment  was  accurate  assessment  of  the 
vent  openings.  Inherent  in  the  aftermath  of  burning  are  the  formation  of  cracks  in  the  structure 
that  must  be  approximated  as  “vent  openings”.  Also  the  mass  loss  rate  of  the  fuel  was 
approximated  based  on  the  use  of  both  diesel  spray  and  heptane  pool  sources. 

The  model  assumes  that  each  compartment  is  a  rectangular  parallelpiped.  Since  none  of  the 
compartments  were  rectangular  parallelpipeds,  their  dimensions  were  adjusted  to  conserve  actual 
compartment  volume.  Also  the  compartment  walls  varied  in  thickness  within  a  given 
compartment.  A  weighted  average  based  on  actual  surface  area  was  used  as  input  to  the  model 
for  these  compartments.  The  model  uses  the  wall  thickness  to  calculate  conductive  heat  transfer 
through  the  walls. 


A-9 


PILOTHOUSE 

01-8M-C 


A-10 


1.1. 5.2  Results 


Prior  to  incorporation  of  the  algorithm  into  CFAST  which  properly  accormts  for 
vertical  heat  transfer  by  conduction,  the  model  predicted  that  the  temperatures  in  the 
compartments  and  decks  above  a  fire  compartment  would  not  increase  above  ambient. 
After  incorporating  the  algorithm,  the  model  predicted  much  more  realistic  results.  These 
predictions  compared  reasonably  well  with  experimental  results  for  the  fire  compartment 
(Berthing  2)  as  well  as  for  the  deck  and  compartment  directly  above  it  (Ricer  2).  The 
model  over-predicted  the  temperatures  of  the  compartments  and  decks  not  directly 
adjacent  to  the  fire  compartment  (CIC  and  Pilot  House).  For  Berthing  2,  both  the  model 
and  experimental  data  results  show  that  almost  the  entire  compartment  was  in  the  upper 
layer  during  the  experiments.  The  temperature  predictions  fi'om  version  2. 1+  were  higher 
than  those  firom  version  2.1.  This  is  because  in  version  2.1  the  model  assumes  the  energy 
is  conducted  through  the  fire  compartment  overhead  into  the  ambient  environment  which 
remains  at  a  constant  temperature.  In  version  2. 1+,  the  heat  is  conducted  into  Ricer  2 
which  heats  up  as  energy  is  accumulated  in  the  compartment. 

In  Ricer  2.  version  2.1  of  CFAST  predicts  that  the  deck  will  remain  at  ambient  (room) 
temperature  throughout  the  experiment.  Clearly,  the  predictions  firom  the  enhanced 
model  (2.1+)  are  much  more  realistic.  Version  2.1+  predictions  for  Ricer  2  air 
temperatures  compared  reasonably  well  at  the  beginning  and  the  end  with  experimental 
data,  but  are  different  by  about  50  degrees  Celsius  in  the  middle  portion.  However,  since 
version  2.1  does  not  contain  the  vertical  heat  transfer  algorithm  the  comparison  between 
model  and  experiment  was  very  poor  throughout  the  entire  time  period. 

The  largest  over-prediction  by  CFAST  2.1+  as  compared  to  experimental  data  was  in 
CIC  deck.  There  was  also  an  over-prediction  in  the  air  temperature  of  CIC  but  not  of  the 
same  magnitude.  This  could  be  a  result  of  the  over-prediction  in  the  Ricer  2  air 
temperature  and  also  the  representative  experimental  temperature  data  for  the  CIC  deck 
(one,  rather  than  several,  thermocouple  readings). 

The  deck  and  air  temperatures  in  the  Pilot  House  are  over-predicted  also,  but  is 
viewed  as  relatively  insignificant  since  at  these  locations  the  increase  above  ambient  is 
quite  small  for  both  the  predictions  and  the  experiments. 

Additional  validations  were  performed  with  the  ISCC  SHAD  WELL  data  and 
CFAST  3. 1 .  In  these  comparisons,  CFAST  3. 1  consistently  under  predicted  the 
temperatures  (Bailey,  et  al.,  1998).  In  the  second  validation,  in  order  to  model  Ricer  2, 
the  CIC  and  the  Pilot  House  as  one  zone  each,  the  SHAFT  option  was  used  (see  Section 
1.1.6). 

1.1.6  Improved  Radiation  Transport  Sub-Model 


One  of  the  assumptions  originally  made  by  CFAST  was  the  use  of  constant  (rather 
than  time-  variant)  extinction  coefficients  in  the  upper  and  lower  layers  to  calculate  the 
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absorbance  or  radiation  effect  of  the  combustion  gases.  In  order  to  address  this 
deficiency,  CFAST  was  modified  to  permit  the  absorbance  of  each  layer  to  be  calculated 
as  a  function  of  the  time- variant  concentration  of  carbon  dioxide,  water  vapor  and  soot 
(Hoover  et  al.,  1997;  Hoover,  1995;  Hoover  and  Bailey,  2000.). 

With  the  new  algorithm,  the  total  absorbance  is  calculated  analytically  from  the  gas 
absorbance  and  various  soot  parameters  and  the  gas  absorbance  is  estimated  by 
interpolation  of  values  stored  in  lookup  tables.  This  absorbance  functionality  was 
incorporated  into  CFAST  version  2.2.1.  Predictions  from  the  modified  CFAST,  known 
in  this  document  as  CFAST  2.2. 1+,  were  compared  to  those  fi'om  CFAST  2.2.1  and, 
where  possible,  with  data  from  experiments  on  ex-USS  SHAD  WELL  fire  tests  in  the 
Integrated  Ship  Conflagration  Control  (ISCC)  test  area.  This  experimental  description 
and  configuration  is  given  in  this  report  under  Section  “Vertical  Heat  Transfer  Algorithm 
(Heat  Conduction  Through  Decks)  (Hoover  et  al.,  1996). 

The  compartments  that  comprised  the  ISCC  test  area  were  Berthing  2.  Ricer  2.  CIC 
and  the  Pilot  House.  Berthing  2  was  handled  as  a  “two-zone”  space,  comparing  the 
experimental  upper  layer  temperatures  to  CFAST  2.2.1  and  CFAST  2.2. 1+.  There  were 
no  experimental  temperature  data  for  the  lower  layer,  so  only  the  predictions  from  the 
two  CFAST  model  versions  could  be  compared.  There  was  very  good  agreement 
between  the  predicted  upper  layer  air  temperatures  for  the  two  CFAST  model  versions 
and  the  corresponding  experimental  values.  The  lower  layer  of  Berthing  2  is  the  zone  in 
which  the  largest  differences  among  the  various  models  are  seen.  The  changes  in  the 
temperature  predictions  between  the  two  models  are  consistent  with  expectations  but, 
since  there  are  no  experimental  data  for  this  zone,  there  is  no  way  to  determine  which  of 
the  models  is  better. 

In  general,  errors  in  the  predictions  accumulate  as  one  moves  farther  from  the 
compartment  of  origin  and  it  becomes  much  harder  to  make  meaningful  comparisons 
with  experimental  data.  However,  some  useful  conclusions  can  be  drawn  from  the  results 
obtained  in  this  particular  validation. 

For  the  air  temperature  in  Ricer  2.  CFAST  2.2. 1+  is  more  accurate  than  CFAST  2.2. 1 , 
while,  in  CIC.  there  is  little  difference  between  the  two,  as  long  as  a  reasonable  value  of 
the  calculated  total  absorbance  is  used.  In  the  case  of  the  Pilot  House,  conditions  are  so 
close  to  ambient  that  there  is  little  effect,  regardless  of  which  model  or  what  value  of  the 
calculated  total  absorbance  is  used.  Inconsistencies  observed  in  the  Berthing  2  upper 
layer  air  temperatures  and  the  Ricer  2  deck  temperature  (common  boundary)  were 
possibly  related  to  the  mechanism  for  calculating  horizontal  heat  transfer  through  the 
bulkheads.  This  phenomenon  is  accounted  for  in  a  higher  version  of  CFAST. 

The  SHAFT  option  was  used  in  this  validation  of  the  gas  absorbance  algorithm. 

When  this  option  is  invoked  on  a  particular  compartment,  the  model  uses  a  single  set  of 
parameters  (i.e.,  temperature  and  gas  concentration)  to  represent  the  entire  compartment 
volume,  rather  than  two  sets  (one  for  each  layer).  This  is  appropriate  whenever  the 
compartment  gases  are  well  mixed  and  may  therefore  be  represented  by  a  single  zone. 
With  the  SHAFT  option  in  effect,  the  temperature  at  the  top  of  a  compartment  is  affected 
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by  heat  conducted  from  the  compartment  below;  therefore,  energy  transfer  from  Ricer  2 
to  CIC  and  CIC  to  the  Pilot  House  is  handled  correctly.  With  this  in  mind,  users  must  be 
aware  that  the  original,  hard-coded  extinction  value  used  for  the  upper  layer  in  CFAST 
2.2. 1  is  inappropriate  for  compartments  which  accumulate  little,  if  any,  combustion 
products;  for  example,  those  compartments  far  from  the  fire.  If  the  SHAFT  option  is 
used  for  one  of  these  compartments,  the  upper  layer  extinction  will  be  applied  to  the 
entire  compartment  and  the  results  for  the  whole  compartment  will  be  in  error.  This 
problem  will,  of  course,  be  moot  if  CFAST  2.2.1+  or  higher  is  used  because,  in  that  case, 
the  extinction  will  be  calculated  on  the  fly. 

One  final  consideration:  when  CFAST  2.2.1+  or  higher  is  used,  it  is  vital  that  some 
reasonable  value  of  the  total  gas  absorbance  be  specified  in  the  input  file  (i.e.,  non-zero). 
If  no  value  is  explicitly  defined,  a  value  of  zero  will  be  assumed  and  this  results  in  the 
soot  in  the  model  being  “turned  off’,  which  is  rarely  a  reasonable  thing  to  do. 


1.1.7  Corridor  Flow  or  Hybrid  Smoke  Movement  (CFD  Correlation)  Sub-Model 


The  Corridor  Flow  or  Hybrid  Smoke  Movement  Sub-model,  developed  by  NIST,  was 
validated  by  NRL  using  an  experiment  from  the  SHAD  WELL/688  test  series  and  CFAST 
versions  4.0.1  and  3.1.3  (Bailey,  et  al.,  2002,  Bailey  and  Tatem,  1998).  Prior  to  addition 
of  the  new  sub-model,  CFAST  did  not  properly  account  for  the  time  it  took  a  ceiling  jet, 
entering  at  one  end  of  a  corridor,  to  travel  down  to  the  other  end.  Previously,  CFAST 
assumed  an  isothermal  upper  layer  would  immediately  cover  the  entire  ceiling  area  along 
the  entire  length  of  the  corridor.  The  new  sub-model  enables  CFAST  to  predict  the  flow 
of  the  ceiling  jet  down  the  corridor  with  time.  CFAST  also  calculates,  at  any  given  time, 
the  temperature  distribution  of  the  gases  inside  the  ceiling  jet.  Results  presented  here  are 
predictions  from  the  enhanced  version  of  CFAST  (4.0.1)  that  contains  the  new  sub¬ 
model  and  to  experimental  data  from  the  SHAD  WELL/688  test  series. 

A  section  of  the  port  wing  wall  on  ex-USS  SHAD  WELL  was  configured  to  represent 
the  forward  compartment  of  an  SSN  688  Class  submarine  (Figure  A-3).  The  test  area 
was  comprised  of  14  compartments,  located  on  four  levels.  The  compartments  were 
connected  with  various  hatches,  watertight  doors  and  a  mechanical  ventilation  system. 
The  experiments  consisted  essentially  of  igniting  a  pan  fire  in  one  of  the  compartments 
and  then  manipulating  variables  to  observe  the  effect  on  the  smoke  flow  within  the  test 
area. 

There  was  only  one  compartment  within  the  test  area,  the  Laundry  Room,  which 
contained  a  corridor.  Both  exhaust  and  supply  terminals  were  located  near  the  ceiling  in 
the  corridor  of  this  compartment.  The  Corridor  Flow  Sub-model  does  not  take  into 
account  any  effects  of  forced  ventilation  on  the  flow  of  the  ceiling  jet.  Therefore,  an 
experiment  chosen  for  validation  purposes  did  not  employ  the  mechanical  ventilation 
system. 


A-13 


The  Laundry  Room  was  actually  subdivided  into  two  sub-compartments.  The  sub¬ 
compartment  which  contained  the  fire  is  the  Fire  Compartment  and  the  remaining  sub¬ 
compartment  is  the  Passageway.  The  Fire  Compartment  contained  a  fire  pan  filled  with  diesel 
fuel,  facilitated  in  its  ignition  with  a  small  amount  of  heptane.  The  temperature  comparisons 
used  for  CFAST  predictions  were  obtained  from  thermocouples  located  in  the  Fire  Compartment 
and  at  two  locations  referred  to  as  A  and  B,  in  the  Passageway. 

The  assumptions  that  had  to  been  made  for  these  comparisons  were  similar  to  the  previous 
validations  reported:  mass  loss  rate  estimated  based  on  temperature  data  from  inside  the  Fire 
Compartment;  translation  of  experimental  temperatures,  measured  at  distinct  locations,  to  an 
overall  zone  temperature  as  predicted  by  CFAST;  and  conservation  of  compartment  volume  due 
to  non-rectangular  geometry  in  the  test  area. 

The  experiment  to  be  modeled  began  with  ignition  of  a  diesel  pan  fire  at  time  zero.  Times 
were  noted  for  (1)  time  of  pan  fire  ignition,  (2)  time  at  which  combustion  products  flowed  into 
the  Passageway  through  the  Fire  Compartment  door,  (3)  time  at  which  ceiling  jet  reached 
Location  A,  (4)  time  at  which  ceiling  jet  reached  Location  B,  and  (5)  time  at  which  ceiling  jet 
reached  the  end  of  the  Passageway.  Experimental  times  for  Events  1, 3,  and  4  can  be  obtained 
from  thermocouple  readings  and  were  0, 27,  and  32  s,  respectively.  CFAST  4.0.1,  without  the 
Corridor  Flow  Sub-model  invoked,  assumes  that  a  thin  hot  gas  layer  forms  along  the  length  of 
the  entire  corridor  ceiling  the  instant  that  the  hot  gases  pass  under  the  Fire  Compartment  door 
soffit,  so  there  is  no  transient  nature  in  the  smoke  movement  along  the  long  Passageway.  The 
enhanced  model  (CFAST  4.0. 1  with  the  Corridor  Flow  Sub-model  invoked)  predicted  transient 
smoke  movement  so  that  Events  3  and  4  occurred  at  29  and  40  s,  respectively.  Therefore,  the 
enhanced  model  provides  a  more  realistic  prediction  of  the  delay  that  occurs  as  the  ceiling  jet 
travels  down  the  Passageway. 

Experimental  temperatures,  as  well  as  time,  were  correlated  with  the  original  and  enhanced 
model  predictions.  There  is  remarkably  good  agreement  between  the  temperatures  of  the 
experiment  and  the  enhanced  model  at  Location  A,  whereas  the  predicted  temperature  from  the 
original  model  deviates  significantly.  At  Location  B,  the  enhanced  model  agrees  better  with  the 
experiment  than  the  original  model,  but  the  agreement  is  not  as  good  as  at  Location  A. 


1.1.8  Horizontal  Heat  Transfer  Algorithm  (Heat  Conduction  Through  Bulkheads) 


The  horizontal  heat  transfer  algorithm  was  added  to  CFAST  to  correctly  account  for 
conductive  heat  transfer  in  the  horizontal  direction  through  metal  bulkheads  in  naval  platforms 
(Jones  et  al.,  2002).  The  addition  of  this  phenomenon  in  shipboard  environments  is  necessary 
since  rising  temperatures  in  adjacent  steel  compartments  as  a  result  of  fire  can  cause  materials  to 
reach  ignition  levels  even  though  no  breach  has  occurred  in  the  initial  fire  compartment. 

Heat  transfer  between  horizontally  connected  compartments  is  modeled  by  merging  the 
connected  surfaces  for  the  connected  horizontal  compartments.  A  heat  conduction  problem  is 
solved  for  the  merged  walls  using  a  temperature  boundary  condition  for  both  the  near  and  far 
wall.  Temperatures  are  determined  by  the  solver  so  that  the  heat  flux  striking  the  wall  surface 
(both  interior  and  exterior)  is  consistent  with  the  temperature  gradient  at  that  surface. 
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For  horizontal  heat  transfer  between  compartments,  the  connections  can  be  between  partial 
wall  surfaces,  expressed  as  a  fraction  of  the  wall  surface.  CFAST  first  estimates  conduction 
fractions  analogous  to  radiative  configuration  factors.  For  example,  a  conduction  fraction 
between  a  rear  wall  in  room  1  and  a  front  wall  in  room  2  is  the  heat  flux  fraction  from  the  room  2 
wall  that  strikes  and  hence  contributes  to  room  1  ’s  heat  transfer. 

CFAST  version  4.0.1  (and,  by  default.  Version  5.0  also)  contains  the  horizontal  heat 
conduction  (wall/wall)  phenomenon  as  well  as  the  Hybrid  Smoke  Movement  Sub-Model,  as 
stated  above.  Although  Version  4.0.1  was  validated  for  the  corridor  flow  algorithm,  there  was 
not  sufficient  experimental  data  in  any  of  the  ex-SHADWELL  fire  test  series  to  provide  surface 
temperature  readings  to  be  used  for  complete  validation  of  the  horizontal  heat  conduction 
algorithm  in  CFAST  Version  4.0. 1 .  A  preliminary  validation  of  the  horizontal  heat  transfer 
algorithm  was  conducted  using  CFAST  4.0.1  and  ISCC  experiments  (See  Section  “Vertical  Heat 
Transfer  Algorithm”).  This  Version  of  CFAST  still  requires  more  extensive  experimental 
validation  for  this  phenomenon. 
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A  METHODOLOGY  FOR  PREDICTING  FIRE  AND 
SMOKE  SPREAD  FOLLOWING  A  WEAPON  HIT 
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1.  Enclosure  (1)  is  forwarded  for  your  information  and  comment. 

2.  The  purpose  of  this  report  is  to  generically  document  the  fire  and  smoke  spread  methodology 
used  for  the  CVNX  Vulnerability  Assessment  Report  (VAR)  and  to  make  recommendations  for 
future  methodology.  Particular  emphasis  is  placed  on  the  definitions  used  for  the  PDA  (Primary 
Damage  Area)  and  the  APDA  (Adjacent  to  Primary  Damage  Area). 

3.  The  definitions  for  the  PDA  and  APDA  have  varied  depending  on  the  application.  We 
recommend  that  the  PDA  for  fire  and  smoke  spread  he  designated  as  PDA(F)  (Primary  Damage 
Area  (Fire)).  Since  fire  and  smoke  spread  and  the  availability  of  oxygen  for  combustion  is 
controlled  by  vents,  we  recommend  that  the  PDA(F)  spaces  include  any  adjoining  spaces  that 
have  free  communication  including  open  doors  and  scuttles.  These  spaces  are  in  addition  to 
those  that  experience  weapon  induced  structural  alternatives  such  as  bulkhead  failure  or  holing  in 
deck,  overhead  and  bulkheads. 

The  APDA(F)  would  then  by  definition  be  all  spaces  (including  above  and  below)  that  share  a 
common  deck,  bulkhead  or  overhead  with  PDA(F)  spaces.  The  definition  for  BAPDA(F)  would 
be  the  compartments  once  removed  from  the  APDA(F). 

4.  The  Naval  Research  Laboratory’s  point  of  contact  is  Dr.  Frederick  W.  Williams,  Code  6180: 
(202)  767-2002/e-mail:  fwilliam@ccs.nrl.navy.mil. 
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A  METHODOLOGY  FOR  PREDICTING  FIRE  AND  SMOKE 
SPREAD  FOLLOWING  A  WEAPON  HIT 

1.0  INTRODUCTION 


A  methodology  has  been  developed  and  refined  to  predict  smoke  and  fire  spread  after  a 
weapon  hit.  This  methodology  was  developed  to  comply  with  the  congressionally  mandated 
Live  Fire  Testing  and  Evaluation  (LFT&E)  program.  Under  the  LFT&E  program  for  each  ship 
platform,  a  vulnerability  assessment  is  undertaken  and  a  Vulnerability  Assessment  Report  (VAR) 
is  prepared.  The  VAR  provides  an  evaluation  of  the  survivability  of  the  ship  platform  and  ship 
mission  against  simulated  weapon  threats.  The  vulnerability  assessment  process  includes  the 
simulation  of  weapon  damage,  characterization  of  secondary  damage  (e.g.,  flooding,  fire  and 
smoke),  and  description  of  the  battle  damage  repair  activities  and  effectiveness.  These  activities 
are  generally  referred  to  as  a  Total  Ship  Survivability  (TSS)  and  Battle  Damage  Repair 
Assessment  (BDRA).  The  overall  objective  of  the  TSS/BDRA  is  to  identify  specific  weaknesses 
and/or  shortcomings  as  well  as  high  payoff  research  and  development  efforts  that  might  increase 
the  survivability  of  the  particular  ship  platform. 

The  current  TSS/BDRA  process  involves  consideration  of  various  weapon  hits.  For  each 
specific  weapon  threat  scenario  (shot  line),  the  initial  damage  is  identified  by  utilizing  a  blast 
damage  model.  Secondary  damage,  such  as  fire  and  smoke  spread  that  results  from  the  blast 
damage,  must  also  be  characterized.  The  methodology  described  in  this  document  outlines  the 
framework  through  which  a  conservative  estimate  of  the  secondary  fire  and  smoke  damage  can 
be  established.  This  document  also  uses  a  refined  set  of  definitions  for  fire  and  smoke  spread  in 
damaged  areas  of  the  ship. 


2.0  APPROACH 


The  current  approach  for  the  characterization  of  secondary  damage  associated  with  a  weapon 
hit  has  evolved  over  several  TSS/BDRA  evaluations.  The  genesis  of  the  current  methodology 
was  rudimentary  fire  spread  times,  developed  by  the  Naval  Research  Laboratory  (NRL)  and 
Hughes  Associates,  Inc.  based  on  USS  Stark  incident  data  and  testing  [1].  This  information  was 
implemented  in  a  TSS  model  developed  by  the  Naval  Surface  Warfare  Center,  Carderock 
Division  (NSWC/CD).  The  fire  spread  methodology  was  improved  and  the  framework  was 
expanded  to  allow  for  smoke  spread  during  the  evaluation  performed  for  the  LPD  17  VAR  2  [2]. 
Most  recently,  the  methodology  was  further  refined  during  the  CVNX  assessment  [3].  This 
report  reflects  the  most  current  methodology  used  for  the  CVNX  assessment.  The  framework 
incorporates  a  combination  of  computer  fire  modeling  using  the  Consolidated  Model  of  Fire 
Growth  and  Smoke  Transport  (CFAST)  program  and  prescriptive  rules  for  fire  spread.  CFAST 
is  specifically  used  to  characterize  the  environment  in  the  primary  damage  area  which  includes 
vent  connected  spaces. 

Estimates  of  the  worst  plausible  fire  spread  scenario  for  each  hit  are  used  to  determine  the 
extent  of  damage  control  (DC)  operations  that  are  required  to  control  and  secure  the  fire.  As  part 
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of  this  conservative  assessment  framework,  the  effects  of  fixed  fire  suppression  systems  or 
manual  firefighter  activities  normally  associated  with  DC  operations  are  not  considered  in  the 
initial  secondary  damage  assessment.  In  some  instances,  intact  suppression  systems  may  be 
evaluated  on  a  case  by  case  basis.  Estimates  of  the  rate  of  smoke  production,  extent  of  smoke 
migration  and  the  temperature  of  the  smoke  layer  are  used  to  identify  the  potential  exposures  to 
shipboard  equipment,  personnel  and  DC  parties. 

Describing  the  methodology  for  predicting  fire  and  smoke  spread  requires  an  understanding 
of  the  shipboard  damage  zones  developed  for  these  evaluations.  The  zones  are  categorized  in 
three  ways; 

Primary  Damage  Area  (PDAI  compartments:  PDA  compartments  are  those  that  incur 
significant  damage  from  the  weapon  hit  (e.g.,  failed  deck/bulkhead  and  shock  holing).  It  should 
be  noted  that  compartments  that  are  only  vent  connected  to  PDA  compartments  as  a  result 
of  door/hatch  failure  have  not  been  considered  part  of  the  PDA  in  certain  contexts.  This 
definition  differs  from  that  used  in  Damage  Control  Automation  for  Reduced  Manning  (DC- 
ARM)  testing  [4-6]  and  Weapons  Effect  Test  (WET)  for  the  ex-Caron  fire  and  smoke  spread  [7], 
where  all  compartments  that  are  vent  connected  as  a  result  of  the  blast  damage  are  considered 
PDA  compartments.  From  the  standpoint  of  fire  and  smoke  spread,  including  the  vent  connected 
spaces  makes  more  sense.  This  is  due  to  the  fact  that  the  vents  provide  a  means  for  spreading 
fire  and  smoke  directly  as  well  as  supplying  additional  oxygen  for  the  fires.  Therefore,  in  the 
future  we  recommend  designating  the  fire  and  smoke  spread  PDA  (noted  as  PDA(F)) 
where  vent  connected  spaces  are  included. 

Adjacent  to  Primary  Damage  Area  (APDAI  compartments:  These  compartments  share  a 
common  boundary  (bulkhead,  overhead  or  deck)  with  a  PDA  compartment.  Historically  these 
compartments  included  spaces  with  a  communicable  opening  (e.g.,  vent  connection)  between 
PDA  compartments  and  the  APDA  compartment,  such  as  a  failed  door  or  hatch.  Our 
recommended  definition  for  fire  and  smoke  spread  for  adjacent  compartments,  APDA(F), 
excludes  any  compartments  that  have  vent  openings  to  the  PDA.  This  is  a  result  of  those  vented 
compartments  becoming  part  of  the  PDA(F). 

Beyond  Adjacent  to  Primary  Damage  Area  (BAPDAI  compartments:  These  compartments 
border  APDA  spaces.  They  can  be  otherwise  described  as  “twice  removed”  from  the  PDA. 

These  spaces  may  be  a  concern  due  to  potential  fire  spread  later  into  the  event.  It  is  typically 
assumed  that  there  are  no  communicable  openings  (i.e.,  vents,  ducts,  open  doors  or  open  hatches) 
between  BAPDA  and  APDA  spaces.  This  definition  for  BAPDA(F)  remains  the  same  for  fire 
and  smoke  concerns  as  weapon  effect  concerns. 

For  illustrative  purposes.  Figure  B-1  shows  an  example  of  the  compartment  categorization  of 
spaces  for  a  hypothetical  weapon  hit  using  the  PDA,  APDA  and  BAPDA  definitions  (as 
compared  to  PDA(F),  APDA(F)  and  BAPDA(F)).  As  indicated  in  the  legend,  PDA 
compartments  are  marked  in  red,  APDA  compartments  are  marked  in  yellow  and  BAPDA 
compartments  are  marked  in  blue.  The  PDA  consists  of  those  compartments  that  are  structurally 
affected  by  the  blast  overpressure.  The  APDA  compartments  on  the  Third  Deck  and  First 
Platform  are  horizontally  adjacent  to  the  PDA  on  those  decks.  In  addition,  there  are  APDA 
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First  Platform 


Fig.  B-1  •  Example  of  compartment  damage  characterization  after  hypothetical 
weapon  hit  using  PDA,  APDA  and  BAPDA  compartment  designations 
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compartments  on  the  Second  Deck  directly  over  the  PDA  compartments  on  the  Third  Deck.  The 
BAPDA  compartments  are  either  horizontally  adjacent  to  or  above  or  below  APDA 
compartments. 

Figure  B-2  shows  how  changing  the  definition  of  the  PDA,  APDA  and  BAPDA  to  PDA(F), 
APDA(F)  and  BAPDA(F)  would  affect  the  definition  of  spaces.  Irrespective  of  the  terminology, 
vent  connected  spaces  have  always  been  included  (and  will  continue  to  be  included)  in  the  fire 
model  used  to  assess  the  spread  of  smoke  and  fire  associated  with  the  fires  that  start  immediately 
following  the  weapon  detonation.  As  a  result,  the  definitions  used  for  PDA  and  APDA  spaces 
will  not  impact  the  results  that  are  obtained  from  the  analysis.  From  this  point  forward,  the 
recommended  definitions  for  PDA(F),  APDA(F)  and  BAPDA(F)  are  used  in  this 
methodology  report. 


3.0  METHODOLOGY 


The  methodology  for  predicting  fire  and  smoke  spread  follows  the  steps  listed  below.  Each 
process  is  described  in  further  detail  in  the  following  sections. 

1 .  Establish  the  primary  damage; 

2.  Review  the  primary  damage  information; 

3.  Characterize  the  secondary  damage  (fire  and  smoke  spread)  in  the  PDA(F); 

4.  Characterize  the  secondary  damage  to  PDA(F)  spaces; 

5.  Characterize  the  secondary  damage  to  APDA(F)  and  BAPDA(F)  compartments;  and 

6.  Document  the  time  dependent  secondary  damage  estimates  in  snapshot  format. 
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First  Platform 


Fig.  B-2  -  Example  of  compartment  damage  characterization  after  hypothetical  weapon 
hit  using  PDA(F),  APDA(F)  and  BAPDA(F)  compartment  designations 


B-11 


3.1  Establishment  of  the  PDA(F) 


The  primary  damage  from  the  weapon  hit  is  established  based  on  output  from  either  the  Ship 
Survivability  Model  (SVM)  or  the  Advanced  Survivability  Assessment  Program  (ASAP)  model 
[8,  9].  Both  models  are  designed  to  simulate  initial  ship  damage  from  a  weapons  attack.  The 
output  from  either  model  provides  information  such  as  the  locations  of  failed  bulkheads  and 
decks,  shock  holing,  failed  doors  and  hatches,  overpressures,  bulkhead  and  deck  deflection, 
fragment  damage  to  vital  components  (VCs)  and  damage  to  distributive  systems.  Equipment 
specific  damage  as  a  result  of  thermal  or  smoke  exposures  during  the  detonation  and  burning  of 
the  residual  fuel  from  the  weapon  is  not  specifically  included  in  this  methodology. 

The  output  from  the  ASAP  or  SVM  model  is  used  to  establish  the  location  of  PDA(F), 
APDA(F)  and  BAPDA(F)  compartments.  Since  the  blast  damage  models  do  not  provide  the 
output  in  terms  of  PDA,  APDA  and  BAPDA  or  PDA(F),  PADA(F)  and  BAPDA(F),  the 
characterization  of  space  type  is  done  as  part  of  the  fire  and  smoke  damage  assessment.  Based 
on  the  model  output  provided,  PDA(F)  compartments  consist  of  those  compartments  where 
bulkheads  or  decks  have  failed,  or  where  shock  holing  has  occurred.  In  addition,  PDA(F) 
spaces  inelude  those  compartments  that  are  vent-connected.  Vent  connections  may  consist  of 
failed  doors  or  hatches,  or  openings  created  by  shock  holing.  The  structural  damage  incurred 
from  a  weapon  detonation  is  translated  into  ventilation  opening  sizes  based  on  advice  from 
NSWC/CD  Code  665.  For  shock  holing,  it  is  assumed  that  the  area  can  be  treated  as  a  single 
opening  with  an  area  equivalent  to  100%  of  the  shock  holing  value. 

In  addition,  it  is  assumed  that  the  fire  insulation  attached  to  boundaries  that  experience  less 
than  a  specific  amount  of  deflection  (as  reported  in  Reference  2)  remains  intact  and  performs  to 
specifications.  For  deflections  greater  than  this  value,  the  boundary  is  modeled  as  a  bare  steel 
bulkhead  or  deck. 


3.2  Review  of  Primary  Damage  Information 


The  assumptions  for  each  assessment  consist  of  general  assumptions  that  apply  for  each 
analysis  and  those  that  are  specific  to  the  hit  (shot  line).  General  assumptions  include  the 
following: 

1 .  The  ship  is  at  general  quarters  prior  to  each  hit  (e.g.,  material  condition  ZEBRA  is  set 
throughout  the  ship).  This  assumption  implies  that  all  ZEBRA  doors,  hatches  and 
fittings  are  closed. 

2.  Unless  otherwise  stated,  heating,  ventilation  and  air  conditioning  systems  (HVAC)  are 
not  considered.  Engineering  judgment  can  be  used  in  specific  instances  to  account  for 
the  effects  associated  with  failure  of  ventilation  ducts.  In  some  cases,  such  as  with 
smoke  ejection  systems,  these  HVAC  systems  may  help  mitigate  smoke  spread  or 
accumulation; 

3.  For  compartments  that  contain  combustible  materials,  the  fuel  load  is  considered  as 
sufficient  to  support  post  flashover  fires  during  the  evaluation; 
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4.  Thermally  thin  combustible  materials  (i.e.,  readily  ignitable)  are  located  such  that  fire 
spread  may  occur  across  an  intact  boundary.  These  materials  may  normally  be  in  this 
arrangement  or  may  be  redistributed  against  bulkheads  and  decks  as  a  result  of  the 
shock  impact  or  over  pressure; 

5.  Transient  opening  and  closing  of  doors  and  hatches  between  compartments  is  not 
considered  (from  people  exiting  the  area  or  investigators  performing  an  assessment); 

6.  The  weapon  hit  precipitates  ignition  in  all  PDA(F)  spaces  that  are  intimate  with 
compartment  containing  the  point  of  detonation  and  that  have  a  combustible  fuel  load 
(due  to  unspent  propellant  or  hot  fragments);  a  fire  growth  rate  curve  is  followed  for  the 
specific  fuel  load.  The  oxygen  concentration  in  the  space  could  have  considerable 
affect  on  the  fire  growth  curve; 

7.  Toxicity  issues  related  to  unspent  missile  propellant  are  not  considered; 

8.  Generally,  firefighting  activities  and  fire  mitigation  effects  are  not  considered; 

9.  For  compartments  that  are  not  included  in  the  fire  model  (APDA(F)  and  BAPDA(F) 
spaces),  smoke  production  or  accumulation  is  not  characterized.  Rather,  it  is  assumed 
that  these  spaces  become  untenable  for  unprotected  personnel  once  fire  spreads  to 
them;  and 

1 0.  The  fuel  load  and  air  supply  in  APDA(F)  and  BAPDA(F)  compartments  are  not 
eonsidered.  Rather,  it  is  assumed  that  there  is  a  sufficient  quantity  of  fuel  and  air  to 
support  flashover  conditions.  In  reality,  these  compartments  would  not  necessarily 
become  fully  involved,  particularly  if  they  are  elosed  (e.g.,  limited  ventilation). 

An  example  of  a  shot  specific  assumption  is  one  where  a  ventilation  system  remains 
operational  after  the  weapon  hit.  For  example,  this  situation  could  occur  if  a  well  deck  and/or 
vital  spaces  are  part  of  the  PDA(F).  Including  the  effects  of  forced  ventilation  systems  could 
result  in  an  increase  in  fire  intensity,  particularly  for  cases  where  the  fire  would  be  oxygen 
starved  otherwise. 


3.3  Characterization  of  Secondary  Damage  in  the  PDA(F) 

3.3.1  Development  of  the  Fire  Model 


A  fire  model  is  used  to  characterize  the  environment  in  PDA(F)  spaces  for  each  hit  scenario. 
APDA(F)  and  BAPDA(F)  spaces  are  not  ineluded  in  the  fire  model.  Currently,  the  multi¬ 
compartment  zone  fire  model  used  for  this  analysis  is  CFAST  version  3.1.7  [10,1 1]. 

In  developing  the  input  for  CFAST,  the  ventilation  conditions  must  be  defined.  These 
conditions  include  those  provided  by  natural  and  forced  ventilation.  Natural  ventilation  may  be 
provided  by  holes  in  the  skin  of  the  ship  (resulting  from  the  weapon  hit)  and  openings  between 
PDA(F)  compartments.  At  a  minimum,  oxygen  for  combustion  is  provided  naturally  through 
any  openings  to  weather  that  result  from  the  weapon  hit.  Forced  ventilation  may  be  provided  in 
vital  spaces  where  ventilation  is  not  secured  tmder  material  condition  ZEBRA. 
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Another  consideration  for  the  model  input  is  the  fuel  load  in  each  space.  The  types  of  fuel 
consist  of  unspent  missile  propellant,  Class  A  material  and  Class  B  fuel.  Examples  of  Class  A 
material  include  paper  products  found  in  storerooms  and  offices,  electrical  cable  insulation,  pipe 
insulation  and  bedding  in  berthing  spaces.  Class  B  fuels  may  be  introduced  by  JP-5  (or 
equivalent)  fuel  storage  tanks  that  are  breached  or  damage  to  fuel  lines.  The  following 
assumptions  are  used  for  characterizing  the  fires  in  these  spaces: 

The  fiiel  load  is  evenly  distributed  through  the  compartment; 

Trunks,  voids  and  enclosed  ladders  do  not  participate  in  fire  spread; 

In  general,  passageways  do  not  have  concentrated  fuel  loads  except  possibly  cabling 
and  pipe  insulation.  For  passageways  that  serve  as  companionways  or  ladderways, 
there  is  no  appreciable  fuel  load  and  fire  cannot  spread  to  or  from  the  space.  For  long 
passageways  (either  athwartship  or  fore-to-aft),  cabling  and  pipe  insulation  is  installed 
in  the  overhead.  Fire  cannot  spread  from  below  to  these  spaces  but  can  spread  from 
adjacent  compartments.  In  addition,  once  a  fire  is  established  in  the  overhead  of  these 
passageways,  fire  can  spread  to  adjacent  compartments  or  the  compartment(s)  above; 
and 

It  is  assumed  that  the  doors  to  washrooms,  water  closets  and  showers  (WRAVC/SH) 
that  serve  living  areas  are  open  when  a  fire  begins.  These  areas  will  become  involved 
in  the  fire  simultaneously  with  the  contiguous  space;  however,  any  resulting  fire  is 
likely  to  be  very  limited  given  the  low  fuel  loads  in  washroom  spaces.  Since 
WRAVC/SH  areas  typically  do  not  have  a  substantial  fuel  load  associated  with  them,  it 
is  assumed  that  they  do  not  propagate  fire  spread. 

For  the  initial  CFAST  run,  only  the  fuel  load  in  PDA(F)  compartments  is  considered.  It  is 
assumed  that  fires  ignite  in  all  PDA(F)  compartments  that  are  intimate  with  the  compartment 
containing  the  burst  point  and  have  available  fuel  at  the  time  of  the  weapon  detonation  .  This  is  a 
very  aggressive  assumption.  Since  specific  details  about  the  fuel  loading  are  not  usually 
available,  fire  growth  rates  associated  with  the  type  of  materials  expected  in  the  spaces  are  used 
for  these  typical  materials.  More  specifically,  it  is  assumed  that  the  fires  follow  “t  squared”  fire 
growth  curves.  This  type  of  fire  growth  curve  is  represented  by  Equation  1 : 


Q  =  ae  (1) 

where  Q  is  the  heat  release  rate  of  the  fire  (kW),  a  is  the  fire  growth  constant  typically  associated 
with  the  type  of  materials  in  the  space  (kW/s^)  and  t  is  the  length  of  time  the  fire  has  burned  (s). 

A  medium  fire  growth  rate  (a  =  0.01 17  kW/s^)  is  used  to  represent  Class  A  materials.  This 
value  was  adopted  after  analysis  of  potential  fire  threats  [12]  and  was  validated  during  a  weapons 
effects  test  (WET)  conducted  on  the  ex-USS  Dale  (CG-19)  [13].  An  ultra-fast  fire  growth  rate 
(a  =  0.1876  kW/s^)  is  used  to  represent  Class  B  pool  fires  involving  JP-5  or  equivalent  diesel  fuel 
[14]. 
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3.3.2  Interpretation  of  Fire  Modeling  Results 


The  results  from  the  fire  modeling  are  used  to  predict  the  temperatures,  interface  heights 
(demarcation  between  smoke  and  relatively  smoke  free  zones)  and  smoke  levels  in  PDA(F) 
compartments.  The  temperature  data  are  used  to  determine  when  or  if  fire  will  spread  to  these 
spaces.  The  interface  height  and  smoke  density  data,  in  addition  to  the  temperature  data,  are 
used  to  determine  how  tenable  these  compartments  would  be  for  protected  and  unprotected 
personnel.  The  criteria  used  for  tenability  are  described  in  Section  3.4.  This  information  is  the 
most  basic  output  provided  by  the  fire  model.  It  has  been  validated  and  verified  extensively 
[15-19]. 


3.4  Fire  and  Smoke  Spread  within  the  PDA(F) 


Fire  may  spread  to  compartments  within  the  PDA(F)  via  either  hot  gases  passing  through  an 
open  vent  or  heating  of  an  intact  boundary  (bulkhead  or  deck).  In  this  situation,  fires  resulting 
from  the  weapon  detonation  in  the  immediate  vicinity  of  the  burst  point  would  spread  to  other 
compartments  within  the  PDA(F).  Both  modes  of  fire  spread  are  examined  when  evaluating  fire 
spread  to  these  PDA(F)  compartments  to  determine  which  mode  results  in  a  faster  (more 
conservative)  fire  spread  time. 


3 .4. 1  Fire  Spread  Through  V ents 


Fire  spread  through  vents  that  connect  compartments  within  the  PDA(F)  may  result  from 
several  different  scenarios.  The  highest  probability  of  spreading  fire  occurs  when  a  PDA(F) 
compartment  reaches  flashover.  In  this  situation,  hot  gases  that  are  in  excess  of  500  C  (932®F) 
will  flow  into  the  adjacent  space  resulting  in  a  radiant  heat  flux  exposure  of  at  least  20  kW/m^. 
These  fluxes  exceed  the  critical  radiant  heat  flux  for  piloted  ignition  of  wood  and  select  plastic 
materials  [20].  The  configuration  or  orientation  of  the  fixel  in  vent  connected  spaces  will  not 
affect  how  readily  it  is  ignited. 

Another  ignition  scenario  occurs  when  the  fire  is  located  in  close  proximity  to  the  door,  hatch 
or  vent  opening.  The  radiation  from  the  flame  can  pass  through  the  opening  and  ignite 
combustible  materials  in  the  adjacent  space.  In  order  to  achieve  ignition,  the  exposure  must 
exceed  the  critical  ignition  heat  flux  associated  with  the  combustible  material.  Critical  ignition 
heat  fluxes  can  be  as  low  as  10-15  kWW  [21].  This  ignition  mechanism  is  difficult  to  apply  in  a 
typical  TSS/BDRA  since  it  is  dependent  on  the  relative  location  and  geometry  of  both  the  fire 
and  the  target  fuel  source. 

The  final  scenario  results  from  convective  heating  produced  by  hot  gases  (below  500‘’C 
(932  F))  flowing  into  the  adjacent  compartment.  These  hot  gases  can  potentially  heat 
combustibles  in  the  vent  connected  PDF(F)  compartments  to  their  auto-ignition  temperature. 

The  presence  of  glowing  embers  or  firebrands  may  cause  ignition  more  quickly.  The  auto¬ 
ignition  temperature  of  materials  can  be  measured  through  the  use  of  the  Setchkin  Furnace  [22]. 
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The  standardized  procedure  for  measuring  these  temperatures,  ASTM  D  1929  [23],  has  been 
used  for  a  wide  range  of  materials.  Review  of  published  test  data  show  that  paper  will  ignite  at 
approximately  230°C  (446"F)  and  various  polymers  will  ignite  between  330-450°C  (626-842°F) 
[24,  25]. 

Based  on  the  review  of  these  potential  ignition  scenarios,  it  is  conservatively  assumed  that 
ignition  via  hot  gases  will  occur  in  a  PDA(F)  compartment  when  the  accumulated  fire  gas 
temperatures  reach  230°C  (446  °F)  in  that  adjacent  space. 


3.4.2  Fire  Spread  through  Intact  Boundaries 


Fire  spread  within  the  PDA(F)  may  also  occur  via  heat  transfer  through  the  thermally  thin 
steel  bulkheads  [24,  26-28].  This  scenario  occurs  globally  during  post  flashover  fires.  It  could 
also  occur  locally  with  smaller  fires  if  the  burning  material  is  in  close  proximity  to  the  boundary. 
This  would  be  considered  fire  spread  across  an  unbreached  boundary. 

When  flames  impinge  directly  on  a  steel  boundary,  the  temperature  of  the  boundary  will 
increase  until  it  begins  to  glow  red  at  approximately  350°C  (572°F).  Depending  on  the  heat 
source,  the  temperatures  could  reach  values  as  high  as  650®C  (1202®F).  The  boundary  serves  as  a 
fairly  efficient  radiant  heat  source  that  will  increase  the  temperature  of  nearby  objects 
significantly  [20].  The  amount  of  energy  that  is  transferred  through  a  bare  steel  boundary 
resulting  from  direct  flame  impingement  can  approach  50  kW/m^  [20].  This  energy  consists  of 
both  convective  (-50%)  and  radiative  (-50%)  components.  As  shown  in  Figure  B-3,  the  worst 
case  heat  flux  exposmes  decrease  from  25  kW/m^  at  the  boundary  surface  to  approximately 
1.0  kW/m^  at  a  distance  of  ten  meters  away  from  the  boundary.  As  points  of  reference, 
spontaneous  ignition  of  combustible  materials  typically  occurs  around  20  kW/m^,  piloted 
ignition  of  combustibles  occurs  around  10  kWW  and  bare  skin  begins  to  blister  for  short 
exposures  (less  than  one  minute)  around  5  kW/m^  [21, 29, 30].  Under  worst  case  conditions, 
spontaneous  ignition  can  occur  as  far  away  as  0.8  m  (2.5  ft)  away  from  the  boundary. 
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Fig.  B-3  “  Worst  case  heat  flux  exposures  as  a  function  of  distance  from  a  bulkhead  or  deck 


For  aluminum  the  heat  transfer  would  be  difficult  and  failure  of  the  bulkhead  could  be 
expected  under  extreme  conditions.  The  melting  point  of  pure  aluminum  is  660°C  while 
aluminum  alloys  can  melt  around  500®C  [31].  Also  aluminum  loses  one  half  of  its  structural 
strength  at  330°C.  There  is  a  significant  difference  for  composites.  As  composites  are 
insulators,  the  heat  transfer  would  be  significantly  delayed.  It  should  be  kept  in  mind  that  most 
composites  in  general  use  today  could  become  fuels  with  time/temperature  insult.  This  would 
eventually  result  in  a  breached  boundary. 

Ignition  can  also  occur  when  combustible  materials  are  in  contact  with  the  hot  boundary  and 
the  temperature  exceeds  the  auto-ignition  temperature  of  the  material  (230°C  (446  °F)). 

Although  both  mechanisms  (contact  or  radiant  heating)  are  realistic  scenarios,  fire  spread  will 
occur  the  most  quickly  with  contact  heating.  However,  this  type  of  scenario  is  difficult  to  model 
without  knowledge  of  the  fuel  load  distribution.  Even  in  cases  where  the  normal  distribution  is 
known,  it  is  likely  that  items  will  be  relocated  as  a  result  of  the  impact  fi'om  the  weapon  hit.  As  a 
result,  it  was  assumed  that  smaller,  pre-flashover  fires  could  cause  ignition  of  ordinary 
combustibles  due  to  thermal  penetration  of  steel  boundaries  on  the  same  time  scale  as  post 
flashover  fires. 
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In  order  to  determine  the  amount  of  time,  Xspread,  that  is  required  to  spread  fire  across  a 
boundary  (e.g.,  heat  it  to  230°C  (446°F)),  Equation  2  is  used: 


tspread  fgrowth  tpenetration  (2) 

where  Xgrowth  is  the  time  for  the  fire  to  reach  flashover  in  the  initiating  compartment  and  Xpenetration 
is  the  time  required  to  transfer  heat  through  the  boundary.  All  times  are  measured  in  minutes. 

This  equation  applies  to  both  vertical  and  horizontal  spread  independent  of  the  fire  source 
(post  flashover  or  localized  heating). 

The  fire  growth  time  {Xgrowth)  is  a  function  of  the  fuel  loading  and  configuration  in  the 
compartment,  the  oxygen  available  for  combustion  (ventilation)  and  the  size  of  the  compartment. 
Since  the  fuel  loading  and  configuration  is  unknown  at  the  time  of  the  hit,  it  is  assumed  that  there 
is  adequate  fuel  to  sustain  a  fully  developed  fire  for  the  duration  of  this  simulation. 

The  fire  growth  times  were  determined  by  modeling  a  range  of  compartment  sizes  using 
CFAST.  The  compartment  footprints  {Aj)  ranged  from  9.3  to  464.5  m^  (100  to  5000  ft^).  The 
compartment  height  was  taken  as  3.0  m  (10  ft),  a  typical  deck-to-deck  height.  It  was  assumed 
that  the  fire  is  unconstrained  (i.e.,  ample  oxygen  to  support  combustion)  in  a  closed 
compartment.  In  addition,  a  Class  A  fire  growth  curve  was  used  (a  =  0.01 17  kW/sec^).  The 
Class  A  growth  curve  was  selected  since  most  shipboard  spaces  will  contain  Class  A  fuels  as 
opposed  to  Class  B  fuels.  The  growth  time  was  defined  as  the  number  of  minutes  required  for 
the  upper  layer  temperature  to  reach  500°C  (932°F). 

The  fire  growth  times  that  were  conservatively  developed  based  on  the  fire  modeling  results 
are  summarized  in  Table  B-1 .  It  should  be  noted  that  the  growth  time  for  vertical  spread  from 
any  compartment  and  horizontal  spread  from  any  PDA(F)  compartment  is  5  minutes,  regardless 
of  the  compartment  size.  Since  fires  can  quickly  expose  the  overhead  to  high  heat  fluxes, 
vertical  fire  spread  to  taller  spaces  (multi-deck)  may  be  handled  on  a  case-by-case  basis  (i.e.,  hit 
specific  assumption). 
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Table  B-1  -  Fire  Growth  Times  (Tgrowth)  For  A  Range  Of  Compartment  Footprints 


Spread  Direction 

Compartment  Footprint  (Aj)  [m^  (ft^)] 

Tgrowth  (minutes) 

Vertical 

any  compartment 

5 

Horizontal 

any  PDA(F)  compartment 

5 

Horizontal 

0  (0)  <  Af<  93  (1000) 

5 

Horizontal 

93  (1000)  <  ^/^<  139  (1500) 

7 

Horizontal 

139  (1500)  <  <186  (2000) 

9 

Horizontal 

186  (2000) <  Af  <232(2500) 

11 

Horizontal 

232  (2500)  <  <  279  (3000) 

13 

Horizontal 

279  (3000)  <  Af  <  325  (3500) 

15 

Horizontal 

325  (3500)  <  Af  <  372  (4000) 

17 

Horizontal 

372  (4000)  <  <  418  (4500) 

19 

Horizontal 

418  (4500)  <Af<  465  (5000) 

21 

Horizontal 

Af  >  465  (5000) 

23 

The  thermal  penetration  time  (tpe„etration)  was  determined  using  the  finite  difference  heat 
conduction  model  -  HEATING  (Version  7.3  [32]).  Steel  plate  thicknesses  between  12.7  mm 
(0.5  in.)  and  101.6  mm  (4.0  in.)  were  modeled.  The  results  in  Table  2  are  based  on  the  following 
conservative  botmdary  conditions  applied  to  the  exposed  steel  surface:  the  surface  temperature 
was  set  at  1000°C  (1832^),  an  emissivity  of  0.7  was  used  and  a  convective  heat  transfer 
coefficient  of  10  W/m^K  was  imposed.  Based  on  numerous  tests  on  Shadwell ,  a  steel 
temperature  of  600“C  is  rare  as  is  a  flame  temperature  above  800°C  in  ship  compartments.  The 
penetration  time  was  defined  as  the  number  of  minutes  required  to  heat  the  unexposed  surface  to 
230°C  (446“F)  by  direct  flame  impingement  on  the  opposite  side.  The  Verification  and 
Validation  (V&V)  of  HEATING  has  been  accomplished  previously  and  is  being  summarized  as 
part  of  the  V&V  process  of  another  ship  platform  VAR.  Results  fi'om  this  analysis  are 
summarized  in  Table  B-2. 


Table  B-2  -  Penetration  times  (tpeneiration)  for  various  plate  thicknesses 


Plate  Thickness  (Pt)  [mm  (in.)] 

^penetration  (ininutes) 

P,<  12.6  (0.49) 

1 

12.7  (0.5)  <Pt<  19.0  (0.74) 

2 

19.1  (0.75)  <Pt<  25.3  (0.99) 

3 

25.4(1.0)  <Pt<  31.7(1.24) 

4 

5 

50.8  (2.0)  <P,<  76.1  (2.99) 

7 

76.2  (3.0)  <  Pt  <101.5(3.99) 

11 

P,  =101.6(4.0) 

15 
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It  should  be  recognized  that  the  times  shown  in  Table  2  are  conservative/generic  guidelines. 
They  do  not  account  for  the  effects  of  insulation  or  bulkhead/deck  treatments  which  could 
significantly  impede  thermal  penetration  or  prevent  it  altogether. 

In  summary,  the  amount  of  time  required  to  spread  fire  through  an  intact  boundary  is 
determined  by  identifying  the  fire  growth  time  and  the  thermal  penetration  time.  The  fire  growth 
time  is  a  function  of  the  type  of  compartment  (PDA(F)  as  compared  to  APDA(F)  or  BAPDA(F)) 
and/or  compartment  size  (Table  1).  The  thermal  penetration  time  is  a  function  of  the  bulkhead  or 
deck  thickness  (Table  2). 


3.4.3  Assessment  of  Tenability 


The  tenability  guidelines  are  divided  into  two  categories:  those  for  unprotected  personnel 
and  those  for  protected  personnel  (DC  parties).  It  is  assumed  that  protected  personnel  are 
outfitted  with  self-contained  breathing  apparatus  (SCBAs)  and  firefighting  ensembles  (FFEs). 

For  unprotected  personnel,  conditions  are  considered  untenable  when  the  smoke  is  too  hot  to 
breathe  or  the  visibility  is  too  poor  for  evacuation  activities.  Untenable  conditions  are  assumed 
to  occur  when  a  smoke  layer  is  below  1.5  m  (4.9  ft),  combined  with  one  or  both  of  the  following 
conditions  [33]: 

Upper  layer  temperature  >  100°C  (212  ®F)  and/or 

Visibility  <  4  m  (13  ft)  for  transported  troops  and  <  1.7  m  (5.6  ft)  for  unprotected  ship 
personnel. 

The  smoke  layer  height  of  1.5  m  (5.0  ft)  is  significant  because  it  is  below  the  normal  head 
height.  This  condition  would  result  in  the  potential  exposure  of  standing  personnel  to  hot  gases 
or  reduced  visibility.  The  temperature  criteria  was  chosen  based  on  the  fact  that  personnel 
exposed  to  gases  above  100°C  (212°F)  will  experience  pain  and  potential  damage  to  their 
respiratory  tract  [34-36].  Burning  of  the  respiratory  tract  can  result  in  severe  injury  or 
incapacitation  for  any  person  forced  to  breathe  these  gases.  The  failure  criteria  of  100°C  is 
conservative  as  compared  to  the  human  tolerance  temperatures  in  the  NSTM  555,  which  indicate 
incapacitation  after  an  exposure  for  5  minutes  to  gases  of  149°C  (300°F)  and  incapacitation  after 
an  exposure  for  35  minutes  to  gases  of  93°C  (199°F)  [21]. 

Visibility  throughout  the  primary  and  secondary  damage  compartments  is  important  from  a 
life  safety/egress  standpoint.  Successful  egress  is  contingent  upon  personnel  being  able  to  locate 
the  ladders,  hatches  and  the  appropriate  egress  paths.  Visibility  degradation  relates  to  the  ability 
of  occupants  to  egress  through  areas  that  may  be  smoke  logged.  It  is  possible  for  personnel  to 
become  disoriented  in  the  event  the  visibility  becomes  significantly  reduced.  The  minimum 
visibility  for  egress  activities  is  dependent  upon  the  occupant.  Since  some  ships  may  transport 
troops  that  are  not  familiar  with  the  ship’s  layout,  more  visibility  would  be  required  for  egress 
activities  than  for  ship  personnel.  This  difference  is  accounted  for  by  assuming  a  visibility 
criterion  for  ship  personnel  that  corresponds  to  the  more  generous  limit  identified  in  the  literature 
of  1.7  m  (5.6  ft)  [33].  This  value  corresponds  to  data  that  indicate  it  is  the  point  where  egressing 
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occupants  will  decrease  their  walking  speed  significantly  [37].  A  more  stringent  criterion  of 
4.0  m  (13  ft)  visibility  is  used  for  troops  that  are  being  transported  [33,  38,  39]. 

Failure  criteria  for  protected  personnel  were  established  based  on  the  eonditions  that  would 
not  permit  necessary  firefighting  operations  to  proceed.  A  smoke  temperature  of  200°C  (326°F) 
was  conservatively  chosen.  Certainly,  firefighters  in  FFEs  could  tolerate  higher  exposure 
temperatures  for  short  durations.  Tests  have  shown  that  a  firefighter  in  an  FFE  may  be  forced  to 
abandon  a  space  after  5  minutes  when  the  temperatures  are  approximately  260°C  (500  F)  at  head 
level.  This  translates  to  a  temperature  gradient  in  the  compartment  of  163°C  (325°F)  at  waist 
level  and  63°C  (145°F)  at  knee  level.  Well-trained  Navy  fire  fighters  that  have  been  through  the 
ex-US  S  Shadwell  experience  can  withstand  hot  gas  layers  above  these  and  at  higher  temperatures 
[40].  For  the  purpose  of  evaluating  the  fire  modeling  results,  it  was  assumed  that  a  layer 
interface  height  of  0.6  m  (2.0  ft)  and  an  upper  layer  temperature  of  200°C  (392°F)  constitute 
untenable  conditions  for  protected  personnel. 


3.5  Characterization  of  Secondary  Damage  to  APDA(r)  and  BAPDA(F)  Spaces 


Fire  spread  fi'om  APDA(F)  spaces  and  BAPDA(F)  spaces  occurs  via  heat  transfer  through 
intact  bulkheads  and  decks.  This  mode  of  fire  spread  is  the  most  predominant  since  these  spaces 
normally  do  not  have  vent  conneetions  to  the  spaces  with  fires.  As  a  result,  the  methodology 
described  in  3.4.2  is  used  to  characterize  secondary  damage  to  these  spaces.  The  tenability  of 
these  spaces  is  not  assessed.  Rather,  it  is  assumed  that  they  are  untenable  for  unprotected 
personnel  as  soon  as  fire  spreads  to  them.  The  analysis  continues  until  the  fire  spread  times  are 
larger  than  the  length  of  time  that  is  specified  for  the  analysis.  For  example,  if  secondary 
damage  is  only  a  concern  for  30  minutes  after  the  weapon  hit,  fire  spread  times  longer  than 
30  minutes  are  not  considered. 

In  many  situations,  the  use  of  these  prescriptive  rules  is  conservative  since  it  is  assumed  that 
there  is  an  unlimited  amount  of  fuel  and  air  available  in  each  compartment.  In  the  event  that 
there  is  a  small  amount  of  fuel  in  the  compartment  or  the  compartment  is  closed,  fire  spread  may 
be  delayed  or  may  not  occur  at  all. 


3.6  Documentation  of  Secondary  Damage  Estimates 


The  secondary  damage  for  a  particular  weapon  hit  is  reported  in  a  snapshot  format.  These 
snapshots  show  which  spaces  contain  fires  at  particular  points  in  time.  In  addition,  the  tenability 
for  PDA(F)  spaces  is  shown.  Fire  suppression  effects  are  shown  in  situations  where  suppression 
systems  are  considered.  Figures  B-4  thru  B-6b  show  examples  of  these  snapshots  5, 10  and  15 
minutes,  respectively,  after  the  hypothetical  weapon  hit  represented  in  Figures  B-1  and  B-2. 
Compartments  containing  fire  are  marked  in  red.  Compartments  that  are  not  tenable  for 
unprotected  personnel  are  designated  with  blue.  Those  compartments  that  are  untenable  for  all 
personnel  (unprotected  and  protected)  are  shown  in  yellow.  It  should  be  noted  that  the  rules- 
based  fire  spread  framework  was  used  for  the  development  of  these  examples.  These  results 
would  be  contingent  on  the  incorporation  of  CFAST  results. 
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Second  Deck 


Main  Deck 


Second  Deck 


Fig.  B-5a  -  Example  of  secondary  damage  10  minutes  after  hypothetical  weapon  hit 
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First  Platform 


Fig.  B-5b  -  Example  of  secondary  damage  10  minutes  after  hypothetical  weapon  hit  (continued) 


Main  Deck 


Second  Deck 


Fig.  B-6a  -  Example  of  secondary  damage  15  minutes  after  hypothetical  weapon  hit 
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Third  Deck 


First  Platform 


Fig.  B-6b  -  Example  of  secondary  damage  1 5  minutes  after  hypothetical  weapon  hit  (continued) 


The  methodology  is  summarized  in  the  flowchart  shown  in  Figure  B-7.  Figure  B-8  is  a 
flowchart  for  tenability.  The  process  involves  modeling  PDA(F)  compartments  with  the  multi¬ 
compartment  zone  fire  model  CFAST.  The  use  of  CFAST  allows  for  a  more  complete 
characterization  of  these  compartments.  The  output  from  CFAST  is  used  to  predict  fire  spread 
times  and  assess  tenability  conditions.  Fire  spread  beyond  these  compartments  to  APDA(F)  and 
BAPDA(F)  compartments  is  predicted  using  a  conservative  set  of  prescriptive  rules. 
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Establish  the  PDA  (F). 
APDA(F)  &  BAPDA(F)  (based 
on  SVM  or  ASAP  output) 
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Fig-  B-7  -  Flowchart  for  fire  and  smoke  spread  methodology 


Assess  tenability 
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4.0  SUMMARY 


A  methodology  has  been  developed  as  part  of  a  congressionally  mandated  LFT&E  programs 
to  predict  smoke  and  fire  spread  after  a  weapon  hit.  The  objective  of  the  current  methodology  is 
to  identify  specific  weaknesses  and/or  shortcomings  that  might  increase  survivability. 


5.0  FUTURE  DEVELOPMENTS 


The  methodology  described  in  this  report  is  currently  being  applied  in  the  VAR  3  process  for 
the  LPD  17  platform  detailed  design.  Any  further  developments  or  improvements  in  the 
methodology  will  be  documented  and  reflected  in  a  future  revision  of  this  report.  Separate 
efforts  are  underway  to  develop  a  network  fire  model  [41].  The  network  fire  model  would 
provide  a  physics  based  model  that  could  quickly  assess  fire  and  smoke  spread  for  ship 
platforms.  Plans  also  include  accounting  for  fire  detection  and  suppression  effects.  This  type  of 
model  would  provide  better  far  field  secondary  damage  characterization  capabilities.  It  is 
unclear  if  the  network  fire  model  will  replace  the  current  methodology  or  if  it  will  be  used  in 
conjimction  with  the  method  or  portions  of  the  method  described  in  this  report.  Significant 
changes  in  direction  for  predicting  secondary  damage  associated  with  LFT&E  efforts  will  be 
documented  in  future  reports. 
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1.0  INTRODUCTION 


A  series  of  eight  benchmark  reference  cases  were  proposed  by  Wickstrom  and  Palsson 
[1999]  for  use  in  verifying  two-dimensional  heat  conduction  computer  programs  [Wickstrom  et 
al.,  1999;  Palsson  et  al.,  2000].  These  entail  thermal  calculations  of  increasing  complexity.  The 
reference  cases  were  developed  to  address  the  following  aspects  of  a  numerical  heat  transfer 
calculation: 

•  Convergence; 

•  Variable  convective  and  radiative  boundary  conditions; 

•  Variable  thermal  material  properties; 

•  Multiple  materials; 

•  Phase  change  of  material  constituents;  and 

•  Radiation  and  convection  in  voids. 


Convergence  is  verified  by  starting  with  a  coarse  grid  in  each  domain  and  then  by  increasing 
the  number  of  nodes  until  a  desired  error  tolerance  is  met.  The  remaining  aspects  are  verified 
when  the  difference  between  the  converged  numerical  results  and  the  touchstone  calculations  is 
below  an  acceptable  error  threshold. 

In  this  study,  the  eight  reference  cases  will  be  solved  with  the  finite  volume  conduction 
program  HEATING  7.3.  HEATING  stands  for  heat  engineering  and  transfer  in  nine  geometries. 
HEATING  was  developed  at  the  Oak  Ridge  National  Laboratory's  (ORNL)  Radiation  Safety 
Information  Computational  Center  (RSICC)  [Childs,  1998].  HEATING  was  created  to  analyze 
the  thermal  impact  of  the  U.S.  Energy  Department's  various  high  energy  projects,  including  all 
aspects  of  radiation  shielding.  It  has  one  of  the  longest  development  histories  of  any 
computational  heat  transfer  software  package:  [Fowler  &  Volk,  1959;  Liguori  &  Stephenson, 
1961;  Turner  &  Crowell,  1969;  Turner  &  Siman-Tov,  1971;  Turner  et  al.,  1977;  Becker,  1977; 
Elrod  et  al.,  1984;  Bryan  et  al.,  1986;  Childs,  1990;  Childs  et  al.,  1990;  Childs,  1991]. 
Verification  studies  and  tools  for  HEATING  can  be  found  in  [Bryan  et  al.,  1986;  Chu,  1989]. 

HEATING  is  a  field  model.  It  solves,  using  finite  volume  techniques  [Patankar,  1980; 
Childs,  1998],  the  heat  equation  [Arpaci,  1966;  Carslaw  &  Jaeger,  1959], 

pc^  =  V-{kVT)+q\  (1) 

ot 
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subject  to  the  general  boundary  condition 


«  =  =  {a,  -T.(-+h\T,^  *t:It,+TJ](T,-TJ 

or  .  (2) 

T,=U^,t) 

The  variables  are  as  follows;  T  is  the  material  internal  temperature,  p  is  the  density,  c  is 
the  specific  heat,  t  is  the  time,  q"  is  the  volumetric  heat  generation  rate,  ql  is  the  heat  flux  for 
surface  ^ ,  n  is  the  coordinate  normal  to  the  surface,  k  is  the  thermal  conductivity,  is  the 
material  surface  temperature,  is  a  temperature  function  associated  with  surface  5 ,  7^  is  the 
ambient  temperature,  is  the  heat  transfer  coefficient  for  forced  convection,  is  the  heat 
transfer  coefficient  for  natural  convection,  \  is  the  exponent  for  natural  convection,  is  the 
radiative  transfer  coefficient,  and  V  and  d  are  partial  derivative  operators. 

HEATING  7.2i  and  7.3  are  the  most  recent  developments  in  a  series  of  heat-transfer  codes 
distributed  by  RSICC  as  SCA-1 /HEATINGS  and  PSR199/HEATING  6.  HEATING7  can  solve 
steady-state  and/or  transient  heat  conduction  problems  in  one-,  two-,  and  three-dimensional 
Cartesian,  cylindrical,  or  spherical  coordinates.  A  model  may  include  multiple  materials,  and  the 
thermal  conductivity,  density,  and  specific  heat  of  each  material  may  be  both  time-  and 
temperature-dependent.  The  thermal  conductivity  may  also  be  anisotropic.  Materials  may 
undergo  change  of  phase.  Thermal  properties  of  materials  may  be  input  or  may  be  extracted 
from  a  material  properties  library.  Heat-generation  rates  may  be  dependent  on  time,  temperature, 
and  position.  Boundary  temperatures  may  be  time-  and  position-dependent.  The  boundary 
conditions,  which  may  be  surface-to-environment  or  surface-to-surface,  may  be  specified 
temperatures  or  any  combination  of  prescribed  heat  flux,  forced  convection,  natural  convection, 
and  radiation.  The  boimdary  condition  parameters  may  be  time-  and/or  temperature-dependent. 
General  gray-body  radiation  problems  may  be  modeled  with  user-defined  factors  for  radiant 
exchange.  Surfaces  that  are  not  in  contact  can  exchange  heat  via  custom  connectors.  The  mesh 
spaeing  may  be  variable  along  each  axis. 


2.0  VERIFICATION 


[Wickstrom  &  PMsson,  1999]  published  a  suite  of  eight  configurations  suitable  for  the 
verification  of  heat  transfer  programs  that  can  predict  the  temperature  in  structures  exposed  to 
fires.  They  used  TASEF  [Sterner  and  Wickstrdm,  1990]  as  the  benchmark  for  verifying  other 
computer  codes.  TASEF  stands  for  "Temperature  Analysis  of  Structures  Exposed  to  Fire."  It 
employs  the  finite  element  method  in  order  to  analyze  heat  transfer  in  two-dimensional  Cartesian 
and  axisymmetric  cylindrical  coordinates.  These  region  may  be  composed  of  multiple  materials 
and  voids.  TASEF  is  general  enough  to  account  for  nonlinearities  in  the  boundary  conditions 
and  in  the  material  properties.  As  the  constituents  of  its  acronym  imply,  TASEF  was  specifically 
developed  to  determine  the  response  of  structural  assemblies  and  components  to  fire.  The  first 
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reference  case  in  the  suite  is  a  linear  problem  that  has  an  exact  analytic  solution.  The  following 
five  reference  cases  are  non-linear  problems  that  do  not  have  exact  solutions  available.  The  last 
two  examples  in  this  suite  deal  with  heat  transfer  in  enclosures.  The  results  produced  by 
HEATING  7.3  are  to  be  evaluated  for  these  eight  cases. 

The  computer  model  INST2D  [Rudolph!  and  Muller,  1980]  has  been  verified  using  the  first 
seven  benchmark  reference  cases  [Hamann  et  al.,  1999].  Hamann  et  al.  demonstrate  reasonable 
agreement  between  TASEF  and  INST2D.  SAFIR,  another  finite  element  heat  transfer  program, 
has  been  used  to  successfully  solve  Reference  Case  7  [Pintea  and  Franssen,  1997].  These  two 
favorable  outcomes  add  credibility  to  the  verification  suite  based  on  results  from  TASEF.  This 
document  summarizes  the  benchmark  verification  process  per  [Wickstrom  and  Palsson,  1999] 
for  the  finite  volume  heat  transfer  computer  code  HEATING,  version  7.3  [Childs,  1998]. 

It  is  worthwhile  to  take  the  opportunity  now  to  clarify  some  points  that  will  make  the 
following  comparison  between  TASEF  and  HEATING  more  sensible.  HEATING  is  a  finite 
volume  program  that  offers  first  and  second  order  accuracy  in  space  and  time,  depending  on  the 
algorithm  that  is  used.  Spatial  accuracy  can  only  be  improved  by  increasing  the  number  of  cells. 
TASEF  is  a  finite  element  program.  Spatial  accuracy  is  likewise  governed  by  the  number  of 
elements  but  it  is  also  affected  by  the  order  of  the  multidimensional  polynomial  serving  as  the 
basis  function  for  each  finite  element.  The  geometrical  shape  of  the  subvolumes  for  each  of  the 
two  methods  also  matters.  Hence  a  finite  element  model  composed  of  20  fourth  order 
multidimensional  polynomial  quadrilaterals  is  expected  to  be  more  accurate  than  a  second  order 
finite  volume  simulation  with  20  cells.  When  percent  differences  are  presented,  the  “actual”  is 
taken  to  be  the  solution  from  an  analytic  expression,  if  available.  Otherwise,  the  answer 
provided  by  TASEF  is  employed.  All  the  specifications  for  the  examples,  unless  otherwise 
noted,  are  taken  fi'om  the  publications  [Wickstrom  and  Palsson,  1999;  Wickstrom,  1999;  Palsson 
&  WickstrOm,  2000].  The  HEATING  calculations  were  done  on  a  Linux  workstation  with  two 
Pentium  III  Xeon  processors,  1  GByte  of  memory,  and  19  Gbytes  of  storage  and  a  Silicon 
Graphics,  Inc.  Octane2  graphics  workstation  with  V6®  graphics,  two  R 12,000 A  processors, 

2  Gbytes  of  RAM  memory,  and  in  excess  of  136  Gbytes  of  storage.  The  October  1999  release  of 
HEATING  7.3  was  used  for  this  investigation. 


2.1  Reference  Case  1:  Comparison  Against  Analytical  Results,  Constant  Material 
Properties 


Reference  Case  1  consists  of  a  square  section  of  a  generic  material  whose  thermal  properties 
are  to  be  determined  from  the  stipulation  that  its  thermal  diffusivity  (a  =A:/yO  c)  and  Biot 
number  (Bi=A£/yfc)are  both  unity.  L  is  the  half-length  of  the  side  of  the  square  section.  The 
material,  initially  at  0°C,  is  subjected  to  a  surface  gas  temperature  of  1,000°C.  Convection  is  the 
means  of  heat  transfer  between  the  square  section  and  the  surrounding  gas.  The  specific  material 
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properties  are  arbitrary  in  this  case.  Refer  to  Figure  C-1  for  a  depiction  of  the  calculation 
domain.  The  implicit  Crank-Nicolson  method  [Childs,  1998;  Smith,  1987;  Jaluria  &  Torrance, 
1986]  is  used  to  solve  for  the  temperature  within  the  square.  The  tolerance  for  convergence  is 
10-^. 

The  temperature  at  the  center  of  the  square  region  is  presented  as  a  function  of  dimensionless 
time  (i.e.,  the  Fourier  number).  The  Fourier  number,  Fo,  is  at!  Il .  Figure  C-2  shows  the 
results  of  HEATING  and  TASEF  for  grids  using  4  through  256  elements.  This  plot  shows  that 
both  TASEF  and  HEATING  converge  to  a  solution  with  fewer  than  64  elements.  Table  C-1 
compares  HEATING  with  the  exact  solution  for  a  64  cell  grid.  Table  C-1  demonstrates  excellent 
agreement  with  the  results  of  the  analytical  solution  reported  in  [Wiclcstrdm  and  PMsson,  1999] 
which  is  based  on  the  formulae  in  [Carslaw  &.  Jaeger,  1959]. 


Fig.  C-1  -  Diagram  of  benchmark  Reference  Case  1. 
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Symbols  -  TASEF  (Benchmark  Model)  Lines  -  HEATING? 


Fig.  C-2  -  Center  temperature  for  Reference  Case  1  as  a  function  of  the  number  of  grid  cells. 


Table  C-1  -  Comparison  of  HEATING  to  Analytical  Solution  for  the  Reference  Case  1. 


Fourier  Number 

Center  Temperature  f°Cl  I 

HEATING  (64  Cells) 

Analytical  Solution 

%  Difference  | 

0.1 

986 

986 

0.2 

904 

904 

0.4 

691 

690 

515 

515 

384 

383 

285 

285 

0.000  1 

2.2  Reference  Case  2:  Non-Linear  Boundary  Conditions  and  Constant  Material 
Properties 


Reference  Case  2  consists  of  a  square  section  of  a  concrete-like  material  with  constant 
properties.  Figure  C-3  shows  the  region.  The  dimensions  and  material  properties  are  as  follows: 


•  Length:  0.1  m 

•  Thermal  conductivity:  1.0W/(m-K) 

•  Heat  capacity:  1,000  J/kg 

•  Density:  2,400  kg/m^. 
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•  Density: 


2,400  kg/m\ 


Fig.  C-3  —  Diagram  of  benchmark  Reference  Case  2 

The  boundary  condition  is  applied  in  the  following  form: 

q-  =  as(T;-n)+h(T.-T.)  (3) 

where  q"  is  the  boundary  heat  flux  (W/m^),  e  is  the  surface  emissivity,  a  is  the  Stefan- 
Boltzman  constant  (5.669x10'^  wW  k'’),  is  the  exposure  temperature  (K),  is  the  surface 
temperature  (K),  and  h  is  the  convective  heat  transfer  coefficient  (W/m^-K).  The  boundary 
condition  parameters  are  as  follows: 

•  Emissivity: 

•  Convection  coefficient: 

•  Initial  temperature  (Tg): 

•  Exposure  temperature  ( 1 ) : 

•  Exposure  temperature  (2): 


0.8 

10  W/(m2-K) 

0°C 
1,000  °C 

Tg  +  345  logic  (8r  +  0  >  where  t  is  in  minutes. 
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The  implicit  Crank-Nicolson  scheme  [Childs,  1998;  Smith,  1987;  Jaluria  &  Torrance,  1986] 
is  used.  The  convergence  tolerance  is  10'®.  Figures  C-4  and  C-5  compare  the  temperature 
predictions  of  TASEF  and  HEATING  for  the  three  locations  shown  in  Figure  C-3;  the  center,  the 
surface,  and  the  comer.  The  time  window  within  which  the  temperatures  are  compared  is 
30  minutes  to  180  minutes.  The  figures  show  that  HEATING  and  TASEF  are  in  agreement 
when  64  or  more  cells  are  used  for  the  reference  case. 


Fig.  C-4  —  HEATING  and  TASEF  temperature  predictions  for  Reference 
Case  2  -  constant  1 ,000  °C  exposure  temperature. 
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Fig.  C-5  —  HEATING  and  TASEF  temperature  predictions  for  Reference 
Case  2  -  transient  exposure  temperature. 

2.3  Reference  Case  3:  Non-Linear  Boundary  Conditions  and  Temperature  Dependent 
Material  Properties 


Reference  Case  3  is  the  same  as  Reference  Case  2  except  that  the  thermal  conductivity  of  the 
square  region  is  temperature-dependent  as  follows  in  Table  C-2. 

Table  C-2  -  Time-dependant  thermal  conductivity  of  the  square  region. 


Temperature  [°C] 

Thermal  Conductivity  [W/(m-K)] 

0 

1.5 

200 

0.7 

1000 

0.5 

(HEATING  linearly  interpolates  the  data  in  Table  C-2.)  Figures  C-6  and  C-7  compare  the 
temperature  predictions  of  TASEF  and  HEATING  for  the  three  locations  indicated  in  Figure  3. 
The  time  interval  when  the  temperatures  are  compared  is  30  minutes  to  180  minutes.  The  figures 
show  that  HEATING  and  TASEF  are  in  agreement  when  64  or  more  cells  are  used. 
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2.4  Reference  Case  4:  Latent  Heat  Due  to  Water  Content  -  A  Concrete  Block  with 
Moisture 


Reference  Case  4  is  the  same  as  Reference  Case  3  except  that  a  five  percent  by  weight 
moisture  content  is  assumed  in  the  material.  The  change  of  phase  of  the  liquid  water  within  the 
concrete-like  material  alters  the  thermal  response  because  the  latent  heat  of  vaporization  of  the 
water  acts  as  a  heat  sink  for  the  whole  composite.  This  effect  is  modeled  by  (discontinuously) 
stepping  the  specific  heat  c  up  from  5650  J/(kg-K)  to  6858  J/(kg-K)  in  the  temperature  range  of 
100°C  to  120°C.  The  implicit  Crank-Nicolson  scheme  [Childs,  1998;  Smith,  1987;  Jaluria  & 
Torrance,  1986]  is  used.  The  convergence  tolerance  is  10'^. 

Figures  C-8  and  C-9  compare  the  temperature  predictions  of  TASEF  with  HEATING  for  the 
three  locations  shown  in  Figure  C-3.  The  times  at  which  the  temperatures  are  compared  vary 
from  30  minutes  to  180  minutes.  The  figures  show  that  HEATING  and  TASEF  are  in  agreement 
when  64  or  more  cells  are  used  for  the  reference  case. 


2.5  Reference  Case  5:  Composite,  Steel  and  Concrete 


Reference  Case  5  is  the  same  as  Referenee  Case  4  except  that  a  10  mm  steel-like  tube 
surrounds  the  concrete-like  core  material.  The  properties  of  the  steel-like  material  are  as 


follows: 

•  Thermal  conductivity  @  20°C: 

54  W/(m  K) 

•  Thermal  conductivity  @  800°C: 

27  W/(m-K) 

•  Thermal  conductivity  @  1,000°C: 

27  W/(m-K) 

•  Heat  capacity: 

600  J/(kgK) 

•  Density: 

7,850  kg/m3 

Reference  Case  5  is  shown  in  Figure  C-10.  The  implicit  Crank-Nicolson  scheme  [Childs, 
1998;  Smith,  1987;  Jaluria  &  Torrance,  1986]  is  used.  The  convergence  tolerance  is  5x10'^. 
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Temperature  Predicted  by  HEATING  (“C)  Temperature  Predicted  by  HEATING  (“C) 


Fig.  C-6  —  HEATING  and  TASEF  temperature  predictions  for  Reference 
Case  3  -  constant  1,000  exposure  temperature. 


Fig.  C-7  —  HEATING  and  TASEF  temperature  predictions  for  Reference 
Case  3  -  transient  exposure  temperature. 
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Temperature  Predicted  by  HEATING  (*C)  Temperature  Predicted  by  HEATING  (*C) 


Temperature  Predicted  by  TASEF  (®C) 


Fig.  C-8  —  HEATING  and  TASEF  temperature  predictions  for  Reference 
Case  4  -  constant  1,000  °C  exposure  temperature. 


Temperature  Predicted  by  TASEF  (®C) 


Fig.  C-9  —  HEATING  and  TASEF  temperature  predictions  for  Reference 
Case  4  -  transient  exposure  temperature. 
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Fig.  C-10  -  Diagram  of  benchmark  Reference  Case  5. 


Figures  C-1 1  and  C-12  compare  the  temperature  predictions  of  TASEF  and  HEATING  for 
the  three  locations  shown  in  Figure  C-IO  in  the  time  window  from  30  minutes  to  180  minutes. 
The  figures  show  that  HEATING  and  TASEF  are  in  agreement  when  64  or  more  cells  are  used. 


Fig.  C-1 1  —  HEATING  and  TASEF  temperature  predictions  for  Reference 
Case  5  -  constant  1,000  °C  exposure  temperature. 
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Fig.  C-12  —  HEATING  and  TASEF  temperature  predictions  for  Reference 
Case  5  -  transient  exposure  temperature. 


2.6  Reference  Case  6:  Composite,  Steel  and  Mineral  Wool 


Reference  Case  6  is  the  same  as  Reference  Case  5  except  that  the  concrete-like  core  is 
replaced  with  a  mineral  wool-like  material.  The  properties  of  the  mineral  wool-like  material  are 
as  follows: 

•  Thermal  conductivity:  0.05  W/(m  K) 

•  Heat  capacity:  1,000  J/(kg-K) 

•  Density:  50  kg/m3 

Reference  Case  6  is  shown  in  Figure  C-13.  The  implicit  Crank-Nicolson  method  [Childs, 
1998;  Smith,  1987;  Jaluria  &  Torrance,  1986]  is  used.  The  convergence  tolerance  is  5x10"^. 
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Fig.  C-13  —  Diagram  of  benchmark  Reference  Case  6. 


Figures  C-14  and  C-15  compare  the  temperature  predictions  of  TASEF  with  HEATING  for 
the  three  locations  shown  in  Figure  C-13  in  the  time  interval  from  15  minutes  to  180  minutes. 
The  figures  show  that  HEATING  and  TASEF  are  in  agreement  when  64  or  more  cells  are  used 
for  the  reference  case. 


0  200  400  600  800  1000 


Temperature  Predicted  by  TASEF  (®C) 

Fig.  C-14  —  HEATING  and  TASEF  temperature  predictions  for  Reference 
Case  6  -  constant  1 ,000  °C  exposure  temperature. 


Fig.  C-15  —  HEATING  and  TASEF  temperature  predictions  for  Reference 
Case  6  -  transient  exposure  temperature. 


2.7  Reference  Case  7:  Radiative  Transfer  Across  an  Oblong  Rectangular  Cavity 


These  examples  demonstrate  the  effect  of  radiation  in  an  enclosure.  HEATING  needs  the 
results  of  external  calculations  in  order  to  determine  the  transfer  factor.  Appendix  A  describes 
the  procedure  used  to  obtain  the  transfer  factor.  The  domain  consists  of  a  rectangular  cavity  that 
is  ten  times  longer  in  one  dimension  than  the  other.  The  rectangular  enclosure  has  a  wall 
thickness  of  5  mm.  The  dimensions  of  the  cavity  are  100  mm  in  the  x -direction  and  10  mm  in 
the  y .  (See  Figure  C-16).  If  a  temperature  difference  is  applied  across  the  10  mm  gap,  the 
transfer  is  essentially  one-dimensional.  When  the  temperature  difference  is  applied  across  the 
100  mm  gap,  the  transfer  is  two-dimensional.  The  thermal  properties  of  the  case  are  k  =  1 
W/(m-K),  c  =  1000  J/(kg'K),  p  =  1000  kg/m^,  and  e  =  0.8.  A  steady-state  solution  is  sought. 

The  temperature  used  to  initiate  the  iterations  is  500°C. 
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Fig.  C-16  —  Contours  of  the  steady-state  temperature  (in  °C)  for  Reference  Case  7a, 


The  tolerance  for  convergence  is  lO"^.  The  solution  methods  that  have  been  used  are  the 
iterative  successive  over-relaxation  (SOR)  and  conjugate  gradient  procedures  as  well  as  direction 
solution  of  the  linear  algebraic  system  of  equations  using  a  numerical  procedure.  All  three  were 
found  to  yield  results  that  were  identical  to  three  significant  figures  or  more.  The  results  of  the 
conjugate  gradient  solution  are  presented. 


2.7. 1  Temperature  Difference  Applied  Across  the  y-Coordinate 


For  this  case  (7a),  the  bottom  temperature  is  1000°C,  the  top  temperature  is  0°C,  and  the  left 
and  right  sides  are  adiabatic.  These  sides  refer  to  the  outer  surface  of  the  solid  enclosure.  The 
results  that  are  presented  in  Table  C-3  clearly  show  (rapid)  convergence  and  excellent  agreement 
with  the  results  of  the  analytical  solution  reported  in  [Wickstrdm  and  Palsson,  1999].  The 
contour  plot  presented  in  Figure  16  confirms  that  the  heat  transfer  is  essentially  one-dimensional. 

Table  C-3  -  Temperature  [°C]  at  two  different  locations  for  heat  transfer  along 
the  narrow  gap  (Reference  Case  7a). 


Number  of  Bottom  Cells 

TASEF 

Analytical 

Lower,  Center 

2 

789 

779.4 

779.6 

0.0257 

4 

774 

779.6 

0.0257 

8 

779 

779.6 

0.0257 

2 

219.4 

Upper,  Center 

4 

220.6 

220.0 

8 

221 

220.0 
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2.7.2  Temperature  Difference  Applied  Across  the  ^-Coordinate 


For  this  case  (7b),  the  left  temperature  is  1000°C,  the  right  temperature  is  0°C,  and  the  top 
and  bottom  sides  are  adiabatic.  The  results  are  plotted  in  Figure  C-17.  Table  C-4  compares  the 
HEATING  results  with  those  of  TASEF  as  well  as  SAFIR  [Pintea  and  Franssen,  1997].  The 
comparison  is  plotted  in  Figure  C-18. 


Fig.  C-17  —  Contours  of  the  steady-state  temperature  (in  °C)  for  Reference  Case  7b. 


Table  C-4  -  Temperature  [°C]  at  several  j:-locations  (in  mm)  for  heat  transfer  across  the  wide  gap 
(Reference  Case  7b).  n  refers  to  the  number  of  bottom  cells. 


MODEL 

n 

X  [mm] 

O 

12.5 

25 

37.5 

50 

62.5 

75 

87.5 

TASEF 

2 

881 

629 

nBi 

HEATING 

860 

663 

lEQ 

%  Difference 

2.38 

5.41 

19.7 

TASEF 

B 

898 

748 

675 

499 

98 

HEATING 

891 

785 

672 

503 

107 

%  Difference 

0.780 

4.95 

0.444 

0.802 

9.18 

SAFIR 

8 

QQ 

820 

777 

662 

586 

484 

344 

90 

TASEF 

820 

779 

662 

587 

485 

335 

89 

HEATING 

904 

837 

785 

730 

667 

591 

490 

340 

QQ 

%  Difference 

0.221 

2.07 

0.777 

0.968 

0.755 

0.681 

1.03 

1.49 

QQI 
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Fig.  C-1 8  —  Comparison  of  the  results  for  Reference  Case  7b  at  the  positions  indicated  in 
Table  4.  n  represents  the  number  of  bottom  cells. 


2.8  Reference  Case  8:  Exposure  of  an  HE200B  Beam  Section 


This  example  demonstrates  the  effects  of  radiative  and  convective  heat  transfer  in  cavities. 
Appendix  A  describes  the  steps  taken  to  accomplish  this  with  HEATING.  For  this  case,  an 
HE200B  beam  section  is  protected  by  a  10  mm  thick  rectangle  of  Promatek  insulation  boards. 
The  outer  dimensions  of  the  section  are  200  mm  x  200  mm.  The  thickness  of  the  web  is  9  mm 
and  that  of  the  flange  is  15  mm.  Four  cavities  are  formed  when  the  beam  is  thus  enclosed.  The 
one  above  the  top  flange  has  a  10  mm  gap  in  the  y -direction.  The  other  dimension  is  100  mm 
long.  The  gap  below  the  bottom  flange  has  identical  dimensions.  The  two  remaining  gaps 
measure  45.5  mm  by  70  mm.  (See  Figure  C-19.)  The  exterior  of  the  Promatek  is  exposed  to  an 
ISO  834  fire  curve: 


'^isom  ~  ^0  +  3451og,o(8/  +  l),  (4) 


C-23 


Fig.  C-19  —  Diagram  of  benchmark  Reference  Case  8. 


where  t  is  in  minutes  and  the  temperature  is  in  °C  (or  K).  The  initial  temperature,  ,  is  20°C. 

For  the  outside  of  the  the  assembly,  e  =  0.8  and  h  =\0  W/(m^  K).  The  thermal  properties  for 
the  steel  are  c  =  600  J/(kg  K),  p  =  7850  kg/m^,  s  =  0.8,  and  the  thermal  conductivity  is  given 
in  Table  C-5. 


Table  C-5  -  Thermal  conductivity  of  the  steel  for  Reference  Case  8. 


Temperature  T 

rc] 

Thermal  Conductivity  k 
[W/(m-K)] 

20 

54.0 

800 

27.3 

1200 

27.3 

The  thermal  properties  of  the  Promatek  are  c  =  1 130  J/(kg-K),  p  -  870  kg/m^,  s  -  0.8,  and 
the  thermal  conductivity  is  given  in  Table  C-6.  As  was  stated  previously,  radiative  and 
convective  heat  transfer  can  occur  within  the  four  gaps.  The  heat  transfer  coefficient  for  the  top 
and  bottom  gaps  is  1 .5  W/(m^-K)  and  it  is  2.0  W/(m^-K)  for  the  left  and  right  cavities.  Four 
different  meshes  are  used  for  the  HEATING  calculations.  (See  Figure  C-20.)  Although 
HEATING  can  match  the  grids  used  by  TASEF,  for  the  sake  of  simplicity  these  meshes  are  more 
uniform  than  those  employed  by  [Wickstrom  &  Palsson,  1999].  A  transient  solution  is  sought. 
The  Crank-Nicolson  implicit  algorithm  [Childs,  1998;  Smith,  1987;  Jaluria  &  Torrance,  1986]  is 
used.  The  tolerance  for  convergence  is  lO'^. 
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Table  C-6  -  Thermal  conductivity  the  Promatek  insulation  boards  for  Reference  Case  8. 


Temperature  T 

rci 

Thermal  Conductivity  k 
[W/(m*K)] 

0 

0.174 

250 

0.188 

1100 

0.188 

0.02  0.04  0.06  0.08 


Fig.  C-20  —  These  are  quarter  domains  showing  the  four  gnds  used  for  Reference  Case  8 
The  upper  left  is  Mesh  0,  the  upper  right  is  Mesh  1,  the  lower  left  is 
Mesh  2,  and  the  lower  right  is  Mesh  3. 


Figure  C-21  shows  the  temperature  contours  with  both  radiative  and  convective  heat  transfer 
within  the  cavities  after  two  hours  of  exposure  to  the  IS0834  fire.  Table  C-7  shows  the  results 
for  the  middle  flange  temperature.  HEATING  and  TASEF  are  in  good  agreement  when  only 
radiation  is  occurring  in  the  cavities.  The  addition  of  convection  to  the  cavities  has  a  more 
pronounced  impact  on  TASEF  than  HEATING.  For  both  cases,  the  agreement  between  the  two 
models  improves  as  t  increases  from  30  min.  Figure  C-22  shows  the  time  histories  of  the 
middle  flange  (  x  =  0  m,  y  =  -0.0925  m)  temperature.  HEATING  is  clearly  converging  to  its 


Fig.  C-2 1  —  Contours  of  the  temperature  (in  °C)  @  t  =  120  min  in  a  quarter  of  the 
domain  for  Reference  Case  8  for  Mesh  3, 


Table  C-7  -  Middle  flange  temperature  [°C]  for  the  two  models  under  different 
conditions  for  Mesh  3  (Reference  Case  8). 


Time 

(min) 

Radiation  Only 

Radiation  +  Convection 

TASEF 

HEATING 

%  Diff. 

TASEF 

HEATING 

%  Diff. 

mm 

226 

230 

1.77 

246 

232 

5.69 

518 

521 

0.579 

536 

523 

2.43 

90 

736 

738 

0.272 

746 

739 

120 

879 

880 

0.114 

885 

881 
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Fig.  C-22  —  Temperature  at  the  center  of  the  bottom  flange  for  Reference  Case  8. 
Heat  transfer  is  by  radiation  and  convection  within  the  cavities. 


own  solution.  This  tends  to  be  lower  than  TASEF’s  in  the  mid-time  range.  However,  the  two 
solutions  approach  each  other  as  the  temperature  begins  to  level  out.  The  temperature  after  one 
hour  of  exposure  to  the  IS0834  fire  along  the  L-shaped  path  traced  from  the  center  of  the  section 
(0  m,  0  m)  to  the  center  of  the  bottom  flange  (0  m,  -0.0925  m)  to  the  bottom  right-edge  of  the 
flange  (0.1  m,  -0.0925  m)  is  shown  in  Figure  C-23.  Once  again  it  is  evident  that  HEATING  is 
converging  to  a  solution  that  is  slightly  lower  than  TASEF’s  during  the  temperature-rise  time 
range. 


3.0  CONCLUSIONS 


HEATING  has  compared  successfully  with  Reference  Cases  1  through  8  of  the  benchmark 
proposed  by  [Wickstrom  et  al,  1999].  It  was  found  that  for  Cases  1  through  6,  64  cells  are 
sufficient  for  satisfactory  convergence  with  either  the  analytical  solution  or  with  the  results 
calculated  by  TASEF.  For  Reference  Case  7,  two  bottom  cells  suffice  to  yield  excellent  results 
for  the  one-dimensional  transfer  Case  (7a).  For  the  two  dimensional  case  (7b),  8  bottom  cells 
result  in  differences  smaller  than  3.2%.  For  the  more  complex  geometry  of  Reference  Case  8, 
the  results  for  the  most  refined  mesh  are  in  excellent  agreement  with  TASEF  when  only  radiative 
transfer  is  occurring  in  the  cavities.  When  convective  heat  transfer  was  added  to  the  cavities, 
HEATING  yielded  moderately  lower  temperatures  than  TASEF. 


C-27 
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Distance  s  (m) 

Fig.  C-23  —  Temperature  @  t  =  60  min  along  the  L-shaped  path  traced  from  the  center  to  the 
bottom  right-edge  of  the  flange  for  Reference  Case  8.  Heat  transfer  is  by 
radiation  and  convection  within  the  cavities. 
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6.0  NOMENCLATURE 


Roman 

A  Area  [m  ] 

c  Specific  heat  [J/(kg-K)] 

Cp  Specific  heat  at  constant  pressure  [J/(kg-K)] 

E  Emissivity-shape  factor  matrix  [] 

F  Shape  factor  matrix  [] 

/  A  generic  function  [varies] 

G  Absorption  factor  matrix  [] 

H  Effective  conductance  matrix  [m^j 

h  Heat  transfer  coefficient  or  exponent  [W/(m^-K)] 

k  Thermal  conductivity  [W/(m-K)] 

L  Distance  [m] 

M  A  number  of  nodes  [] 
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N  The  number  of  nodes  in  a  cavity  [] 

n  Normal  coordinate  [m] 

q  Heat  (energy)  [J] 

R  Reflectivity-shape  factor  matrix  [] 

S  Distance  [m] 

s  Distance  [m] 

T  Temperature  [°C] 

/  Time  [s] 

V  Volume  [m^] 

w  Width  [m] 

X  The  Cartesian  position  vector,  equal  to  (x,y  ,z)  [(m,m,m)] 

X  The  first  Cartesian  coordinate  [m] 

y  The  second  Cartesian  coordinate  [m] 

z  The  third  Cartesian  coordinate  [m] 

Greek 


a  Thermal  diffusivity  [m  /s] 

y  Ratio  of  the  specific  heats  [] 

A  Difference  operator  [] 

<1>  Transfer  factor  matrix  [] 

p  Density  [kg/m^] 

p  Reflectivity  [] 

cr  Stefan-Boltzmann  constant  [5.669x  1  O'*  W/(m^-K'*)] 

Superscripts 

•  Quantity  per  unit  time  [s'*] 

'  First  derivative 

'  Quantity  per  unit  length  [m'*] 

"  Quantity  per  unit  area  [m'  ] 

Quantity  per  unit  volume  [m'^J 

Subscripts 

0  Denotes  an  initial  quantity 
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Denotes  a  natural  convection  exponent 
An  index 
An  index 
An  index 
An  index 

An  index ;  denotes  a  natural  convection  quantity 

Denotes  a  forced  convection  quantity 

Denotes  a  gas  quantity 

Denotes  a  radiation  quantity 

Denotes  a  solid  quantity 

Denotes  a  surface  quantity 

Denotes  an  x  -coordinate  quantity 

Denotes  a  y  -coordinate  quantity 

Denotes  a  z  -coordinate  quantity 

Denotes  an  ambient  quantity 


APPENDIX  A:  DEVELOPMENT  OF  CAVITY  HEAT  TRANSFER  ALGORITHMS 

INTRODUCTION 

This  appendix  details  the  methods  employed  to  facilitate  the  calculation  of  heat  transfer  in 
voids  using  HEATING  7.3.  The  term  “connectors”  is  used  to  describe  data  blocks  for  the 
radiative  and  convective  conductance  between  any  two  nodes.  The  results  have  been  verified 
using  the  two  relevant  examples  recommended  in  [Wickstrom  &  Palsson,  1999;  Wickstrom, 
1999;  Palsson  8c  Wickstrdm,  2000]. 


DIFFUSE,  GRAY  SURFACE  RADIATIVE  HEAT  TRANSFER 

It  is  assumed  that  the  surfaces  participating  in  radiative  exchange  have  emissivities  that 
essentially  do  not  vary  with  the  wavelength  of  the  radiation  (i.e.,  they  are  gray)  and  that 
reflections  off  of  a  surface  are  diffuse.  The  ratio  of  the  mean  surface  roughness  to  the  median  (or 
bounding)  radiation  wavelength  serves  as  a  measure  for  the  type  of  reflection  [Siegel  &  Howell, 
1992].  If  it  is  less  than  one,  the  reflections  are  specular.  If  it  is  greater  than  one,  it  is  diffuse.  A 
FORTRAN  95  program,  h7_3radconn,  has  been  written  to  create  the  CONNECTORS  file  for 
two  dimensional,  diffuse-gray  radiative  transfer.  As  its  name  implies,  h7_3radconn  is  intended 
to  work  in  conjunction  with  version  7.3  of  HEATING.  It  performs  two  basic  tasks;  it  calculates 
the  node-to-node  shape  and  transfer  factors. 


SHAPE  FACTORS 

The  method  of  crossed  strings  [Siegel  &  Howell,  1992]  is  used  to  determine  the  shape  factors 
between  cells.  The  area  for  the  i  *  cell  is  denoted  .  The  shape  factor  from  cell  /  to  cell  j  is 

Fy .  The  method  of  crossed  strings  makes  use  of  the  two  shortest  distances  from  cell  edges,  Syy , 
(these  are  the  “straight  strings”),  and  the  two  longest  distances  that  can  be  measured  from 
one  cell’s  edge  to  another,  Ly^ ,  Ly2  (these  are  the  “crossed  strings”).  The  formula  for  the  shape 
factor  is 


(5) 

HEATING  7.3  actually  needs  values  at  nodes  so  the  area  for  each  node  is  the  sum  of  half  the 
area  of  the  two  cells  surrounding  the  node.  (The  emissivity  for  the  node  adjoining  two  disparate 
surfaces  is  determined  in  a  similar  fashion.  The  half-cell-area  average  of  the  emissivities  of  the 
two  abutting  materials  is  used.)  However,  for  internal  comers,  the  projected  area  (i.e.,  the 
hypotenuse)  is  used. 
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ABSORPTION  FACTORS 


HEATING’S  method  for  calculating  the  radiative  heat  transfer  for  node  i ,  ,  is 

(6) 

y=i 

0y  is  the  (Hottel)  transfer  factor.  It  is  the  fraction  of  the  black  body  thermal  energy  between 
cells  i  and  j  that  can  be  exchanged  due  to  the  geometry  and  reflectivity  of  ALL  the 
participating  surfaces  (nodes  in  this  case).  0y  can  be  related  to  the  absorption  factor,  Gy ,  via 
the  formula 


SiGy. 


(7) 


Si  is  the  emissivity  for  surface  i .  The  reflectivity  for  surface  k,  Pi,  is 

Pi={l-Si).  (8) 


Gy  is  the  fraction  of  the  emission  from  surface  /  that  reaches  surface  j  and  is  absorbed.  Like 
the  shape  factor,  it  sums  to  unity  for  each  surface.  I.e., 

Vi.  (9) 

;=i 

The  reciprocity  relation  is 

MGy^AySjGj,  (10) 


Note  that  this  implies  that 


A0  =A0. 


(11) 


which  is  completely  analogous  with  the  shape  factor  reciprocity  relation.  Although  for  a  flat 
surface  F),  =  0,  this  is  generally  not  so  for  G,.,.  since  energy  emitted  by  surface  i  can  be  reflected 

back  to  surface  i  and  be  subsequently  (fractionally)  absorbed.  An  analysis  of  the  energy  that  is 
absorbed  by  each  surface  [Gebhart,  1971;  Siegel  &  Howell,  1992]  results  in  the  matrix  equation 

RG  =  E.  (12) 

The  components  of  the  reflectivity-shape  factor  matrix,  R ,  are 
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-PA 

i>j 

1 

i  =  j 

(13) 

-Pipji 

j>i 

and  for  the  emissivity-shape  factor  matrix,  E ,  they  are 

E,=s,F,.  (14) 

Experience  has  shown  that  if  all  the  <1,  then  the  matrix  R  is  generally  dense,  strictly 

row-and-column  diagonally  dominant,  and  non-symmetric.  If  all  the  =1,  then  the  matrix  R 
is  generally  sparse,  strictly  row-and-column  diagonally  dominant,  and  symmetric. 

The  program  that  determines  the  shape  factor  also  solves  Equation  (Error!  Reference 
source  not  found.).  The  user  can  choose  different  numeric  methods  to  achieve  the  solution. 
The  available  methods  are  LU  decomposition  [Anderson  et  al.,  1999],  Gauss-Jordan  elimination 
[Press  et  al.,  1996;  Press  et  al.,  1993],  diagonal  pivoting  [Anderson  et  al.,  1999],  QR 
decomposition  [Anderson  et  al.,  1999],  and  singular  value  decomposition  [Anderson  et  al., 
1999].  It  is  recommended  that  for  each  new  geometry,  the  HEATING  7.3  connections  file  be 
determined  by  at  least  two  of  these  methods. 


CONVECTIVE  HEAT  TRANSFER  IN  A  CAVITY 

HEATING  7.3  can  readily  model  convective  heat  transfer  at  external  solid  boundaries 
because  the  ambient  temperature  is  assumed  to  be  known.  If  the  boundary  is  part  of  a  closed 
cavity,  the  background  temperature  will  be  changing.  A  method  is  presented  in  this  appendix 
that  addresses  this  issue  of  determining  T„{t) ,  the  background  temperature  as  a  function  of  time. 
As  before,  HEATING  7.3  needs  values  at  nodes  so  the  area  for  each  node  is  the  sum  of  half  the 
area  of  the  two  cells  surrounding  the  node.  For  any  comer,  the  half  area  of  each  branch  is  used 
since  convection  depends  on  the  “wetted”  area. 

[Wickstrom  &  Palsson,  1999]  suggest  letting  T„{t)  be  equal  to  the  average  temperature  for 
all  the  cells  in  the  cavity.  In  order  to  impliment  this  in  HEATING  using  connections,  first 
consider  how  HEATING  models  convective  heat  transfer  for  node  i : 

g,=A,h,{T,-T,(0).  (15) 

For  now  it  is  assumed  that  the  heat  transfer  coefficient,  h. ,  is  constant.  HEATING  also 
allows  node  i  to  participate  with  any  number  (say,  M )  of  nodes, 

-T,l  (16) 

M 
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where  Hy  is  the  effective  conductance  matrix  associated  with  nodes  i  and  j .  The  effective 
conductance  is 


h,Hy.  (17) 

Applying  Wickstrom  and  PMsson’s  suggestion  to  Equation  (Error!  Reference  source  not 
found.)  results  in 


qi  =  A,h. 


^  r(5,t)ds 


(18) 


where  A  is  the  total  length  (per  unit  depth)  of  the  cavity  perimeter  and  s  is  the  coordinate  about 
the  closed  loop  of  the  cavity.  Replacing  the  integral  with  the  middle  Riemann  sum  (i.e.,  by  the 
mid-point  method  [Burden  &  Faires,  1989])  yields 


^  N 


q,  =  A,h, 


T,- 


y=i 


(19) 


For  the  present  two-dimensional  Cartesian  coordinate  system,  Aj  is  the  cell  length  associated 
with  node  j .  Equation  (Error!  Reference  source  not  found.)  can  be  written  as 


q,  =  4/i, 


f  N  N 

y=i  y=i 


(20) 


because  ^Aj 

M 


=  A .  Further  rearrangement  results  in 


^  i  \ 


(21) 


Here,  A.Aj  /  A  is,  by  comparison  with  Equation  (Error!  Reference  source  not  found.),  the 
conductance  matrix.  This  conductance  matrix. 


AA, 


(22) 
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is  calculated  as  the  input  for  HEATING  7.3  CONNECTOR  lines.  Note  that  the  procedure 
presented  in  this  section  could  also  have  been  carried  out  by  using  the  arithmetic  mean.  The 
result  would  replace  the  term  Aj  I  A  m  Equation  (Error!  Reference  source  not  found.)  with 

M  N .  In  fact,  if  the  step-size  around  the  cavity  is  constant,  Aj!  A  =  \l  N .  However,  for 
nonuniform  gridding,  the  area-averaged  approach  is  clearly  the  correct  one. 
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HEATING  7.3  VALIDATION  FOR  FIRE  RESISTIVE  APPLICATIONS 


1.0  INTRODUCTION 

A  series  of  eight  validation  cases  were  selected  to  validate  the  thermal  computer  model 
HEATING.  These  entail  thermal  calculations  of  increasing  complexity.  The  algorithms  used  by 
HEATING  have  been  previously  verified  in  Trelles  et  al.  [2003]  and  Bryan  et  al.  [1986].  The 
focus  of  these  validation  cases  is  to  compare  the  calculated  results  of  HEATING  with  test  data 
for  various  configurations. 

The  validation  cases  were  selected  to  address  the  following  aspects  of  a  numerical  heat 
transfer  calculation: 

•  Variable  convective  and  thermal  radiation  boundary  conditions; 

•  Variable  thermal  material  properties; 

•  Multiple  materials;  and 

•  Radiation  and  convection  in  voids. 


In  addition,  the  sensitivity  of  the  results  to  the  boundary  condition  parameters,  material 
properties,  convergence  criteria,  and  the  numerical  solution  method  were  investigated.  This  was 
accomplished  by  varying  the  boundary  condition  parameters  and  material  properties  for  which 
there  are  degrees  of  uncertainty,  increasing  the  resolution  of  the  solution  grid,  using  different 
numerical  methods,  and  changing  the  solution  error  tolerances. 

The  validation  cases  were  selected  fi'om  available  test  data  involving  fire  tests  of  building 
materials,  insulation  materials,  and  composite  materials.  The  availability  of  material  properties, 
the  type  of  exposure  conditions,  and  the  reliability  of  the  exposure  and  sample  temperature  data 
were  also  factors  in  selecting  the  validation  cases. 


2.0  BACKGROUND 


HEATING  was  developed  at  the  Oak  Ridge  National  Laboratory's  (ORNL)  Radiation  Safety 
Information  Computational  Center  (RSICC)  [Childs,  1998].  The  model  was  created  to  analyze 
the  thermal  impact  of  the  U.S.  Energy  Department's  various  high  energy  projects,  including  all 
aspects  of  radiation  shielding.  It  has  one  of  the  longest  development  histories  of  any 
computational  heat  transfer  software  package  [Fowler  et  al.,  1959;  Liguori  et  al.,  1961;  Turner  et 
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al.,  1969;  Turner  et  al.,  1971;  Turner  et  al.,  1977;  Becker,  1977;  Elrod  et  al.,  1984;  Bryan  et  al, 
1986;  Childs,  1990;  Childs  et  al.,  1990;  and  Childs,  1991].  Verification  studies  and  tools  for 
HEATING  can  be  found  in  Bryan  et  al.  [1986],  Chu,  1989],  and  more  recently  in  Trelles  et  al. 
[2003]. HEATING  is  a  field  model.  It  solves,  using  finite  volume  techniques  [Patankar,  1980; 
Childs,  1998],  the  following  heat  equation  [Arpaci,  1966;  Carslaw  etal.,  1959]: 

yOC^  =  V-(/:Vr)+^",  23) 

subject  to  the  general  boundary  condition: 

?:  =  =  {*/ +*,(7;"  +ri)(r, +2-j)(2;  -r.) 

on 

or  (24) 

7:=/.(V) 

where  T  is  the  material  internal  temperature  (K),  p  is  the  density  (kg/m^),  c  is  the  specific  heat 
(kJ/kg),  t  is  the  time  (s),  q"  is  the  volumetric  heat  generation  rate  (kJ/kg-s),  q]  is  the  heat  flux 
for  surface  5'  (kW/m^),  n  is  the  coordinate  normal  to  the  surface,  k  is  the  thermal  conductivity 
(kW/m-K),  is  the  material  surface  temperature  (K),  X  a  temperature  function  associated 

with  surface  ^ is  the  ambient  temperature  (K),  hj.  is  the  heat  transfer  coefficient  for  forced 
convection  (kW/m^-K),  is  the  heat  transfer  coefficient  for  natural  convection  (kW/m^-K*'), 
is  the  exponent  for  natural  convection,  is  the  thermal  radiation  heat  transfer  coefficient 
(kW/m^-K^). 

HEATING  7.2i  and  7.3  are  the  most  recent  developments  in  a  series  of  heat-transfer  codes 
distributed  by  RSICC  as  SCA-1 /HEATINGS  and  PSRl 99/HEATING  6.  HEATING7  can  solve 
steady-state  and/or  transient  heat  conduction  problems  in  one-,  two-,  and  three-dimensional 
Cartesian,  cylindrical,  or  spherical  coordinates.  The  following  are  some  notable  features  of 
HEATING  7. 

•  Multiple  materials  may  be  used. 

•  Material  properties  may  be  both  time-  and/or  temperature-dependent.  Thermal 
conductivity  may  be  anisotropic  and  materials  may  undergo  a  phase  change. 

•  Material  properties  may  be  inputed  directly  or  extracted  fi'om  a  material  properties  library 
provided  with  the  model. 

•  Heat  generation  rates  may  be  a  function  of  time,  temperature,  and/or  position. 

•  Boundary  temperatures  may  be  time-  and/or  position-dependent. 

•  Boundary  conditions  may  be  surface-to-environment  or  surface-to-surface.  They  may 
consist  of  a  specified  temperature  or  a  combination  of  prescribed  heat  flux,  forced 
convection,  natural  convection,  and  thermal  radiation. 

•  Boundary  condition  parameters  may  be  time-  and/or  temperature  dependent. 
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•  General  gray-body  radiation  problems  may  be  modeled  with  user-defined  factors  for 
radiant  exchange.  Surfaces  that  are  not  in  contact  can  exchange  heat  via  custom 
connectors.  The  mesh  spacing  may  be  variable  along  each  axis. 

The  HEATING  calculations  were  performed  on  the  following  computer  systems; 

Linux  workstation  with  two  Pentium  HI  Xeon  processors,  1  GB  of  RAM  memory,  and 
19  GB  of  storage; 

Silicon  Graphics,  Inc.  Octane2  graphics  workstation  with  V6®  graphics,  two  R12,000A 
processors,  2  GB  of  RAM  memory,  and  more  than  136  GB  of  storage;  and; 

Dell  workstation  with  a  single  Pentium  4  processor,  1  GB  of  RAM  memory,  and  80  GB 
of  storage. 

3.0  VALroATION 


The  validation  set  for  HEATING  consists  of  eight  cases  of  increasing  complexity.  The  basis 
for  comparison  between  a  validation  case  and  HEATING  includes  temperatures  at  a  specific 
time  and  location,  transient  temperatures  at  a  specific  position,  and  average  temperature  at  a 
specific  time.  The  particular  metrics  vary  among  the  validation  cases  due  to  the  nature  of  the 
data  available  for  comparison.  The  validation  cases  and  the  comparison  basis  are  summarized  in 
Table  D-1  below. 


Table  D-1  —  Summary  of  HEATING  Validation  Cases 


Validation 

Case 

Description 

Key  Elements 

Validation  Basis 

1 

Bare  steel  plate  heated  by 
known  transient  temperature 
profile. 

1-D  thermal  analysis,  material 
properties  of  steel  well 
documented,  boundary  conditions 
deduced  fi'om  measurements. 

Transient  surface 
temperatures. 

2 

Concrete  section  heated  by 
ASTM  El  19  Standard  Time- 
Temperature  profile. 

1-D  thermal  analysis,  properties 
for  concrete  documented,  includes 
water  vaporization  in  concrete, 
thermocouple  pads. 

Fire  resistance  rating, 
unexposed  surface 
temperature  at  a  fixed 
time. 

3 

Aluminum  protected  with 
Unifrax  high  temperature 
mineral  fiber  blanket; 
assembly  exposed  to  test 
conditions  specified  by  IMO 
A.754(18). 

1-D  thermal  analysis  with  multiple 
materials.  Thermal  material 
properties  of  insulation  material 
are  proprietary,  thermocouple 
pads. 

Fire  resistance  rating, 
transient  surface 
temperatures. 

4 

Steel  protected  with  mineral 
fiber  insulation  -  exposed  to 
an  empirical  furnace  curve 

1-D  thermal  analysis,  material 
properties  of  steel  and  Structo- 
gard  are  well  documented. 

Transient  surface 
temperatures. 
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Table  D-1  —  Summary  of  HEATING  Validation  Cases  (Continued) 


Validation 

Case 

Description 

Key  Elements 

Validation  Basis 

5 

Unprotected  steel  column 
filled  with  concrete  in  web 
area  exposed  to  ASTM  El  19 
Standard  Time-Temperature 
profile. 

2-D  thermal  analysis  analogous  to 
Validation 

Case  2. 

Transient  average 
temperatures,  fire 
resistance  rating. 

6 

Steel  beam  protected  with 
Blazeshield  D-C/F  insulation 
exposed  to  ASTM  El  19 
Standard  Time-Temperature 
profile  (UL  Design  1907). 

2-D  thermal  analysis  involving 
steel  W8X28  beam  protected 
with  Blazeshield  D-C/F  spray 
applied  mineral  fiber  insulation 
supporting  improtected  concrete 
floor,  material  properties 
documented. 

Fire  resistance  rating, 
average  steel 
temperature  at  a  fixed 
time. 

7 

Steel  beam  protected  with 
Blazeshield  D-C/F  insulation 
exposed  to  ASTM  El  19 
Standard  Time-Temperature 
profile  (UL  Design  P809). 

2-D  thermal  analysis  involving 
steel  W8X28  beam  and  concrete 
floor  protected  with  Blazeshield 
D-C/F  spray  applied  mineral 
fiber  insulation,  material 
properties  documented. 

Fire  resistance  rating, 
average  steel 
temperature  at  a  fixed 
time. 

8 

Gypsum-steel  wall  assembly 
exposed  to  CAN/ULC-SIOI- 
M89  Standard  Time- 
Temperature  profile. 

2-D  thermal  analysis  involving  a 
radiation  and  convection  heat 
transfer  in  a  cavity  space, 
material  properties  of  gypsum 
and  steel  well  documented. 

Transient  surface 
temperatures. 

Convergence  was  verified  for  each  validation  case  by  increasing  the  number  of  nodes, 
decreasing  the  time  step,  decreasing  the  error  tolerance,  and/or  altering  the  solution  method. 


3.1  Validation  Case  1:  Bare  Steel  Deck  Exposed  to  Known  Transient  Temperature 


Validation  Case  1  consists  of  bare  steel  deck  exposed  to  a  post-flashover  fire  environment. 
The  post-flashover  fire  was  approximated  in  a  cubical  enclosure  measuring  2.4  m  by  2.4  m  by 
2.4  m  using  a  JP-5  spray  fire  [Leonard  et  ai,  1991].  The  enclosure  was  constructed  of  0.95  cm 
thick  steel  plates.  The  temperature  in  the  compartment  as  well  as  the  exposed  and  unexposed 
surface  temperatures  of  the  steel  were  measured  as  a  function  of  time. 

The  validation  basis  is  the  average  deck  exposed  and  unexposed  temperature  as  a  function  of 
time.  The  input  parameters  include  the  thermal  material  properties  of  the  steel  and  the  exposed 
and  unexposed  surface  boundary  conditions.  A  one-dimensional  thermal  analysis  is  used  to 
calculate  the  exposed  and  unexposed  surface  temperature  of  the  steel  plate  as  a  function  of  time. 
Figure  D-1  depicts  the  configuration  modeled  by  HEATING.  The  thermal  material  properties 
for  steel  were  obtained  from  Abrams  [1979].  Figure  D-2  shows  the  thermal  conductivity  and 
heat  capacity  as  a  function  of  temperature  for  A36  steel.  The  spike  in  the  heat  capacity 
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corresponds  to  the  perlite-austenite  transformation,  which  typically  occurs  between  700°C  and 
850°C  [Abrams,  1979].  A  density  of  7,800  kg/m^  is  assumed  for  the  steel  [Holman,  1990]. 


p —  0.095  m - 

Fig.  D-1  —  Validation  Case  1  (Not  to  Scale) 
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Temperature  (X) 

Fig.  D-2  — Thermal  Material  Properties  for  A36  Steel  [Abrams,  1979] 

Boundary  Conditions 

A  natural  convection  boundary  condition  was  applied  to  the  exposed  and  unexposed  surface 
of  the  steel  in  the  following  form  [Holman,  1990]: 

«  =  *.-(|7;-7’.|r  (3) 

where  q1  is  the  convective  heat  flux  boimdary  condition  (W/m^),  hn  is  the  effective  natural 
convection  coefficient  (W/m^-°C'‘'^),  Ts  is  the  steel  surface  temperature  (K),  and  Ta  is  the 
ambient/exposure  temperature  (K).  Equation  3  is  applicable  to  heated  plates  facing  up  (on  the 
unexposed  boundary)  or  cooled  plates  facing  down  (on  the  exposed  boundary  condition)  when 
the  convection  coefficient  is  1.52  W/m^-^C"*^  [Holman,  1990]. 

A  thermal  radiation  boundary  condition  was  applied  to  the  exposed  and  unexposed  surface  of 
the  steel  in  the  following  form  assuming  a  gray  gas  and/or  surfaces  that  have  the  same  radiosity 
[Holman,  1990;  Bird  et  al.  I960]: 

j;=£,<TF(r,*-r;)  (4) 

where  ql  is  the  thermal  radiation  heat  flux  boundary  condition  (W/m\  is  the  resultant 
(effective)  emissivity  between  the  steel  surface  and  the  ambient  environment,  and  F  is  the  view 
factor  between  the  exposure  and  the  surface.  Note  that  the  heat  flux  would  be  negative  {i.e.  the 
steel  gains  energy)  if  the  ambient  temperature  is  greater  than  the  steel  temperature  and  vice 
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versa.  The  resultant  effective  emissivity  was  determined  using  the  following  equation  [Barnett, 
1989;  Bird  et  al.  1960;  Pettersson  et  al.,  1976]: 


where  is  the  emissivity  of  the  ambient  environment  and  is  the  emissivity  of  the  steel.  The 

emissivity  of  steel  is  0.8  and  the  emissivity  of  the  ambient  environment  on  the  unexposed  side  of 
the  steel  is  assumed  to  be  0.95  [Holman,  1990;  Sparrow  et  al,  1978;  Incropera  et  al.,  1985].  The 
emissivity  of  the  fire  exposure/compartment  surfaces  was  determined  fi-om  the  test  data  to  be 
approximately  0.7,  which  takes  into  account  surfaces  which  may  be  at  different  temperatures  due 
to  the  formation  of  a  smoke  layer  in  the  test  compartment.  The  view  factor  was  unity  between 
ambient  environment  and  the  steel  on  the  exposed  and  unexposed  surfaces.  As  such,  the 
effective  emissivity  on  the  exposed  steel  surface  was  0.6  and  0.77  on  the  unexposed  steel 
surface. 
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Results 


The  implicit  Crank-Nicolson  method  [Childs,  1998;  Smith,  1987;  Jaluria  et  ah,  1986]  was 
used  to  solve  for  the  steel  plate  surface  temperature.  The  tolerance  for  convergence  was  1X10'® 
and  the  solution  time  step  was  0.5  seconds.  Decreasing  the  tolerance,  decreasing  the  time  step, 
or  using  the  Classical  Explicit  Procedure  or  Classical  Implicit  Procedure  produced  the  same 
results. 

The  validation  results  are  shown  in  Figure  D-3  for  a  model  using  twenty  cells.  Also  shown 
in  Figure  D-3  is  the  measured  exposure  temperature  profile.  The  results  did  not  change  when  the 
number  of  cells  was  increased  to  one  hundred. 


Fig.  D-3  —  Validation  Case  1  Temperature  Results  -  Twenty  Cell  Model 


Figure  D-3  shows  good  agreement  between  the  measured  and  calculated  exposed  steel 
surface  temperature.  However,  the  measured  unexposed  surface  temperature  is  considerably 
lower  than  the  calculated  unexposed  surface  temperature.  The  test  data  indicate  that  there  is  an 
approximately  90  °C  temperature  difference  between  the  exposed  and  unexposed  steel  surfaces 
whereas  the  calculated  temperature  difference  is  approximately  20  °C.  The  measured 
temperature  difference  likely  includes  heat  loss  effects  around  the  thermocouples  that  are  not 
included  in  the  thermal  model.  This  may  be  seen  when  considering  the  quantity  of  energy  that 
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the  steel  transmits  via  conduction  from  one  surface  to  the  other.  The  following  equation  from 
Holman  [1990]  can  be  used  to  calculate  this  effect: 

(6) 

where  is  the  energy  conducted  through  the  steel  (kW/m^),  k  is  the  thermal  conductivity  of  the 

steel  (kW/m-°C),  S  is  the  thickness  of  the  steel  (m),  and  Tj  and  T2  are  the  exposed  and 
unexposed  surface  temperatures  (°C),  respectively.  The  thermal  conductivity  of  steel  is  about 
0.03  kW/m-°C  when  the  steel  temperature  is  800  °C,  as  seen  in  Figure  D-2.  If  the  temperature 
difference  between  the  exposed  and  unexposed  surfaces  were  90  °C  as  measured  25  minutes 
after  initial  exposure,  then  the  energy  conducted  through  the  steel  would  be  about  284  kW/m^,  as 
calculated  using  Equation  6  for  9.5  mm  thick  steel.  Because  the  steel  temperature  is  increasing 
slowly  at  25  minutes,  most  of  this  energy  would  have  to  be  radiated  and  convected  into  the 
ambient  area.  This  would  require  a  blackbody  surface  temperature  around  1,200°C,  which 
greatly  exceeds  the  measured  surface  temperatures,  exposed  or  xmexposed.  If  the  temperature 
difference  between  the  exposed  and  unexposed  surfaces  is  closer  to  20°C,  as  calculated  by 
HEATING,  then  the  energy  conducted  through  the  steel  would  be  about  60  kW/m^,  which  is 
closer  to  the  expected  boundary  heat  loss  for  a  blackbody  surface  temperature  in  the  800°C  range 
and  an  ambient  temperature  of  20°C. 

Sensitivity  of  Results  to  Input  Parameter  Assumptions 

The  sensitivity  of  the  calculated  temperatures  to  the  assumed  boundary  conditions  was 
investigated  by  increasing  and  decreasing  the  effective  emissivity  of  the  exposed  and  unexposed 
surfaces.  The  convection  coefficient  was  also  increased  and  decreased  by  twenty-five  percent. 
The  results  are  shown  in  Figure  D-4  for  the  exposed  steel  surface  temperature.  The  figure 
indicates  that  the  calculated  steel  surface  temperature  may  vary  by  plus  or  minus  50°C  depending 
on  the  selection  of  the  effective  emissivity.  The  calculated  temperature  is  less  sensitive  to  the 
convection  coefficient.  The  figure  also  indicates  that  increasing  the  effective  emissivity  has 
opposite  impacts  on  the  calculated  temperature  when  applied  to  the  exposed  versus  unexposed 
surface.  The  same  is  true  for  decreasing  the  effective  emissivity.  This  suggests  that  there  is 
some  compensating  effects  such  that  the  combined  impact  of  an  incorrect  effective  emissivity 
would  be  less  than  that  shown  in  Figure  D-4,  as  it  is  unlikely  that  the  effective  emissivity  would 
be  increased  on  one  surface  and  decreased  on  another. 
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3.2  Validation  Case  2:  Concrete  Wall  Section  Exposed  to  the  ASTM  E119  Standard 
Time-Temperature  Fire  Exposure 


Validation  Case  2  involves  a  concrete  wall  section  exposed  to  the  ASTM  El  19  Standard 
Time-Temperature  Curve.  The  ASTM  El  19  Standard  Time-Temperature  fire  exposure  is 
intended  to  bound  a  compartment  fire  temperature  exposure  during  post-flashover  fire  conditions 
[ASTM  El  19-98,  1999;  Lie,  1995].  The  wall  is  assumed  to  consist  of  one  of  three  types  of 
concrete:  carbonate  aggregate  normal  weight  concrete,  siliceous  aggregate  normal  weight 
concrete,  and  lightweight  aggregate  concrete.  Carbonate  aggregate  normal  weight  concrete  is 
concrete  made  with  aggregate  consisting  mainly  of  calcium  carbonate  or  magnesium  carbonate. 


. . . . . .  Predicted  temperature  (baseline) 

•  Measured  temperature 

-  Exposure  surface  effective  emissivity  increased  to  0.75 

—  —  —  Exposure  surface  effective  emissivity  decreased  to  0.5 

- Unexposed  surface  effective  emissivity  increased  to  0.9 

- —  Unexposed  surface  effective  emissivity  decreased  to  0.65 

.  Convection  increased  by  25% 

-  Convection  decreased  by  25% 


Time  (minutes) 

Fig.  D-4  —  Sensitivity  of  the  Calculated  Exposed  Steel  Surface  Temperature  Predictions 
to  Uncertainty  in  the  Boundary  Condition  Parameters 


typically  limestone  or  dolomite,  and  has  a  dry  unit  weight  between  2,200  kg/m3  and  2,400  kg/m3 
[ASCE  29-99,  1998;  Harmathy,  1983;  Lie,  1972].  Siliceous  ag^egate  normal  weight  concrete 
also  has  a  dry  unit  weight  between  2,200  kg/m^  and  2,400  kg/rn^  and  is  made  with  aggregate 
consisting  mainly  of  silica  or  compounds  other  than  calcium  carbonate  or  magnesium  carbonate 
[ASCE  29-99, 1998;  Harmathy,  1983;  Lie,  1972].  Lightweight  aggregate  concrete  is  made  with 
aggregates  of  expanded  clay,  shale,  slag,  slate  or  sintered  fly  ash  and  has  a  dry  unit  weight  of 
1,360  -  1,840  kg/m^  [ASCE  29-99,  1998;  Lie,  1972].  Because  the  thermal  material  properties 
vary  with  the  nature  of  the  aggregate,  the  fire  resistance  rating  of  a  concrete  wall  varies  with  the 
type  of  aggregate.  As  such,  this  section  considers  wall  sections  made  with  each  of  the  three 
aggregates  separately. 
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The  basis  for  the  validation  case  is  the  fire  resistance  rating  of  the  concrete  wall,  which  is 
defined  as  a  temperature  rise  of  139°C  on  the  unexposed  side  of  the  wall  per  ASTM  El  19 
[1999].  A  one-dimensional  thermal  analysis  is  used  to  calculate  the  fire  resistance  rating  of  the 
concrete  wall  sections.  Figure  D-5  depicts  the  configuration  modeled  by  HEATING.  Fire 
resistance  ratings  between  1  hr  and  4  hrs  for  each  of  the  three  types  of  concrete  noted  above  are 
used  for  the  validation  basis  [ASCE  29-99,  1998;  Fleischmann,  1995].  The  fire  resistance 
ratings  for  the  wall  sections  are  based  on  tests  reported  by  Abrams  et  al.  [1968].  Per  ASTM 


El  19,  thermocouples  used  to  measure  the  unexposed  wall  surface  temperature  are  covered 
with  pads  in  an  attempt  to  simulate  contact  with  a  combustible  material.  Because  the 
thermocouple  pads  affect  the  temperature  profile  in  the  wall,  they  are  included  in  the  thermal 
model  used  in  this  validation  case. 


ASTM  E119  RATING; 
T3(t=rating)=  139°C  +  T. 


-T(t) 


AMBIENT  CONDITIONS,  \ 

T.  J 


-T,(t) 

->  CONVECTION  =  hj(Tj-T.)‘ 


->  THERMAL  RADIATION  =  eE,6(T3^-T.‘) 


K 


CONCRETE 

(THICKNESS 

VARIES) 


9.5mm  THICK 

THERMOCOUPLE 

PAD 


Fig.  D-5  —  Validation  Case  2  (Not  to  Scale) 
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Material  Properties 

The  material  properties  that  are  necessary  for  calculating  the  fire  resistance  rating  of  the 
concrete  wall  section  are  the  thermal  conductivity,  the  heat  capacity,  the  density,  and  the  percent 
moisture  content  of  the  concrete  and  the  thermal  conductivity,  heat  capacity,  and  the  density  of 
the  thermocouple  pad.  Figures  D-6  and  D-7  show  the  thermal  conductivity  and  heat  capacity  as 
a  function  of  temperature  for  carbonate  and  siliceous  aggregate  normal  weight  concretes  and 
lightweight  concrete  assumed  in  this  evaluation  [Harmathy,  1983;  Lie,  1972].  Material 
properties  vary  considerably  amongst  concretes  [Flynn,  1999]  and  were  the  primary  focus  of  the 
sensitivity  evaluation.  The  actual  material  properties  of  the  concrete  that  was  the  basis  for  the 
fire  resistance  ratings  was  unknown.  A  general  assumption  applied  to  this  validation  case  was 
that  the  upper  bound  thermal  conductivity  values  for  the  particular  aggregate  apply.  The  density 
of  each  concrete  is  as  follows: 

•  Carbonate  aggregate,  normal  weight:  2,136  kg/m  ; 

•  Siliceous  aggregate,  normal  weight:  2,300  kg/m  ;  and 

•  Lightweight:  1,150  kg/m^. 


Carbonate  aggregate  normal  weight  concrete 
Siliceous  aggregate  normal  weight  concrete 
Lightweight  expanded  shale  concrete 


Fig.  D-6  —  Thermal  Conductivity  of  Various  Concrete  Aggregates 
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Carbonate  aggregate  normal  weight  concrete 
Siliceous  aggregate  normal  weight  concrete 
Lightweight  expanded  shale  concrete 


Fig.  D-7  —  Heat  Capacity  of  Various  Concrete  Aggregates 


The  moisture  content  of  the  concrete  may  vary  depending  on  the  local  environmental 
conditions,  the  cure  time,  and  other  factors.  The  moisture  content  may  range  from  zero  percent 
to  twelve  percent  or  greater  [Flyim,  1999].  A  1 14  percent  by  mass  moisture  content  is  the 
minimum  moisture  content  of  the  samples  tested  by  Harmathy  [1983]  and  is  assumed  for  this 
validation  case.  The  vaporization  of  the  moisture  is  included  in  the  model  by  applying  the  total 
vaporization  energy  to  the  heat  capacity  term  in  HEATING  over  a  10°C  range  starting  at  100°C. 
The  heat  of  vaporization  of  water  at  100°C  and  1  atm  pressure  is  2,258  kJ/kg  [Jones  et  al.,  1986] 
The  resulting  vaporization  energy  and  values  added  to  the  heat  capacity  terms  in  HEATING  are 
summarized  in  Table  D-2. 

Table  D-2  —  Moisture  Vaporization  Energy  for  Concrete 


Concrete 

Density  (kg/m3) 

Moisture  Mass 
Density  (kg/m^) 

Vaporization 
Energy  (kJ/kg) 

Heat  Capacity 
Increase 
(kJ/kg-°C) 

Carbonate 

2,136 

32 

72,260 

34.5 

17.3 
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The  thermal  conductivity  and  (indirectly)  the  density  of  the  thermocouple  pads  are  prescribed 
by  ASTM  El  19-98  [1999].  The  thermal  conductivity  of  refractory  fiber  pads  should  be  0.053 
plus  or  minus  0.004  W/m-K  and  the  density  should  be  305  kg/m^  but  may  fall  between  159 
kg/m^  and  519  kg/m^.  The  large  density  range  is  the  result  of  tolerances  in  the  pad  area, 
thickness,  and  dry  weight.  Due  to  the  uncertainty  in  the  density  of  the  pad  material,  an  average 
heat  capacity  for  mineral  fiber  over  a  600  °C  temperature  range  [Harmathy,  1983]  was  assumed 
(960  kJ/kg-K). 

Boundary  Conditions 

A  natural  convection  boundary  condition  was  applied  to  the  exposed  concrete  surface  and 
unexposed  thermocouple  pad  surface  as  described  by  Equation  3.  The  convection  coefficient  for 
vertical  planes  was  assumed  to  be  1.31  W/m^-°C*^^  [Holman,  1990].  A  thermal  radiation 
boundary  condition  was  applied  to  the  exposed  and  unexposed  surface  of  the  concrete  as 
described  by  Equation  4.  An  emissivity  of  0.91  was  assumed  for  the  concrete  surfaces  and  the 
thermocouple  pad  and  a  value  of  0.95  was  assumed  for  all  ambient  surroundings  (exposed  and 
unexposed)  [Holman,  1990;  Sparrow  et  al,  1978;  Incropera  et  al.,  1985].  The  resultant  or 
effective  emissivity  for  the  exposed  and  unexposed  surfaces,  as  calculated  using  Equation  5,  was 
0.87.  The  view  factor  was  unity  between  ambient  environment  and  the  concrete  and  the 
thermocouple  pad  surfaces. 

The  exposure  temperature  was  the  ASTM  El  19-98  Standard  Time-Temperature  profile 
[1999],  which  maybe  estimated  using  the  following  curve  fit  equations: 

^<120 

(7) 

120  <t  <240 

where  7/ is  the  exposure  temperature  (°C)  and  t  is  the  time  after  the  start  of  the  exposure 
(minutes).  The  ambient  exposure  temperature  is  23.9°C  such  that  the  critical  unexposed  surface 
temperature  is  162.8°C. 

Results 

The  implicit  Crank-Nicolson  method  [Childs,  1998;  Smith,  1987;  Jaluria  et  al.,  1986]  was 
used  to  solve  for  the  unexposed  concrete  surface  temperature.  The  tolerance  for  convergence 
was  1  X  10‘®  and  the  solution  time  step  was  0.5  seconds.  Decreasing  the  tolerance,  decreasing 
the  time  step,  or  using  the  Classical  Explicit  Procedure  or  Classical  Implicit  Procedure  produced 
the  same  results.  A  total  of  one  hundred  cells  was  used  in  the  concrete  and  twenty  in  the 
thermocouple  pad.  The  results  did  not  change  when  the  number  of  cells  was  increased  in  either 
material. 

The  thickness  of  the  concrete  varies  with  the  fire  resistance  rating  and  the  type  of  concrete 
aggregate  [ASCE  29-99,  1998;  Fleischmann,  1995].  Table  D-3  summarizes  the  wall  thickness 


r(^)=_^+^l^_(i7.8) 
3.14  +  t  490  +  t  '  ' 


r.  (t)  =  1010  +  ( ^  --]  •  (83) 

I  120  j  ^  ' 


and  ASTM  El  19  fire  resistance  rating  for  the  cases  considered.  The  thickness  of  the 
thermocouple  pads  is  9.5  mm  as  specified  by  ASTM  El  19-98  [1999].  The  calculated  fire 
resistance  rating,  which  is  the  time  at  which  the  unexposed  concrete  surface  reaches  163°C 
assuming  an  initial  temperature  of  24°C,  is  compared  to  the  listed  fire  resistance  rating  in 
Table  D-3. 

Table  D-3  —  Concrete  Wall  Fire  Resistance  -  120  Cell  Model 


Type  of  Concrete 

Thickness  (mm) 

Listed  Fire  Resistance 
(min) 

Calculated  Fire 
Resistance  (min) 

Normal  Weight 
Carbonate  Aggregate 

81 

60 

64  (+  6.7%) 

102 

90 

95  (+  5.6%) 

117 

120 

121  (+1.0%) 

145 

180 

179  (- 1.0%) 

168 

240 

235  (-  2.0%) 

Normal  Weight  Siliceous 
Aggregate 

89 

60 

67  (+  11.1%) 

109 

90 

127 

120 

157 

176  (-2.2%) 

178 

240 

220  (-  8.3%) 

Light  Weight  Expanded 
Shale 

64 

60 

64  (+  6.7%) 

79 

90 

93  (+  3.3%) 

91 

120 

120  (-) 

112 

180 

175  (-  2.8%) 

130 

240 

231  (-3.8%) 

Table  D-3  indicates  that  the  calculated  fire  resistance  is  within  eleven  percent  and  typically 
within  five  percent  or  less  of  the  actual  fire  resistance  rating. 

Sensitivity  of  Results  to  Input  Parameter  Assumptions 

A  sensitivity  analysis  of  the  material  properties  and  boundary  conditions  was  performed  to 
determine  critical  parameters.  Figure  D-8  shows  the  impact  of  uncertainty  in  the  thermal 
material  properties  of  the  concrete  for  the  carbonate  aggregate.  Figure  D-8  indicates  that  there  is 
a  considerable  range  in  the  predicted  fire  resistance  with  a  material  property  uncertainty  of 
twenty-five  percent.  This  uncertainty  is  within  the  range  of  data  for  various  concrete  aggregates 
[Harmathy,  1983;  Lie,  1972;  Flynn,  1999],  thus  identification  of  the  correct  aggregate  and 
thermal  properties  is  essential  when  modeling  concrete  materials  with  HEATING.  Similar 
results  are  observed  for  siliceous  and  expanded  shale  aggregates. 

Figure  D-9  shows  the  impact  of  boundary  condition  uncertainty  and  Figure  D-10  shows  the 
impact  of  variations  in  the  thermocouple  pad  material  properties  and  physical  dimensions.  Also 
shown  in  Figure  D-10  is  the  temperature  of  the  concrete  without  thermocouple  pads,  indicating  a 
significant  improvement  in  performance.  This  is  consistent  with  observations  of  temperature 
measurements  made  with  and  without  thermocouple  pads  that  demonstrated  an  increase  of  up  to 
30  minutes  in  the  fire  resistance  rating  and  an  overall  temperature  difference  of  44°C  between 
the  bare  and  covered  thermocouples  [Carlson  et  al.,  1962].  Collectively,  Figures  D-9  and  D-10 
suggest  that  the  boundary  conditions  and  the  material  properties  of  the  thermocouple  pad 
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(provided  it  is  included  in  the  model)  are  not  significant  parameters  when  modeling  the 
performance  of  concrete  with  HEATING. 


#  Actual  fire  resistance  rating  [ASCE  29-99, 1998] 

Predicted  fire  resistance  rating  (baseline) 

-  Concrete  thermal  conductivity  decreased  by  25% 

- Concrete  thermal  conductivdty  Increased  by  25% 

- Concrete  thermal  heat  capacity  decreased  by  25% 

- —  Concrete  thermal  heat  capacity  increased  by  25% 

.  Concrete  moisture  content  reduced  to  0% 

-  Concrete  moisture  content  Increased  to  3% 


Concrete  Thickness  (mm) 

Fig.  D-8  —  Sensitivity  of  Carbonate  Aggregate  Normal  Weight  Concrete  Predicted  Temperature  to 
Uncertainty  in  the  Thermal  Material  Property  Data 


#  Actual  fire  resistance  rating  [ASCE  29-99, 1998] 

- A -  Predicted  fire  resistance  rating  (baseline) 

-  Effective  surface  emissivity  decreased  to  0.7  on  all  surfaces 

- Effective  surface  emissivity  increased  to  0.95  on  all  surfaces 

- Convection  increased  by  25%  on  all  surfaces 

- Convection  decreased  by  25%  on  all  surfaces 


Concrete  Thickness  (mm) 

Fig.  D-9  —  Sensitivity  of  Carbonate  Aggregate  Normal  Weight  Concrete  Predicted 
Temperature  to  Uncertainty  in  the  Boundaiy  Condition  Parameters 
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•  Actual  fire  resistance  rating  [ASCE  29-99, 1998] 

-  Predicted  fire  resistance  rating  (baseline) 

-  TC  pad  den^  decreased  to  159  kg/m^  thickness  increased  to  1 1.1  mm 

- TC  pad  density  increased  to  51 9kg/m*;  thickness  decreased  to  7.9  mm 

- jC  pad  theimal  conductivity  decreased  to  0.049  W/m-"C 

- TC  pad  thermal  conductivity  Increased  to  0.057  W/m-*C 

-  Bare  thermocouple  (pad  removed) 


Concrete  Thickness  (mm) 

Fig.  D-IO  —  Sensitivity  of  Carbonate  Aggregate  Normal  Weight  Concrete  to  Variations  in 
Thermocouple  Pad  Dimensions  and  Thermal  Material  Properties 


3.3  Validation  Case  3:  Aluminum  Plate  Protected  with  Mineral  Fiber  Blanket 
Insulation  Exposed  to  ISO  6944  Test  Conditions 


Validation  Case  3  is  a  bulkhead/deck  assembly  that  involves  a  6  mm  thick  aluminum  plate 
and  Unifrax  Insulfrax  mineral  fiber  blanket  fire  insulation.  The  test  data  available  for  this 
assembly  includes  two  designs  tested  in  accordance  with  A.754(18)  [1993],  which  uses  the  ISO 
834  Standard  Time-Temperature  Profile  [IMO  A.754(18),  1993;  ISO  834, 1999]].  The  two 
assembly  ratings  are  A30  (30  minutes)  and  A60  (60  minutes).  The  material  properties  of  the 
Insulfrax  insulation  are  proprietary  to  Unifrax;  however  the  test  results  (exposure  temperature 
and  aluminum  temperature  versus  time)  and  the  rating  of  the  assembly  may  be  used  and  form  the 
comparison  basis  for  this  validation  case  [DIFST,  2002a;  DIFST,  2002b;  DIFST,  2002c,  and 
DIFST,  2002d]. 

The  specific  basis  for  the  validation  case  is  the  fire  resistance  rating  of  the  aluminum- 
insulation  assembly  and  the  transient  aluminum  temperature.  A  one-dimensional  thermal 
analysis  is  used  to  calculate  the  aluminum  temperature  as  a  function  of  time.  Figure  D-1 1 
depicts  the  configuration  modeled  by  HEATING.  Per  IMO  A.754(18)  [1993],  thermocouples 
used  to  measure  the  unexposed  aluminum  surface  temperature  are  covered  with  pads  in  an 
attempt  to  simulate  contact  with  a  combustible  material.  As  noted  in  Section  3.2,  the  presence  of 
thermocouple  pads  impacts  the  fire  resistance  rating  and  the  unexposed  surface  temperature  of 
the  assembly.  As  such,  they  are  included  in  the  thermal  model  used  in  this  validation  case. 
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Material  Properties 

The  material  properties  that  are  necessary  for  calculating  the  fire  resistance  rating  of  the 
Insulfrax-aluminum  bulkhead/deck  assembly  are  the  thermal  conductivity,  heat  capacity,  and 
density  of  all  materials  involved,  including  the  thermocouple  pads.  The  thermal  conductivity 
and  the  heat  capacity  of  the  mineral  fiber  blanket  were  provided  by  Insulffax  but  are  proprietary 
and  therefore  are  not  plotted  in  this  document.  The  density  for  each  material  is  as  follows. 

•  Aluminum:  2,787  kg/m^  [Holman,  1990]; 

.  Insulfrax:  96  kg/m^  [DIFST,  2002a;  DIFST,  2002c];  and 

.  Thermocouple  pad:  300  kg/m^  [ASTM  El  19-98, 1999]. 

•  The  heat  capacity  for  the  aluminum  and  thermocouple  pad  is  as  follows: 

•  Aluminum:  883  J/kg-°C  [Holman,  1990];  and 

•  Thermocouple  pad:  960  J/kg-°C,  assumed  (typical  mineral  fiber)  [Harmathy,  1983]. 
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The  thermal  conductivity  for  aluminum  is  a  function  of  temperature  as  shown  in  Table  D-4 
[Holman,  1990].  The  thermal  conductivity  for  the  thermocouple  pad  is  0.55  W/m-°C.  Note  that 
this  is  greater  than  the  value  prescribed  by  ASTM  El  19.  Due  to  the  considerable  uncertainty  in 
the  thermocouple  pad  thickness,  density,  and  thermal  conductivity,  this  parameter,  along  with  the 
emissivity  of  the  thermocouple  pad,  was  used  to  correlate  the  data  and  the  value  is  thus  not  likely 
the  true  thermal  conductivity. 

Table  D-4  —  Thermal  Conductivity  of  Aluminum  [Holman,  1990] 


Temperature  (°C) 

0 

159 

100 

182 

200 

194 

500 

194 

Boundary  Conditions 

A  forced  convection  boundary  condition  was  applied  to  the  exposed  surface  of  the 
bulkhead/deck  assembly  and  a  natural  convection  boundary  condition  was  applied  to  the 
unexposed  surface.  The  following  form  was  used  in  both  cases  [Holman,  1990]: 

(8) 

where  n  is  an  exponent  that  ranges  from  unity  (forced  convection)  to  about  1.33  (natural 
convection)  and  all  other  terms  have  been  defined  in  Section  3.1.  A  convection  coefficient  of 
15  W/m^-°C  and  an  exponent  of  unity  were  applied  to  the  exposed  surface  of  the  assembly  and  a 
convection  coefficient  of  1  wW-°C’  and  an  exponent  of  1.25  were  applied  to  the  unexposed 
thermocouple  pad  surface.  These  values  are  consistent  with  furnace  exposures  [Barnett,  1989; 
Issen,  1974]. 

A  thermal  radiation  boundary  condition  was  applied  to  the  exposed  and  unexposed  surfaces 
as  described  by  Equation  4.  The  emissivity  of  the  Insulfrax  surface  is  also  proprietary  and  is  thus 
not  provided.  The  emissivity  of  the  exposure  source  was  0.686  and  the  effective  emissivity  of 
the  ambient  thermocouple  pad  surface  was  0.53.  The  view  factor  was  unity  between  the  ambient 
environment  and  the  Insulfrax  and  the  thermocouple  pad  surfaces. 

The  exposure  fire  was  the  ISO  834  Standard  Time-Temperature  profile.  The  temperature  as  a 
function  of  time  is  as  follows  [ISO  834,  1999]: 

7>(0  =  7;+3451og8r  (9) 

where  is  the  initial  ambient  temperature  (°C)  and  t  is  the  time  in  minutes. 


D-25 


Results 


The  implicit  Crank-Nicolson  method  [Childs,  1998;  Smith,  1987;  Jaluria  et  al.,  1986]  was 
used  to  solve  for  the  aluminum  temperature.  The  tolerance  for  convergence  was  1X10^  and  the 
solution  time  step  as  1.0  second.  Decreasing  the  tolerance,  decreasing  the  time  step,  or  using  the 
Classical  Explicit  Procedure  or  Classical  Implicit  Procedure  resulted  in  a  less  than  3°C 
fluctuation  in  the  maximum  temperature  predictions.  Twenty-five  cells  were  used  in  the  Unifrax 
Insulfrax,  four  cells  were  used  in  the  aluminum,  and  five  cells  were  used  in  the  thermocouple 
pad.  The  maximum  temperature  predictions  changed  by  less  than  1°C  when  the  number  of  cells 
was  quadrupled  in  each  of  the  three  materials. 


The  thickness  of  the  insulation  for  the  A60  assembly  is  50  mm  [DIFST,  2002a;  DIFST, 
2002b]  and  for  the  A30  bulkhead/deck  assembly  it  is  38  mm  [DIFST,  2002c;  DIFST,  2002d]. 
The  calculated  average  aluminum  temperature  is  compared  to  the  measured  test  data  in  Figures 
D-12  and  D-13  for  the  A60  and  A30  assemblies.  The  figures  indicate  that  the  calculated  results 
are  in  agreement  with  the  test  data,  particularly  10-20  minutes  after  test  initiation. 


-  Average  unexposed  surface  temperature  increase  (measured) 

- Peak  unexposed  surface  temperature  increase  (measured) 

- Average  temperature  increase,  center  of  aluminum  sheet  (measured) 

— - — .  ^  Maximum  cold  face  temperature  increase  per  IMO  A.754(1 8) 

•  Predicted  aluminum  cold  face  temperature  increase 


Time  (minutes) 

Fig.  D-12  —  Predicted  versus  Measured  Aluminum  Temperature  for  A60  Bulkhead/Deck  Assembly  -  34  Cell  Model 


-  Peak  aluminum  temperature  increase  (measured) 

- Average  aluminum  temperature  inaease  (measured) 

- — ,  Maximum  aluminum  temperature  increase  per  IMO  A.754(18) 

•  Predicted  average  aluminum  temperature  increase 


Fig.  D- 13  —  Predicted  versus  Measured  Aluminum  Temperature  for  A30  Bulkhead/Deck  Assembly  -  34  Cell  Model 

The  calculated  time  to  reach  a  140°C  temperature  rise  on  the  unexposed  surface  is  compared 
in  Table  D-5.  The  rating  of  the  bulkhead/deck  assemblies  is  not  strictly  based  on  the  time  to 
reach  a  140°C  temperature  increase.  The  aluminum  temperature  must  not  exceed  200°C  at  60 
minutes  for  load  bearing  sections.  This  is  confirmed  in  Table  D-5  for  both  the  A60  and  A30 
bulkhead/deck  assemblies.  As  such,  HEATING  accurately  predicts  the  rating  of  the  assemblies 
relative  to  the  IMO  A.754(18)  test  criteria. 

Table  D-5  —  Time  Aluminum  Reaches  Critical  Temperature  for  Bulkhead/Deck 
Test  Assemblies  —  34  Cell  Model 


Test 

Measured  Time  to  140°C 
Aluminum  Temperature  Rise 
(minutes) 

Calculated  Time  to  140°C 
Aluminum  Temperature  Rise 
(minutes) 

A60 

60 

60 

A30 

44 

44 

Sensitivity  of  Results  to  Input  Parameter  Assumptions 

The  sensitivity  of  the  calculation  results  to  the  material  properties  and  boundary  conditions 
was  determined.  Figure  D-14  shows  the  impact  of  uncertainty  in  the  thermal  material  properties 
of  the  insulation  and  the  thermocouple  pad  for  the  A30  bulkhead/deck  assembly.  Figure  D-15 
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shows  the  sensitivity  of  the  temperature  predictions  to  the  boundary  condition  parameters  for  the 
A30  bulkhead/deck  assembly.  Similar  results  are  observed  for  the  A60  bulkhead/deck  assembly. 

Figures  D-14  and  D-15  suggest  that  there  are  two  critical  parameters:  the  thermal 
conductivity  of  the  insulation  material  and  the  effective  emissivity  of  the  exposed  surface.  The 
other  parameters  impact  the  results  by  less  than  five  percent. 


3.4  Validation  Case  4:  Steel  Protected  with  Mineral  Fiber  Insulation  Exposed  to  UL 
1709  Rapid  Temperature  Rise  Fire  Exposure 


Validation  Case  4  is  a  steel  plate  protected  with  Stucto-gard  mineral  fiber  insulation  that  is 
exposed  to  the  Underwriters  Laboratories  (UL)  1709  [1998]  Rapid  Temperature  Rise  fire 
exposure.  The  test  data  available  for  this  assembly  are  the  fire  exposure  temperature  as  a 
function  of  time  and  the  temperature  at  the  Structo-gard-steel  interface  as  a  function  of  time. 
The  total  test  duration  was  seventy  five  minutes. 

The  specific  basis  for  the  validation  case  is  transient  temperature  at  the  Structo-gard-steel 
interface  of  the  assembly.  A  one-dimensional  thermal  analysis  is  used  to  calculate  the 
temperature  profile.  Figure  D-16  depicts  the  configuration  modeled  by  HEATING.  Although 
the  UL  1709  test  requires  thermocouple  pads  on  the  unexposed  surface,  the  location  where  the 
interface  temperature  was  measured  was  not  near  areas  in  which  there  were  thermocouple  pads 
on  the  unexposed  steel  [Leonard  et  al.,  1993;  Soriathia  et  al.,  1998]. 


Predicted  aluminum  temperature  (baseline) 
insutfrax  thermal  conductivity  increased  by  25% 

Insulfrax  thermal  conductivity  decreased  by  25% 

Insulfirax  heat  capacity  increased  by  25% 

Insulfrax  heat  capacity  decreased  by  25% 

Thermocouple  pad  thermal  conductivity  increased  by  25% 
Thermocouple  pad  thermal  conductivity  decreased  by  25% 


Time  (minutes) 

Fig.  D-14  —  Sensitivity  of  the  A30  Bulkhead/Deck  Predicted  Temperature  to  Uncertainty  in  the  Material  Property  Data 
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•  Predicted  aluminum  temperature  (baseline) 

-  Exposed  surface  effective  emissivity  Increased  by  25% 

.  Exposed  surface  effective  emissivity  defreased  by  25% 

- Unexposed  surface  effective  emissh/ity  increased  by  25% 

-  Unexposed  surface  effective  emissivity  decreased  by  25% 

- Exposed  surface  convection  increased  by  25% 

-  Exposed  surface  convection  decreased  by  25% 

- Unexposed  surface  convection  increased  by  25% 

.  Unexposed  surface  convection  decreased  by  25% 


Time  (minutes) 

Fig.  D-15  —  Sensitivity  of  the  A30  Bulkhead/Deck  Predicted  Temperature  to  the  Uncertainty 
in  the  Boundary  Condition  Parameters 
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Fig.  D-16  —  Validation  Case  4  (Not  to  Scale) 


Material  Properties 

The  material  properties  necessary  for  this  validation  case  are  the  thermal  conductivity,  the 
heat  capacity,  and  the  density  of  the  Structo-gard  and  steel.  The  properties  for  structo-gard  were 
obtained  from  the  manufacturer  [Thermal  Ceramics,  2000]  and  the  properties  for  the  steel  were 
obtained  from  the  HEATING  material  library  database  for  carbon  steel  [Childs,  1998]. 

Figure  D-17  shows  the  thermal  conductivity  and  heat  capacity  for  the  insulation  and  steel  as  a 
function  of  temperature.  The  density  of  the  two  materials  is  as  follows; 

.  Steel  -  7,860  kg/m^  [Childs,  1998];  and 

•  Structo-Gard  -  102.5  kgW  [Thermal  Ceramics,  2000]. 
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Thermal  conductivity  (Structo-gard) 
Thermal  conductivity  (Carbon  Steel,  X100) 


Heat  capactiy  (Structo-gard) 
Heat  capacity  (Steel) 
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Fig,  D-17  —  Thermal  Material  Properties  of  Structo-gard  and  Steel 


Boundary  Conditions 

A  forced  convection  boundary  condition  was  applied  to  the  exposed  and  unexposed  surfaces 
of  the  Structo-gard  -  steel  assembly  as  described  by  Equation  8.  A  convection  coefficient  of 
15  W/m^-°C  was  used  for  both  surfaces. 

A  thermal  radiation  boundary  condition  was  applied  to  the  exposed  and  unexposed  surfaces 
as  described  by  Equation  4.  The  ambient  emissivity  was  assumed  to  be  0.95.  The  emissivity  of 
the  Structo-gard  was  assumed  to  be  0.95  and  the  emissivity  of  the  steel  was  assumed  to  be  0.8. 
As  such,  the  effective  emissivity  calculated  using  Equation  5  of  the  exposed  (Structo-gard) 
surface  was  0.9  and  the  effective  emissivity  of  the  unexposed  (steel)  surface  was  0.77.  The  view 
factor  was  unity  between  the  ambient  environment  (exposed  and  unexposed)  and  the  all  surfaces. 

The  exposure  fire  for  the  tests  was  the  UL  1709  as  shown  in  Figure  D-18.  The  actual  average 
temperature  data,  also  shown  in  Figure  D-18,  was  used  as  the  input  for  the  validation  model. 
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Fig,  D‘18  —  Fire  Exposure  for  Leonard  et  al.  [1993]  Test  #  1 


Results 

The  Classical  Implicit  Scheme  [Childs,  1998;  Smith,  1987;  Jaluria  et  al.,  1986]  was  used  to 
solve  for  the  unexposed  concrete  surface  temperature.  The  tolerance  for  convergence  was  1  X 
10'^  and  the  solution  time  step  was  2.0  seconds.  Decreasing  the  tolerance,  decreasing  the  time 
step,  or  using  the  Classical  Explicit  Procedure  or  Crank-Nicolson  Method  resulted  in  a  less  than 
1°C  fluctuation  in  the  maximum  temperature  predictions  at  the  Structo-gard  -  steel  interface. 
Eighty-  cells  were  used  in  the  Structo-gard  and  sixteen  cells  were  used  in  the  steel.  Increasing 
the  number  of  cells  in  either  material  did  not  alter  the  temperature  results. 

The  measured  temperature  at  the  Stucto-gard  -  steel  interface  is  compared  with  the 
calculated  temperature  in  Figure  D-19  for  the  first  thirty  minutes  of  Test  #1  in  Leonard  et  al. 
[1993].  The  results  are  conservative  {i.e.,  overestimate  the  interface  temperature)  and  are  within 
ten  percent  of  the  measured  data. 
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Fig.  D-19  —  Predicted  versus  Measured  Structo-gard  -  Steel  Interface  Temperature  (96  Cell  Model) 


Sensitivity  of  Results  to  Input  Parameter  Assumptions 

The  sensitivity  of  the  calculation  results  to  the  material  properties  and  boundary  conditions 
was  determined.  Figure  D-20  shows  the  impact  of  uncertainty  in  the  thermal  material  properties 
of  the  Structo-gard  insulation.  The  material  properties  of  steel  are  well  established  and  were  not 
varied.  Figure  D-21  shows  the  sensitivity  of  the  temperature  predictions  to  the  boundary 
condition  parameters  on  the  exposed  and  unexposed  surfaces. 

Figures  D-20  and  D-21  indicate  that  the  thermal  conductivity  of  the  Structo-gard  is  the  most 
critical  parameter.  The  other  parameters  impact  the  results  by  less  than  five  -  ten  percent. 


3.5  Validation  Case  5:  Unprotected  Steel  Column  Filled  with  Concrete  in  Web  Area 
Exposed  to  ASTME119  Test  Conditions 


Validation  Case  5  involves  fire  tests  on  two  unprotected  steel  W-shape  columns  with  the  area 
between  flanges  filled  with  concrete  [Underwriters  Laboratories,  1919].  The  columns  are 
exposed  to  the  ASTM  El  19  Standard  Time-Temperature  profile  and  the  steel  temperatures  are 
measured  at  several  locations  and  the  failure  time  of  the  columns  noted.  The  flange  area  for  the 
column  in  Test#  14  was  filled  with  carbonate  aggregate  normal  weight  concrete;  siliceous 
aggregate  normal  weight  concrete  was  used  in  the  same  manner  for  Test  #15. 
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Structo-gard  -  Steel  Interface 
Temperature  (®C) 


Predicted  (baseline) 

Structo-gard  thermal  conductivity  Increased  by  25% 
Structo-gard  thermal  conductivity  decreased  by  25% 
Structo-gard  heat  capacity  increased  by  25% 
Structo-gard  heat  capacity  decreased  by  25% 


Fig.  D-20  —  Sensitivity  of  Structo-gard  -  Steel  Interface  Predicted  Temperature  to  Uncertainty  in  the  Material  Property  Data 
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Predicted  (baseline) 

-  Exposed  surface  effective  emissivity  increased  to  1 .0 

.  Exposed  surface  effective  emissivity  decreased  to  0.8 

- Unexposed  surface  effective  emissivity  increased  to  0.85 

-  Unexposed  surface  effective  emissivity  decreased  to  0.65 

- Convection  Increased  by  25%  on  exposed/unexposed  surfaces 

-  Convection  decreased  by  25%  on  exposed/unexposed  surfaces 


Fig.  D-21  —  Sensitivity  of  Structo-gard  -  Steel  Interface  Predicted  Temperature  to  Uncertainty 
in  the  Boundary  Condition  Parameters 


The  specific  basis  for  the  validation  case  is  transient  temperature  at  various  locations  on  the 
steel  and  the  time  that  the  steel  columns  exceed  the  ASTM  El  19  critical  steel  temperatures.  A 
two-dimensional  thermal  analysis  is  used  to  calculate  the  temperature  profiles.  Figure  D-22 
depicts  the  configuration  modeled  by  HEATING.  The  ASTM  El  19  criteria  for  rating  steel  beams 
is  an  average  temperature  of  538°C  or  a  peak  temperamre  of  649°C,  whichever  occurs  first 
[ASTM  El  19-98, 1999].  The  peak  and  average  temperature  are  derived  fi-om  measurements 
taken  at  three  locations,  which  are  shown  in  Figure  D-22  [ASTM  El  19-98,  1999]. 
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Material  Properties 


The  material  properties  for  steel  are  as  described  for  Validation  Case  1.  The  material 
properties  for  the  concrete  are  as  described  for  Validation  Case  2  for  the  carbonate  and  siliceous 
aggregate  normal  weight  concretes. 

Boundary  Conditions 

The  column  was  modeled  as  a  quarter  section  due  to  symmetry,  as  shown  in  Figure  D-22.  An 
adiabatic  boundary  condition  was  thus  applied  to  all  symmetry  axes. 

A  forced  convection  boimdary  condition  was  applied  to  the  exposed  surfaces  of  the  steel  and 
concrete  as  described  by  Equation  8.  A  convection  coefficient  of  15  W/m^-°C  was  used  for  both 
surfaces. 

A  thermal  radiation  boundary  condition  was  applied  to  the  exposed  steel  and  concrete  as 
described  by  Equation  4.  The  emissivity  of  the  steel  was  assumed  to  be  0.8  and  the  concrete 
0.91  [Holman,  1990;  Sparrow  et  al.,  1978;  Incropera  et  al,  1985].  The  emissivity  of  the 
exposure  source  and  surroundings  was  determined  from  the  test  data  to  be  approximately  0.8; 
thus  the  effective  emissivity  for  the  steel  surface  was  0.63  and  for  the  concrete  surface  0.7.  The 
view  factor  was  unity  between  the  ambient  exposure  environment  and  the  all  surfaces. 

The  exposure  fire  for  the  tests  was  the  ASTM  El  19  Standard  Time-Temperature  Profile 
[1998]  as  shown  in  Figure  D-23.  The  actual  average  temperature  data  for  Test  #14  and  Test  #15, 
also  shown  in  Figure  D-23,  was  used  as  input  for  the  validation  model. 

Results 

Levy’s  Explicit  Procedure  [Childs,  1998;  Levy,  1968]  was  used  to  solve  for  the  temperature 
distribution  in  the  column  assemblies.  The  solution  time  step  is  determined  fi:om  the  stability 
criteria  and  a  multiplier,  which  was  set  to  2.0.  Decreasing  the  time  step  multiplier,  or  using  the 
implicit  Crank-Nicolson  method  or  the  Classical  Implicit  Procedure  resulted  in  maximum 
temperature  predictions  that  varied  by  less  then  4°C.  Thirty  cells  were  used  in  each  dimension  in 
the  concrete  and  five  cells  were  used  in  the  narrow  dimension  of  the  steel  for  a  total  of  1,225 
cells.  The  temperature  predictions  did  not  significantly  change  when  the  number  of  total  cells 
was  increased  to  5,041  with  increases  in  each  dimension  and  each  material. 

The  model  was  validated  by  comparing  the  transient  temperature  profile  at  three  locations  on 
the  steel  in  each  test  (shown  in  Figure  D-22).  Location  #1  corresponds  to  the  edge  of  the  flange; 
location  #2  corresponds  to  the  middle  of  the  flange;  and  location  #3  corresponds  to  the  middle  of 
the  web.  The  results  are  shown  in  Figures  D-24  and  D-25  for  Test  #14  and  Test  #15, 
respectively.  Figure  D-26  shows  the  temperature  contours  in  the  Test  #14  column  at  45  minutes, 
the  approximate  time  at  which  the  beam  exceeds  the  ASTM  El  19  critical  peak  temperature  of 
649°C. 
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Fig.  D-23  —  Measured  Exposure  Temperature  for  Test  #14  and  Test  #15 


Fig.  D-24  —  Predicted  versus  Measured  Temperatures  for  Test  #14  Column  -  1,225  Cell  Model 
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Fig.  D-25  —  Predicted  versus  Measured  Temperatures  for  Test  #15  Column  -  1,225  Cell  Model 


Fig.  D-26  —  Predicted  Isotherms  (®C)  in  Test  #14  Column  45  minutes  after  Test  Start  -  1,225  Cell  Model 
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The  time  to  exceed  the  ASTM  El  19  critical  steel  temperatures  for  columns  (viz.,  an  average 
temperature  of  538  °C  or  a  peak  temperature  of  649°C,  whichever  occurs  first  [ASTM  El  19, 
1998])  was  also  used  as  a  validation  criterion.  The  average  temperature  is  based  on  four 
locations  per  ASTM  El  19  [1998]:  three  on  the  upper  and  lower  flanges  and  one  in  the  web.  The 
predicted  time  to  exceed  the  ASTM  El  19  temperature  criteria  is  compared  with  the  measured 
time  in  Table  D-6. 

Table  D-6  —  Calculated  and  Predicted  Time  to  Exceed  ASTM  El  19  Critical  Steel  Temperatures 

for  Columns  -  1,225  Cell  Model 


Test 

Time  to  Exceed  Average  Temperature  of 
538°C  (min) 

Time  to  Exceed  Peak  Temperature  of 
649°C  (min) 

Measured 

Predicted 

Measured 

Predicted 

14 

49 

_ 49B _ 

47 

45  (+4.4%) 

15 

38 

44 

37  (-  16%) 

Table  D-6  indicates  that  HEATING  predicted  the  average  and  peak  temperature  for  Test  #14. 
The  predicted  fire  resistance  rating  for  this  test  was  45  minutes,  as  compared  to  47  minutes 
measured.  The  average  temperature  prediction  for  Test  #15  was  also  predicted  well,  however  the 
peak  temperature  was  over-predicted  in  the  model.  The  predicted  fire  resistance  rating  for  the 
assembly  is  37  minutes,  as  compared  to  38  minutes  measured.  Note  that  the  limiting  criteria  is 
different  in  the  predicted  and  measured  cases  for  Test  #15,  however  the  predicted  versus 
measured  fire  resistance  rating  is  nearly  identical. 

Sensitivity  of  Results  to  Input  Parameter  Assumptions 

The  sensitivity  of  the  calculation  results  to  the  material  properties  and  boundary  conditions 
was  investigated  for  Test  #14,  Location  #1.  Similar  results  are  observed  for  other  locations  and 
for  Test  #15. 

The  sensitivity  of  the  transient  temperature  of  Test  #14,  Location  #1  to  uncertainty  in 
material  property  data  is  shown  in  Figure  D-27.  Likewise,  Figure  D-28  shows  the  sensitivity  of 
the  transient  temperature  at  the  same  location  to  boundary  condition  uncertainty.  Table  D-7 
shows  the  variation  in  the  predicted  time  for  the  Test  #14  column  to  reach  an  average 
temperature  of  538°C  and  a  peak  temperature  of  649°C  as  a  result  of  material  property  and 
boundary  condition  uncertainty. 

Figures  D-27,  D-28  and  Table  D-7  indicate  that  the  column-concrete  system  is  not  especially 
sensitive  to  uncertainty  in  the  material  property  or  boundary  condition  parameters.  The  effective 
emissivity  may  impact  the  maximum  predicted  temperatures  by  plus  or  minus  50°C;  the 
remaining  parameters  have  significantly  less  of  an  impact. 
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— Predicted  tempeature  at  Test  #14,  Location  #1  (baseline) 

-  Concrete  thermal  conductivity  Increased  by  25% 

.  Concrete  thermal  conductivity  decreased  by  25% 

- Concrete  heat  capacity  increased  by  25% 

-  Concrete  heat  capacity  decreased  by  25% 

- Concrete  moisture  content  reduced  to  0  percent 

-  Concrete  moisture  content  Increased  to  3  percent 


Time  (minutes) 

Fig.  D-27  —  Sensitivity  of  Test  #14,  Location  #1  Temperature  Predictions  to  Material  Property  Parameter  Uncertainty 


Time  (minutes) 

Fig.  D-28  —  Sensitivily  of  Test  #14,  Location  #1  Temperature  Predictions  to  Boundary  Condition  Parameter  Uncertainty 
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Table  D-7  —  Sensitivity  of  the  Predicted  Time  to  Exceed  ASTM  El  19  Critical  Steel  Temperatures  to 
Material  Property  and  Boundary  Condition  Parameter  Uncertainty 


Sensitivity  Parameter  Modification 

Predicted  time  to  Reach 
538°C  Average  (min) 

Predicted  time  to  Reach 

649°C  Peak 

(min) 

No  modifications  (baseline) 

49.4 

45.5 

Concrete  conductivity  increased  by  25% 

46.1  (+ 1.3%) 

Concrete  conductivity  decreased  by  25% 

48.5  (- 1.8%) 

44.7  (- 1.8%) 

Concrete  heat  capacity  increased  by  25% 

52.1  (-1-  5.5%) 

47.4  (+4.2%) 

Concrete  heat  capacity  increased  by  25% 

41.1  (-9.7%) 

Concrete  moisture  reduced  to  0% 

44.4  (-  2.4%) 

Concrete  moisture  increased  to  3% 

51.1  (+3.4%) 

46.4  (+  2.0%) 

Effective  emissivity  increased  by  ~25% 

45.2  (-  8.5%) 

Effective  emissivity  decreased  by  ~25% 

54.8  (+ 10.9%) 

5 1.9  (+14.1%) 

Convection  increased  by  25% 

48.5  (- 1.8%) 

44.5  (-  2.2%) 

Convection  decreased  by  25% 

50.3  (+  1.8%) 

46.4  (+2.0%) 

3.6  Validation  Case  6:  Steel  Beam  Protected  with  Blazeshield  Exposed  to  ASTME119 
Test  Conditions  (UL  Design  J907) 


Validation  Case  6  is  Underwriters  Laboratories  (UL)  Design  Number  P907  [UL,  2003].  The 
UL  “P”  designs  feature  roof  or  ceiling  assemblies  that  are  unprotected,  however  the  steel  framing 
that  supports  the  assemblies  is  typically  protected  with  insulation  or  a  membrane.  UL  Design 
No.  P907  involves  a  W8X28  beam  that  supports  fluted  galvanized  steel-form  insulating  concrete 
(i.e.,  lightweight  concrete)  floor  units.  The  UL  Design  No.  P907  roof-ceiling  assembly  is  rated 
for  2  hours  fire  resistance  when  exposed  to  the  ASTM  El  19  Standard  Time-Temperature  profile 
provided  the  steel  is  protected  with  a  31.7  mm  thick  coating  of  Blazeshield  D-C/F  sprayed  fire 
resistive  material  [UL,  2003]. 

The  specific  basis  for  the  validation  case  is  the  unrestrained  beam  fire  resistance  rating  for 
the  ceiling-roof  assembly  using  the  ASTM  El  19  critical  steel  temperature  criteria  [ASTM  El  19- 
98,  1999].  A  two-dimensional  thermal  analysis  is  used  to  calculate  the  temperature  profiles 
within  the  assembly.  Figure  D-29  depicts  the  configuration  that  was  modeled  with  HEATING. 
The  ASTM  El  19  criteria  for  rating  steel  beams  is  an  average  temperature  of  593°C  or  a  peak 
temperature  of  704°C,  whichever  occurs  first  [ASTM  El  19-98,  1999].  The  peak  and  average 
temperature  are  derived  from  measurements  taken  at  four  locations,  which  are  shown  in  Figure 
D-29  [ASTM  El  19-98,  1999]. 

Material  Properties 

The  material  properties  for  lightweight  concrete  are  as  described  for  the  lightweight 
expanded  shale  aggregate  concrete  used  in  Validation  Case  2  in  Section  3.2.  The  material 
properties  for  steel  are  as  described  for  Validation  Case  1  in  Section  3.1. 
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THERMAL  RADIATION  =  Ge,6(Tc‘-T.‘^ 


CONVECTION  =  hjCTc-TJ 


ASTM  E119  FIRE  EXPOSURE.  T^t) 

Fig.  D-29  —  Validation  Case  6  (Not  to  Scale) 


The  thermal  conductivity  and  heat  capacity  of  Blazeshield  D-C/F  vary  with  temperature  and 
are  shown  in  Figure  D-30  [Harmathy,  1983].  The  insulation  density  is  309  kg/m^  [Harmathy, 
1983]. 

Boundary  Conditions 

A  forced  convection  boundary  condition  was  applied  to  all  exposed  and  unexposed  surfaces 
of  the  concrete  and  Blazeshield  D-C/F  insulation  as  described  by  Equation  8.  A  convection 
coefficient  of  15  kW/m^  was  assumed  on  all  surfaces  where  convection  was  applied. 
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Fig.  D-30  —  Thennal  Conductivity  and  Heat  Capacity  of  Blazeshield  D-C/F  [Harmathy,  1983] 


A  thermal  radiation  boundary  condition  was  applied  to  the  exposed  and  unexposed  surfaces 
as  described  by  Equation  4.  The  emissivity  of  the  insulation  was  assumed  to  be  0.9  and  the 
concrete  emissivity  was  assumed  to  be  0.91,  as  described  for  Validation  Case  1  [Holman,  1990; 
Sparrow  et  ai,  1978;  Incropera  et  al.,  1985].  The  emissivity  of  the  exposed  and  unexposed 
ambient  environment  was  assumed  to  be  0.95.  Thus  the  effective  emissivity  for  the  exposed 
surfaces  was  0.86  and  for  imexposed  surfaces  0.87  as  calculated  using  Equation  5.  The  thermal 
radiation  configuration  factor  was  assumed  unity  on  all  surfaces  except  for  those  associated  with 
the  inner  flange  and  the  web  areas.  The  following  equation  was  used  to  calculate  the 
configuration  factor  between  the  web  and  the  exposure  source  [Howell,  1982]: 


F  = 


h 

w 


(10) 


where  F  is  the  configuration  factor  between  the  web  surface  and  the  exposure  source  (i.e., 
fraction  of  radiant  energy  incident  on  the  web  surface  originating  from  the  exposure  source),  h  is 
the  exposed  height  of  the  web  surface,  equal  to  the  web  height  minus  the  insulation  thickness  on 
the  upper  and  lower  flanges  (m),  and  w  is  the  exposed  length  of  the  inner  flange,  roughly  equal  to 
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one-half  the  flange  width  (m).  Likewise,  the  configuration  factor  between  the  inner  flange 
surfaces  and  the  exposure  source  is  given  by  the  following  [Howell,  1983]: 


F  =  0.5 


(11) 


where  F  is  the  configuration  factor  between  the  inner  flange  surfaces  and  the  exposure  source,  s 
is  the  exposed  height  of  the  web  surface,  equal  to  the  web  height  minus  the  insulation  thickness 
on  the  upper  and  lower  flanges  (m),  and  d  is  the  exposed  length  of  the  inner  flange,  roughly 
equal  to  one-half  the  flange  width  (m).  The  resulting  configuration  factors  for  a  W8X28  beam 
protected  with  3 1.7  mm  of  insulation  were  0.323  between  the  web  surfaces  and  the  exposure 
source  and  0.339  between  the  inner  flange  surfaces  and  the  exposure  source. 

The  use  of  equations  10  and  1 1  to  calculate  a  configuration  factor  between  the  surface  and 
the  exposure  source  implies  that  the  insulation  surfaces  on  the  unexposed  web  in  flange  areas  are 
approximately  at  the  same  temperature.  As  such  there  would  be  no  net  radiation  heat  exchange 
between  these  surfaces. 

The  exposure  source  temperature  is  the  ASTM  El  19-98  [1999]  given  by  Equation  6  in 
Section  3.2. 

Results 

The  implicit  Crank-Nicolson  method  [Childs,  1998;  Smith,  1987;  Jaluria  et  al.,  1986]  was 
used  to  solve  for  the  temperature  distribution  in  the  beam  assembly.  The  tolerance  for 
convergence  was  1  X  10"^  and  the  solution  time  step  was  2.0  seconds.  Decreasing  the  tolerance, 
decreasing  the  time  step,  or  using  the  Classical  Explicit  Procedure  or  Classical  Implicit 
Procedure  did  not  result  in  significant  changes  in  the  temperature  predictions.  Typically,  four  or 
five  cells  were  used  in  the  insulation  through  the  dimension  with  the  greatest  temperature 
gradient.  The  total  number  of  cells  was  1,064.  The  flutes  in  the  galvanized  steel  used  to  form  the 
concrete  were  simulated  as  a  1 10  mm  by  56  mm  adiabatic  void  positioned  above  the  center  of 
the  beam.  The  temperature  predictions  and  predicted  fire  resistance  rating  did  not  significantly 
change  when  the  number  of  total  cells  was  increased  to  4,024,  with  resolution  increases  in  each 
dimension  and  in  each  material. 

The  model  was  validated  by  comparing  the  calculated  fire  resistance  rating  using  the  ASTM 
El  19-98  temperature  criteria  for  steel  beams  to  the  listed  fire  resistance  rating  for  UL  Design  No. 
P907  with  3 1 .7  mm  Blazeshield  D-C/F  insulation.  The  listed  fire  resistance  rating  for  this 
assembly  is  2  hours.  The  predicted  fire  resistance  is  as  follows: 

Predicted  time  average  temperature  exceeds  593°C:  123.5  minutes 

Predicted  time  peak  temperature  exceeds  704°C:  132.8  minutes 

Thus,  the  predicted  fire  resistance  of  the  assembly  is  123.5  minutes,  3.5  minutes,  or  about 
three  percent,  greater  than  the  actual  listing  for  the  assembly  modeled. 
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The  predicted  temperature  at  the  four  locations  used  to  estimate  the  fire  resistance  rating  is 
shown  in  Figure  D-3 1.  Figure  D-32  shows  the  predicted  temperature  distribution  within  the 
entire  two-dimensional  model  domain  at  120  minutes  after  start  of  ASTM  El  19  temperature 
exposure. 


Time  (minutes) 

Fig.  D-3 1  —  Predicted  Temperature  at  Locations  Used  to  Obtain  Fire  Resistance 
Rating  for  UL  Design  No.  P907  -  1 ,064  Cell  Model 


Fig.  D-32  —  Calculated  Temperature  Distribution  (°C)  in  UL  Design  No.  P907  at  120  minutes 
after  Start  of  Exposure  -  1,064  Cell  Model 
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Sensitivity  of  Results  to  Input  Parameter  Assumptions 


The  sensitivity  of  the  calculation  results  to  uncertainty  in  the  material  properties  and 
boundary  conditions  is  assessed  relative  to  the  predicted  fire  resistance  rating  of  the  floor  or 
ceiling-roof  assembly.  The  material  properties  of  the  Blazeshield  D-C/F  and  concrete  were 
varied  by  plus  or  minus  25  percent.  Similarly,  the  radiation  and  convection  boundary  conditions 
were  increased  and  decreased  in  order  to  determine  the  change  in  the  predicted  fire  resistance 
rating.  The  results  are  summarized  in  Table  D-8. 

Table  D-8  —  Sensitivity  of  the  Calculated  Fire  Resistance  of  UL  Design  No.  P907  to  Uncertainty  in  the 
Material  Property  Data  and  Boundary  Condition  Parameters 


Modified  Parameter 

Average  temperature  of  593 °C 

Peak  Temperature  of  704°C  j 

Time  (min) 

■E5SH 

Time  (min) 

Predicted  (unmodified  baseline) 

123.5 

- 

132.8 

-  1 

Blazeshield  D-C/F  conductivity 
increased  by  25% 

106.9 

-  16.6 
(- 13.4%) 

114.0 

Blazeshield  D-C/F  conductivity 
decreased  by  25% 

149.8 

162.9 

Blazeshield  D-C/F  heat  capacity 
increased  by  25% 

129.3 

+  5.8 
(+4.7%) 

138.9 

+  6.1 
(+4.6%) 

Blazeshield  D-C/F  heat  capacity 
increased  by  25% 

117.8 

mSmIM 

126.7 

-6.1 

(-4.6%) 

Concrete  conductivity  increased  by 
25% 

123.7 

WaW 

132.8 

_ B _ 

Concrete  conductivity  decreased  by 
25% 

123.4 

-0.1 

(-0.1%) 

132.8 

(-) 

Effective  emissivity  of  insulation 
and  concrete  increased  by  -25% 

122.8 

132.1 

Effective  emissivity  of  insulation 
and  concrete  decreased  by  -25% 

124.7 

mSmIM 

133.8 

BiiiB 

Convection  increased  by  25% 

123.3 

132.7 

-0.1 

(-0.1%) 

Convection  increased  by  25% 

123.8 

+  0.3 
(+  0.2%) 

133.0 

+  0.2 
(+  0.2%) 

Table  D-8  indicates  that  the  calculated  fire  resistance  rating  for  UL  Design  No.  P907  is 
relatively  insensitive  to  uncertainty  in  the  material  properties  and  boundary  condition  parameters 
with  the  notable  exception  of  the  thermal  conductivity  of  the  Blazeshield  D-C/F  insulation 
material.  In  this  case,  a  twenty  five  percent  increase  or  decrease  in  the  thermal  conductivity 
value  could  result  in  a  15  -  25  percent  increase  or  decrease  in  the  predicted  fire  resistance  rating 
using  the  ASTM  El  19-98  [1999]  critical  steel  temperature  criteria. 
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3.7  Validation  Case  7:  Steel  Beam  Protected  with  Blazeshield  Exposed  to  ASTME119 
Test  Conditions  (UL  Design  P809) 


Validation  Case  7  is  Underwriters  Laboratories  (UL)  Design  Number  J809  [UL,  2003].  The 
UL  “J”  designs  feature  conerete  floor,  roof,  or  ceiling  assemblies  where  the  concrete  and  steel,  if 
present,  is  protected  with  a  sprayed  fire  resistive  material.  UL  Design  No.  J809  involves  a 
W8X28  beam  that  supports  normal  weight  siliceous  or  carbonate  aggregate  concrete  floor  or 
roof-ceiling  assembly.  The  UL  Design  No.  J809  floor  or  roof-ceiling  assembly  is  rated  for  2 
hours  fire  resistance  when  exposed  to  the  ASTM  El  19  Standard  Time-Temperature  profile 
provided  the  steel  is  protected  with  a  25.4  mm  thick  coating  of  Blazeshield  D-C/F  sprayed  fire 
resistive  material  and  the  exposed  portions  of  the  concrete  is  protected  with  a  14.2  mm  thick 
coating  of  Blazeshield  D-C/F  [Underwriters  Laboratories,  2003].  The  concrete  thickness  for  this 
particular  rating  and  insulation  thickness  is  88.9  mm  [Underwriters  Laboratories,  2003]. 

The  specific  basis  for  the  validation  case  is  the  unrestrained  beam  fire  resistance  rating  for 
the  ceiling-roof  assembly  using  the  ASTM  El  19  critical  steel  temperature  criteria  [ASTM  El  19- 
98, 1999].  A  two-dimensional  thermal  analysis  is  used  to  calculate  the  temperature  profiles 
within  the  assembly.  Figure  D-33  depicts  the  configuration  that  was  modeled  with  HEATING. 
The  ASTM  El  19  criteria  for  rating  steel  beams  is  an  average  temperature  of  593°C  or  a  peak 
temperature  of  704°C,  whichever  occurs  first  [ASTM  El  19-98,  1999].  The  peak  and  average 
temperature  are  derived  from  measurements  taken  at  four  locations,  which  are  shown  in  Figure 
D-33  [ASTM  El  19-98,  1999]. 
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CONVECTION  =  h^CTc-TJ 

THERMAL  RADIATION  =  Ee26CTc'-TJ) 


ASTM  E1 1 9  FIRE  EXPOSURE,  Tp(t)  ) 


Fig.  D-33  —  Validation  Case  7  (Not  to  Scale) 


Material  Properties 

The  material  properties  for  the  steel  and  Blazeshield  D-C/F  sprayed  fire  resistive  material  are 
as  described  in  Section  3.6  for  Validation  Case  6.  The  material  properties  for  siliceous  aggregate 
normal  weight  are  as  described  in  Section  3.2  for  Validation  Case  2. 

Boundary  Conditions 

The  boundary  conditions  are  as  described  for  Validation  Case  6  with  the  exception  of  the 
radiation  configuration  factor.  Because  there  is  6.3  mm  less  Blazeshield  D-C/F  insulation  on  the 
steel  than  Validation  Case  6,  the  radiation  configuration  factors,  calculated  using  Equations  10 
and  1 1,  are  different  firom  those  calculated  in  Section  3.6.  The  resulting  configuration  factors  for 
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a  W8X28  beam  protected  with  25.4  mm  of  insulation  are  0.28 1  between  the  web  surfaces  and  the 
exposure  source  and  0.36  between  the  iimer  flange  surfaces  and  the  exposure  source. 

Results 

The  implicit  Crank-Nicolson  method  [Childs,  1998;  Smith,  1987;  Jaluria  etal.,  1986]  was 
used  to  solve  for  the  temperature  distribution  in  the  beam  assembly.  The  tolerance  for 
convergence  was  1X10'®  and  the  solution  time  step  was  2.0  seconds.  Decreasing  the  tolerance, 
decreasing  the  time  step,  or  using  the  Classical  Explicit  Procedure  or  Classical  Implicit 
Procedure  did  not  result  in  significant  changes  in  the  temperature  predictions.Typically,  four  or 
five  cells  were  used  in  the  insulation  through  the  dimension  with  the  greatest  temperature 
gradient.  The  total  number  of  cells  was  1,436.  The  temperature  predictions  and  predicted  fire 
resistance  rating  did  not  significantly  change  when  the  number  of  total  cells  was  increased  to 
5,529,  with  resolution  increases  in  each  dimension  and  in  each  material. 

The  model  was  validated  by  comparing  estimating  the  fire  resistance  rating  using  the  ASTM 
El  19-98  [1999]  temperature  criteria  for  steel  beams  to  the  listed  fire  resistance  rating  for  UL 
Design  No.  J809  with  25.4  mm  Blazeshield  D-C/F  insulation  on  the  steel  beam  and  14.7  mm  on 
the  exposed  portions  of  the  concrete  floor  or  ceiling-roof  slab.  The  listed  fire  resistance  rating 
for  this  assembly  is  2  hours.  The  predicted  fire  resistance  is  as  follows: 

Predicted  time  average  temperature  exceeds  593 °C:  130.5  minutes 

Predicted  time  peak  temperature  exceeds  704°C:  1 1 6.3  minutes 

Thus,  the  predicted  fire  resistance  of  the  assembly  is  1 16.7  minutes,  3.7  minutes,  or  about 
three  percent,  less  than  the  actual  listing  for  the  assembly  modeled. 

The  predicted  temperature  at  the  four  locations  used  to  estimate  the  fire  resistance  rating  is 
shown  in  Figure  D-34.  Figure  D-35  shows  the  predicted  temperature  distribution  within  the 
entire  two-dimensional  model  domain  at  120  minutes  after  start  of  ASTM  El  19  temperature 
exposure. 
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Fig.  D-34  —  Predicted  Temperature  at  Locations  Used  to  Obtain  Fire  Resistance  Rating 
for  UL  Design  No.  J809  -  1,436  Cell  Model 
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Fig.  D-35  —  Calculated  Temperature  Distribution  (°C)  in  UL  Design  No.  J809  at  120  minutes 
after  Start  of  Exposure  -  1,436  Cell  Model 


Sensitivity  of  Results  to  Input  Parameter  Assumptions 

The  sensitivity  of  the  calculation  results  to  uncertainty  in  the  material  properties  and 
boundary  conditions  is  assessed  relative  to  the  predicted  fire  resistance  rating  of  the  floor  or 
ceiling-roof  assembly.  The  material  properties  of  the  Blazeshield  D-C/F  and  concrete  were 
varied  by  plus  or  minus  25  percent.  Similarly,  the  radiation  and  convection  boundary  conditions 
were  increased  and  decreased  in  order  to  determine  the  change  in  the  predicted  fire  resistance 
rating.  The  results  are  summarized  in  Table  D-9. 
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Table  D-9  —  Sensitivity  of  the  Calculated  Fire  Resistance  of  UL  Design  No.  J809  to 
Uncertainty  in  the  Material  Property  Data  and  Boundary  Condition  Parameters 


Modified  Parameter 

1  Average  temperature  of  593°C 

1  Peak  Temperature  of  704°C  | 

Time  (min) 

Time  (min) 

Predicted  (unmodified  baseline) 

130.5 

116.3 

-  1 

Blazeshield  D-C/F  conductivity 
increased  by  25% 

111.3 

99.3 

BHIIIl 

Blazeshield  D-C/F  conductivity 
decreased  by  25% 

163.2 

145.0 

+  28.7 
(+  24.7%) 

Blazeshield  D-C/F  heat  capacity 
increased  by  25% 

134.5 

120.5 

Blazeshield  D-C/F  heat  capacity 
increased  by  25% 

126.5 

■|||K2^l||||| 

112.2 

■Kil 

Concrete  conductivity  increased  by 
25% 

133.5 

■iiiBgQmi 

116.8 

^^D^IBi 

Concrete  conductivity  decreased  by 
25% 

127.0 

115.8 

-0.5 
(-  4.2%) 

Effective  emissivity  of  insulation 
and  concrete  increased  by  -25% 

128.8 

■mi 

115.2 

Effective  emissivity  of  insulation 
and  concrete  decreased  by  -25% 

133.5 

■■ngQiiiiillll 

118.7 

+  2.4 
(+2.1%) 

Convection  increased  by  25% 

130.3 

116.2 

-0.1 
(-  0.1%) 

Convection  increased  by  25% 

130.8 

+  0.3 
(+  2.2%) 

116.7 

+  0.4 
(+0.3%) 

Table  D-9  indicates  that  the  calculated  fire  resistance  rating  for  UL  Design  No.  J809  is 
relatively  insensitive  to  uncertainty  in  the  material  properties  and  boundary  condition  parameters 
with  the  notable  exception  of  the  thermal  conductivity  of  the  Blazeshield  D-C/F  insulation 
material.  In  this  case,  a  twenty  five  percent  increase  or  decrease  in  the  thermal  conductivity 
value  could  result  in  a  15  -  25  percent  increase  or  decrease  in  the  predicted  fire  resistance  rating 
using  the  ASTM  El  19-98  [1999]  critical  steel  temperature  criteria. 


3.8  Validation  Case  8:  Gypsum-Steel  Stud  Wall  Assembly  Exposed  to  CAN/ULL-SIOI- 
M89  Test  Conditions 


Validation  Case  8  consists  of  an  uninsulated  steel-stud  wall  assembly  constructed  and  tested 
in  accordance  with  CAN/ULC-S101-M89  [1989;  Sultan,  1996;  Alfawakhiri  et  al,  2000].  The 
wall  assembly  of  interest  consists  of  90  mm  by  30  mm  by  0.46  mm  thick  steel  studs  spaced  600 
mm  on  center.  Type  X  gypsum  15.9  mm  was  attached  to  the  exposed  and  unexposed  sides  of  the 
steel  studs.  The  wall  assembly  was  exposed  to  the  CAN/ULC-S101-M89  time-temperature 
profile  until  the  average  temperature  of  the  imexposed  wall  surface  (measured  using 
thermocouple  pads)  increased  to  140°C  above  ambient  [CAN/ULC-S101-M89, 1989;  Sultan, 
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1996],  the  fire  resistance  rating  of  the  assembly.  This  occurred  at  54  minutes  for  the  test  used  to 
validate  HEATING. 

The  specific  basis  for  the  validation  case  is  the  transient  temperature  within  the  cavity  space 
on  the  exposed  and  unexposed  gypsum  surfaces  and  the  average  temperature  on  the  unexposed 
side  of  the  wall.  The  latter  measurements  were  made  with  and  without  thermocouple  pads,  thus 
both  measurements  techniques  are  used  in  this  validation  case.  The  fire  resistance  rating  of  the 
assembly  was  also  used  as  a  validation  basis,  defined  by  CAN/ULC-S101-M89  [1989]  as  an 
average  temperature  increase  on  the  unexposed  surface  of  140°C  or  a  peak  single  point 
temperature  increase  of  180°C.  Both  measurements  include  thermocouple  pads. 

A  two-dimensional  thermal  analysis  is  used  to  calculate  the  temperatures  at  the  various 
locations  as  well  as  determining  the  fire  resistance  rating  of  the  wall  assembly.  A  single  wall 
cavity  section  was  modeled  due  to  symmetry.  Heat  transfer  within  the  void  space  was  included 
in  the  model.  Figure  D-36  depicts  the  configuration  modeled  where  there  is  no  thermocouple 
pad.  The  same  configuration  was  also  modeled  with  a  thermocouple  pad  present  as  noted  above. 
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Material  Properties 

The  material  properties  for  steel  and  the  thermocouple  pad  are  as  described  for  Validation 
Case  1. 

The  thermal  conductivity,  heat  capacity,  and  density  of  gypsum  vary  with  temperature 
[Harmathy,  1993;  Sultan,  1996;  Alfawakhiri,  2000].  Figure  D-37  shows  the  thermal 
conductivity  and  heat  capacity  of  Type  X  gypsum  as  a  function  of  temperature  [Sultan,  1996; 
Alfawakhiri,  2000].  The  density  is  as  follows  [Sultan,  1996]: 

•  698  kg/m^  -  20°C  to  80°C;  and 

•  576  kg/m^  -  over  80°C. 
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Fig,  D-37  —  Thermal  Conductivity  and  Heat  Capacity  of  Type  X  Gypsum  [Sultan,  1996;  Alfawakhiri,  2000], 


Boundary  Conditions 

A  natural  convection  boundary  condition  was  applied  to  all  exposed  and  unexposed  wall 
surfaces,  as  described  by  Equation  3,  except  for  inside  the  cavity  and  the  axes  of  symmetry, 
which  are  by  definition  adiabatic.  A  forced  convection  boundary  condition  was  applied  inside 
the  cavity.  The  natural  convection  coefficient  and  exponent  for  the  non-cavity  surfaces  was  1.31 
WW-°C’  and  1.33,  respectively  [Holman,  1990],  as  described  in  Section  3.2  for  wall  surface 
boundary  conditions.  The  forced  convection  coefficient  applied  to  the  cavity  boundaries  was 
derived  from  natural  convection  information  provided  by  Sultan  [1996].  Based  on  this 
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information,  the  equivalent  forced  convection  coefficient,  assuming  a  100°C  temperature 
difference  within  the  cavity  space,  a  cavity  height  of  3.05  m,  and  a  cavity  width  of  90  mm,  is  7.3 
W/m^-°C.  Trelles  et  al.  [2003]  describes  the  manner  in  which  convection  boundary  conditions 
are  applied  to  cavity  space  boundaries  using  HEATING. 

A  thermal  radiation  boundary  condition  was  applied  to  all  exposed  and  xmexposed  surfaces, 
including  inside  the  cavity,  of  the  gypsum  and  steel  wall  assembly  as  described  by  Equations  4 
and  5,  except  for  the  axes  of  symmetry,  which  are  by  definition  adiabatic.  The  emissivity  of  the 
steel  was  assumed  to  be  0.8,  the  emissivity  of  the  thermocouple  pad  was  assumed  to  be  0.91,  and 
the  emissivity  of  the  gypsum  was  assumed  to  be  0.84  [Sparrow  et  al.,  1978;  Holman,  1990;  and 
Incropera  et  al.,  1985;  Sultan,  1996].  The  ambient  emissivity  was  assumed  to  be  0.95.  The 
resultant  or  effective  emissivity  for  the  exposed  surface,  as  calculated  using  Equation  5,  is  thus 
0.8  and  for  the  unexposed  surface  it  is  0.87  (thermocouple  pad  present)  or  0.8  (no  thermocouple 
pad).  The  view  factor  was  unity  between  ambient  environment  and  the  concrete  and  the 
thermocouple  pad  surfaces.  The  effective  emissivity  and  view  factor  vary  from  node  to  node 
inside  the  wall  cavity.  Refer  to  Trelles  et  al.  [2003]  for  a  description  of  the  calculation  procedure 
for  thermal  radiation  heat  transfer  in  cavity  spaces. 

The  exposure  temperature  was  the  Standard  CANAJLC-S101-M89  [1989]  Time-Temperature 
profile,  a  similar  exposure  to  the  ASTM  El  19-98  [1999]  Standard  Time-Temperature  profile. 

The  CAN/ULC-S101-M89  temperature  profile  maybe  determined  using  the  following  equation 
[Sultan,  1996]: 


Tf(t)  =  T„  +750[l-exp{-3.79553V7)J+170.4l77  (12) 

where  T/t)  is  the  exposure  temperature  (°C)  at  time  t  (hr). 


Results 

The  implicit  Crank-Nicolson  method  [Childs,  1998;  Smith,  1987;  Jaluria  et  al.,  1986]  was 
used  to  solve  for  the  temperature  distribution  in  the  wall  assembly.  The  tolerance  for 
convergence  was  1  X  10'^  and  the  solution  time  step  was  1.0  seconds.  Decreasing  the  tolerance, 
decreasing  the  time  step,  or  using  the  Classical  Explicit  Procedure  or  Classical  Implicit 
Procedure  produced  the  same  results.  Two  distinct  domain  models  were  prepared:  one  in  which 
thermocouple  pads  were  absent  and  one  in  which  thermocouple  pads  were  present.  The  width  of 
the  thermocouple  pads  was  assumed  to  be  0.15  m,  as  prescribed  by  CAN/ULC-S101-M89 
[1989].  A  total  of  714  cells  were  used  to  model  the  configuration  lacking  the  thermocouple  pads 
and  1,029  cells  were  used  to  model  the  configuration  that  included  the  thermocouple  pads. 
Increasing  the  number  of  cells  to  2,668  and  3,868  for  the  two  configurations  did  not  significantly 
change  the  temperature  predictions. 

The  model  was  validated  by  comparing  the  predicted  and  measured  average  temperature  of 
the  gypsum  at  various  locations  within  the  wall  assembly,  shown  in  Figure  D-36.  The  tranient 
temperatures  at  these  locations  are  shown  in  Figure  D-38.  Figure  D-39  and  Figure  D-40  show 
the  temperature  profile  within  the  wall  assemblies  after  45  minutes  exposure  to  the  CAN/ULC- 
S101-M89  [1989]  temperature  profile. 
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The  fire  resistance  rating,  based  on  an  average  temperature  increase  of  140°C  or  a  peak 
temperature  increase  of  1 80°C  on  the  unexposed  surface  when  measured  using  thermocouple 
pads,  was  also  used  as  a  comparison  basis.  In  both  cases,  the  temperature  measurement  includes 
a  thermocouple  pad.  The  results  are  as  follows  [Sultan,  1996]: 

•  Predicted  time  average  temperature  rise  exceeds  140°C  above  ambient:  48.5 

•  Measured  time  average  temperature  rise  exceeds  140°C  above  ambient:  54.0 

•  Predicted  time  peak  temperature  rise  exceeds  180°C  above  ambient:  49.3 

•  Measured  time  peak  temperature  rise  exceeds  180°C  above  ambient:  52.0 

The  predicted  fire  resistance  rating  of  the  assembly  is  thus  48.5  minutes  and  the  measured 
fire  resistance  rating  was  52  minutes. 


Exposed  gypsum  panel,  inside  cavity  (measured) 

Une)q}osed  gypsum  pan^,  inside  cavity  (measured) 

Unexposed  gypsum  panel,  external  side,  thermocouple  under  pad  (measured) 
Unexposed  gypsum  panel,  external  side,  bare  thermocouple  (measured) 
Exposed  gypsum  panel,  inside  cavity  (predicted) 

Unexposed  gypsum  panel,  inside  cavity  (predicted) 

Unexposed  gypsum  panel,  external  side,  thermocouple  under  pad  (predicted) 
Unexposed  gypsum  panel,  external  side,  bare  thermocouple  (predicted) 


Fig.  D-38  —  Predicted  versus  Measured  Gypsum  Surface  Temperatures  for  Wall  Assembly  - 
714  and  1,029  Cell  Models  [Sultan,  1996;  Alfawakhiri,  2000]. 
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Fig.  D-39  —  Temperature  Profile  in  Wall  Assembly  45  Minutes  after  Start  of  CAN/ULC-S101-M89  [1989] 
Exposure  -  714  Cell  Model,  No  Thermocouple  Pads  (Note:  x  an  y  Axes  Use  Different  Scales) 
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Fig.  D-40  —  Temperature  Profile  in  Wall  Assembly  45  Minutes  after  Start  of  CAN/ULC-S101-M89  [1989] 
Exposure  -  1,029  Cell  Model,  Thermocouple  Pads  Included  (Note;  x  my  Axes  Use  Different  Scales) 


Sensitivity  of  Results  to  Input  Parameter  Assumptions 

The  sensitivity  of  the  calculation  results  to  uncertainties  in  the  material  properties  and  the 
boundary  conditions  is  assessed  relative  to  predicted  transient  average  temperature  at  specific 
locations  and  to  the  predicted  fire  resistance  rating  of  the  wall  assembly.  The  material  properties 
of  the  gypsum  were  varied  by  plus  or  minus  25  percent.  Similarly,  the  radiation  and  convection 
boundary  conditions  on  all  surfaces,  including  the  cavity  boundaries,  were  increased  and 
decreased  in  order  to  determine  the  change  in  the  predicted  fire  resistance  rating.  Figures  D-41 
and  42  show  the  sensitivity  of  the  transient  average  temperature  at  Location  #2  to  uncertainty  in 
the  material  properties  and  boundary  conditions,  respectively.  Likewise,  Figures  D-43  and  44 
show  the  sensitivity  of  the  temperature  predictions  at  Location  #3  (bare  thermocouple).  Similar 
results  are  observed  for  Location  #1  and  Location  #3  (with  thermocouple  pads).  Table  D-10 
summarizes  the  sensitivity  of  the  predicted  fire  resistance  rating  to  parameter  uncertainty. 
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•  Measured  temperature 

-  Predicted  temperature  (baseline) 

.  Thermal  conductivity  of  gypsum  increased  by  25% 

- Thermal  conductivity  of  gypsum  decreased  by  25% 

-  Heat  capacity  of  gypsum  increased  by  25% 

- Heat  capacity  of  gypsum  decreased  by  25% 


Fig.  D-41  —  Sensitivity  of  Location  #2  Temperature  Predictions  to  Material  Property  Uncertainty 
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Temperature  (®C) 


•  Measured  temperature 

-  Predicted  temperature  (baseline) 

.  Effective  emissivity  increased  by  -25%  on  gypsum  surfaces 

- Effective  emissivity  decreased  by  -25%  on  gypsum  surfaces 

-  Convection  increased  by  25% 

-  Convection  decreased  by  25% 


Fig.  D-42  —  Sensitivity  of  Location  #2  Temperature  Predictions  to  Boundary  Condition  Parameter  Uncertainty 
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•  Measured  temperature 

-  Predicted  temperature  (baseline) 

.  Thermal  conductivity  of  gypsum  increased  by  25% 

- Thermal  conductivity  of  gypsum  decreased  by  25% 

-  Heat  capacity  of  gypsum  increased  by  25% 

- Heat  capacity  of  gypsum  decreased  by  25% 


Fig.  D-43  —  Sensitivity  of  Location  #3  (Bare  Thermocouple)  Temperature  Predictions  to  Material  Property  Uncertainty 
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•  Measured  temperature 

-  Predicted  temperature  (baseline) 

.  Effective  emissivity  increased  by  -25%  on  gypsum  surfaces 

- Effective  emissivity  decreased  by  -25%  on  gypsum  surfaces 

-  Convection  increased  by  25% 

-  Convection  decreased  by  25% 


Fig,  D-44  —  Sensitivity  of  Location  #3  (Bare  Thermocouple)  Temperature  Predictions  to 
Boundary  Condition  Parameter  Uncertainty 


Table  D-10  indicates  that  the  calculated  fire  resistance  rating  for  gypsum  wall  assembly  is 
somewhat  sensitive  to  uncertainty  in  the  material  properties  of  the  gypsum  material.  A  twenty 
five  percent  increase  or  decrease  in  the  material  property  values  could  result  in  a  10  -  25  percent 
increase  or  decrease  in  the  predicted  fire  resistance  rating  using  the  CAN/ULC-S101-M89  [1989] 
criteria  for  wall  assemblies. 
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Table  D-10  —  Sensitivity  of  the  Calculated  Fire  Resistance  of  a  Gypsum  Wall  Assembly  to  Uncertainty 
in  the  Material  Property  Data  and  Boundary  Condition  Parameters 


Modified  Parameter 

1  Average  temperature  of  140°C 

1  Peak  Temperature  of  1 80^C  | 

Time  (min) 

Time  (min) 

Predicted  (unmodified  baseline) 

48.5 

- 

49.3 

- 

Gj^sum  conductivity  increased  by 
25% 

42.0 

42.5 

-6.8 

(- 13.8%) 

Gypsum  conductivity  decreased  by 
25% 

59.1 

+  10.6 
(+21.9%) 

60.2 

+  10.9 
(+22.1%) 

Gypsum  heat  capacity  increased  by 
25% 

57.2 

+  8.7 
(+  17.9%) 

57.8 

+  8.5 
(+  17.2%) 

Gypsum  heat  capacity  increased  by 
25% 

39.5 

-9.0 

(-  18.6%) 

40.5 

-8.8 

(+  17.8%) 

Effective  emissivity  increased  by 
~25%  on  all  gypsum  siufaces 

47.8 

49.5 

mm 

Effective  emissivity  decreased  by 
~25%  on  all  gypsum  surfaces 

49.8 

50.8 

fsgBmmm 

Convection  increased  by  25% 

48.5 

1  X 

49.0 

Convection  decreased  by  25% 

48.7 

49.7 

‘Note  that  an  increase  in  the  predicted  fire  resistance  rating  is  observed  when  the  effective  emissivity  is  increased 
and  decreased.  The  reason  this  occurs  is  because  the  emissivity  of  the  thermocouple  pad  remains  constant. 


4.0  CONCLUSIONS 


HEATING  has  compared  successfully  with  eight  validation  cases  of  increasing  complexity. 
Validation  cases  were  selected  based  on  the  availability  of  test  data  and  material  properties  for 
the  configurations  tested.  The  validation  basis  varied  from  case  to  case  depending  on  the 
particular  test  available,  and  included  fire  resistance  ratings,  temperatures  at  fixed  locations  at 
fixed  times,  and  transient  temperatures  at  fixed  locations.  Convergence  was  demonstrated  in 
each  case  by  increasing  the  number  of  nodes,  reducing  the  time  step,  reducing  the  tolerance,  and 
using  alternative  solution  methods. 

The  sensitivity  of  each  validation  case  to  the  input  parameters  was  also  investigated  because 
of  the  inherent  uncertainty  in  the  specific  material  properties,  especially  at  elevated  temperatures, 
as  well  as  the  boundary  condition  parameters.  An  uncertainty  of  plus  or  minus  25  percent  was 
assumed  where  applicable.  It  was  found  that  the  validation  cases  generally  separated  into  three 
categories:  those  sensitive  to  boundary  condition  parameters,  namely  the  effective  emissivity  of 
the  exposed  surface;  those  sensitive  to  the  material  properties,  in  particular  the  thermal 
conductivity  of  the  dominant  material,  and  those  that  were  insensitive  to  uncertainty  in  either 
parameter  group. 
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6.0  NOMENCLATURE 


Roman 

c 

f 

F 

h 

k 

n 

9 

t 

X 

X 

T 

w 

y 


Thermal  heat  capacity  (J/kg-°C  or  J/kg-m-K) 

A  generic  function  (units  vary) 

Radiation  configuration  factor 

Convection  coefficient  (W/m-°C  or  W/m-K),  length  of  inner  web  (m) 
Thermal  conductivity  (W/m-°C  or  W/m-K) 

Normal  coordinate  (m) 

Heat,  energy  (J  or  kJ) 

Time  (s,  min,  or  hr) 

The  first  Cartesian  coordinate  (m) 

The  Cartesian  position  vector,  equal  to  x,y  ,z  (m,m,m) 

Temperature  (°C  or  K) 
width  of  inner  flange  (m) 

The  second  Cartesian  coordinate  (m) 
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Greek 


S  Thickness  (m) 

s  Emissivity  of  a  surface 

p  Density  (kg/m  ) 

a  Stefan-Boltzmann  constant  (5.669E-8  W/m^-K'*) 

Superscript 

he  Denotes  natural  convection  exponent 

n  Natural  convection  exponent  (same  as  he) 

•  Quantity  per  unit  time  (s'^) 

"  Quantity  per  unit  area  (m'  ) 

' ' '  Quantity  per  unit  volume  (m'^) 


Subscript 

1,  2,  etc. 

Denotes  distinct  location  or  quantity 

A 

Denotes  ambient  (exposed  or  unexposed) 

B 

Denotes  boundary 

c 

Denotes  convection 

C 

Denotes  concrete  surface 

E 

Denotes  effective  or  resultant  quantity 

f 

Denotes  forced  convection,  furnace 

I 

Denotes  insulation  surface 

n 

Denotes  natural  convection 

r 

Denotes  thermal  radiation 

s 

Denotes  surface 

S 

Denotes  steel  surface,  general  surface 

00 

Denotes  an  ambient  quantity 
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